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BASIC PROPERTIES AND VARIABILITY 


Frangois /?. Querci 


INTRODUCTION TO GIANT AND 

SUPERGIANT M, S, AND C STARS 

The most diverse of any spectral class, M 
stars range from tiny objects barely able to burn 
core H — and some which probably cannot — 
to dust-enshrouded supergiants. In mass, they 
range from a few hundredths of a solar mass 
to several tens of solar masses and, in radius, 
from hundreds of thousands of km to several 
AU. Spanning the evolutionary spectrum, M 
stars range from pre-main-sequence dwarfs to 
stars ready to end their asymptotic giant branch 
(AGB) stage with a transition to a supernova 
or planetary nebula. Atmospheres of M stars 
likewise include a wide range of conditions and 
phenomena described in this and subsequent 
chapters. As unevolved low-mass stars, main- 
sequence M stars form a class quite distinct 
from the red-giant stars that occupy most of 
our attention, and these dwarfs are treated ex- 
clusively in Chapters 9 and 10. 

Because of the unusual variety of structures 
and processes that appear near the end of nu- 
clear burning and the additional complications 
of binarity, it is not easy to find common char- 
acteristics of the M giant and supergiant stars 
and their numerous relatives of peculiar chem- 
ical composition. The authors prefer to reserve 


the term “cool” for these objects, even though 
that term is also sometimes applied to the F, 
G, and K stars, which are so different that we 
would call them by comparison “warm” or 
“intermediate” stars. Among these cool giant 
and supergiant stars — distinguished firstly by 
spectroscopy into the classes M, MS, S, SC, C, 
R, N-type, CH, Ba, mild (or marginal) Ba, J 
stars, and more — are objects with a wide range 
in mass, composition, size, and evolutionary 
history. Collectively they are called red giants. 
They are red because they are cool (among the 
coolest stars known); their photospheres range 
in effective temperature from perhaps 5000 K 
for the warmest R stars to less than 2000 K for 
the coolest Mira-type variables, and some dust- 
enshrouded stars may be cooler still. Many 
reasons for their becoming giants have been 
elaborated in the scientific literature, but most 
of these have been shown to be wrong, and the 
discussion of that phenomenon continues (cf. 
Eggleton and Faulkner, 1981; Renzini, 1984). 

Red-giant stars are of considerable interest 
for several reasons. Their unusual evolutionary 
states provide both a unique means to test, and 
a unique reason to generalize, stellar evolu- 
tionary theory (cf. Scalo, 1981; Iben and Ren- 
zini, 1983; Wood, 1985; Iben, 1984, 1985). 
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Their chemical composition is important to our 
understanding of nucleosynthesis and mixing 
and, coupled with mass loss, to the enrichment 
of the interstellar medium. Their outer atmo- 
spheres contain unusual regimes of temperature 
and pressure that provide a series of nice 
challenges to theories of stellar atmospheres, 
both static and dynamic. Eventually, such pro- 
cesses as grain formation, chromospheric heat- 
ing, circumstellar chemical reactions, and mass 
loss must be quantitatively understood; the cur- 
rent state of affairs in these areas is reviewed 
in Chapters 2, 3, 4, 5, 6, and 8. 

Ultimately, one would hope to be able to 
point to a red-giant star— Betelgeuse or x Cyg 
or IRC + 10216— and describe in detail its 
characteristics and evolutionary history. Al- 
though this desirable goal is presently still out 
of reach, some broad outlines and details of 
specific processes and phases are available. 
Sometime after hydrogen core exhaustion, al- 
most all stars, except the most massive, become 
red giants or supergiants. Although the exact 
cause of this transformation is still a matter of 
some dispute, there is no doubt that such stars 
develop an inflated atmosphere (including en- 
velope) around a compact core. 

Most or all of the red giants of interest in 
this book are intermediate-mass (1 = M = 9 
Mq ) stars of Population I, and these are also 
the best studied. As core hydrogen is exhausted 
and the main-sequence stage ends, a compact 
helium core develops, and stars move rapidly 
to the right in the HR diagram (toward lower 
temperatures at approximately constant lumi- 
nosity) and then upward on the red-giant 
branch (RGB). For the first time, some mixing 
to the surface of nuclearly processed material 
occurs in all stars (the “first dredge-up”). 
Rather specific predictions of the results of this 
first mixing episode on the observed surface 
abundance are available (Iben and Renzini, 
1983; Iben, 1984), and observations confirm the 
general correctness of these predictions 
(Lambert, 1981). 

At some epoch, helium ignites in the core, 
either by a “helium core flash” in the degen- 
erate cores of lower mass stars ( M less than 2.3 


M o ) or quietly in the nondegenerate core of 
more massive objects. Additional mixing (sec- 
ond dredge-up) may occur upon core helium ig- 
nition in more massive (4 to 9 M q ) stars. In 
any case, core helium ignition lifts the degen- 
eracy, if any, and the star moves rapidly left- 
ward on the HR diagram to a “helium burn- 
ing main sequence.” After core helium is 
burned and a compact, electron degenerate core 
of carbon and oxygen is formed, helium burn- 
ing continues in a shell, and the star again 
moves coolward and then sharply upward 
toward higher luminosities — the AGB phase of 
a star’s life. It is now powered by double-shell 
burning— an outer shell of hydrogen burning 
and an inner shell of helium burning. As the 
helium shell burns outward, it narrows and 
flickers (pulses), driving intermittent convection 
and leading to a third dredge-up. Cool S- and 
N-type carbon stars are presumably formed 
during this AGB phase (Iben and Renzini, 1983; 
Iben, 1984; Wood, 1985), which is also marked 
by dredge-up of heavy metals (s-process ele- 
ments) in many stars. Inferring the details of 
these remarkable events by careful studies of 
elemental and isotopic abundances for all types 
of peculiar red-giant stars is presently an area 
of intense interest (Lambert, 1981, 1985;Tsuji, 
1985b). 

Of more direct connection to this series, 
mass loss becomes general in AGB stars (see 
Goldberg, this volume; Dupree, 1986). Finally, 
part or all of the AGB stars develop a super- 
wind (Iben and Renzini, 1983) that leads to 
rapid mass loss, and the star becomes a plane- 
tary nebula. Exactly which stars become plan- 
etaries is, however, presently unclear. After the 
ejection of a planetary nebula, the remnant 
compact core is a white dwarf that simply con- 
tinues to cool (forever), although some energy 
is released by hydrogen burning in the base of 
the atmosphere. A small fraction (25 percent) 
of white dwarfs undergo a final helium pulse 
and briefly return to the red-giant region of the 
HR diagram, perhaps as RCB stars (Iben, 1984, 
1985). 

Certain red-giant stars do not fit the above 
simple scheme for single intermediate-mass 


Population I stars. The most notable mysteries 
are the barium (or Ba II) stars, so-called be- 
cause of the enhanced line of barium at 4554 A 
in their spectrum, for their low luminosity in- 
dicates an insufficiently advanced evolutionary 
stage to have mixed to the surface the observed 
excess of carbon and s-process elements. The 
discovery of their binary nature (McClure, 
1983) neatly solved the problem; the barium 
stars were contaminated with mass thrown on- 
to the star by its present white-dwarf compan- 
ion at the time that companion was itself an 
AGB star. But at least some of the white dwarf 
companions have cooling times longer than the 
red-giant lifetime of these stars, and the situa- 
tion is presently unclear. R stars pose a similar 
but tougher problem (Iben, 1984 and discussion 
by Wing). As with the barium stars, R stars are 
insufficiently luminous to have dredged up suf- 
ficient carbon to become carbon rich (C/O > 
1). Yet these show no more than average binari- 
ty (McClure, 1985), which would appear to 
eliminate that otherwise attractive hypothesis. 
Certain S and MS stars also may prove to be 
challenges to the present theory (Iben, 1985). 

Important questions about red giants are: 
What is the source of neutrons for heavy ele- 
ment formation? What are the various correla- 
tions between these physical properties and 
others such as their galactic distribution, their 
spatial velocity, their age, their periods of 
variability, and the shape of their light curve? 
Is their mass loss steady, periodic, or sporadic? 
What mechanisms are responsible? How can 
mass loss be measured? What are the main 
criticisms of such measurements? How to pro- 
gress about the evolution of these various types 
of stars that are concentrated on the right side 
of the top of the HR diagram? What fun- 
damental physical data are necessary for im- 
proving the model of evolution? Are these data 
sufficiently accurate? What confidence can we 
have in such models? What is the place of these 
stars in present-day astrophysics? 

Clearly, answers to all of these questions are 
not yet available, and of the remainder, we 
must carefully select our material. In this chap- 
ter, after a general overview, we stress the vari- 


ability because virtually all the stars that are 
discussed are variable, and all the more since 
the variability is due to nonthermal processes, 
the subject of this series of Monographs. This 
chapter will mainly discuss photometric vari- 
ability ; spectroscopic variability is presented in 
the following chapter by M. Querci. 

Particular questions that are addressed on 
variability include the following: What are the 
various characteristic times of the observed 
variations? Are they correlated with any other 
characteristic period intrinsic to the given ob- 
ject? To what radial or nonradial motions of 
matter do such variations correspond? How are 
the irreversible changes explained? Are the lat- 
ter a sure sign of the stellar evolution of these 
stars which suffer them? 

What should we recommend to people who 
want to help to disentangle these fundamental 
questions? In the present exploratory state of 
our knowledge of the red giants and supergiants 
and because the observed variations are so com- 
plex, every feature found to be variable must 
be noticed and analyzed, and all correlations 
must be tested in an attempt to list better the 
various motions of matter suffered by the dif- 
ferent stellar layers. Simultaneously, we must 
endeavor to find the pulsational frequencies 
with methods adapted to the necessarily inter- 
mittent ground-based observations. Also, the- 
oretical studies of the periodic phenomena of 
Miras, which seem to be replaced by chaotic or 
aperiodic phenomena in the semiregular and ir- 
regular variables, are to be pursued. To analyze 
the variability both observationally by monitor- 
ing some characteristic stars and theoretically 
by developing red-giant hydrodynamics Gust 
beginning) would help to clarify the evolu- 
tionary problems raised in this part of the HR 
diagram, where many apparently diverse ob- 
jects are found. 

Basic properties of the M, S, and C vari- 
ables— more specifically, abundance classifica- 
tion, luminosity, mass, radius, surface gravi- 
ty, extension of the outer layers, effective 
temperature, abundances, space distribution, 
period or pseudoperiod, and shape of the light 
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curves and their variations — have been de- 
scribed in the past by many reviews: Bidelman 
(1956), Blanco (1965), Hack (1967), Beeckmans 
(1969), Spinrad and Wing (1969), Fujita (1970 
and references therein), Vardya (1970), Alksne 
and Ikaunieks (1971, updated by Baumert in 
1981), Baumert (1972), Rybski (1973), Waller- 
stein (1973), Barbaro and Dallaporta (1974), 
Boyarchuk (1974), Ikaunieks (1975), Keenan 
and McNeil (1976), Johnson (1978), Keenan 
and Boeshaar (1980), de Jager (1980), Fujita 
(1980), Feast (1981), Scalo (1981), Duerbeck 
and Seitter (1982), Alksne et al. (1983), 
McClure (1984), Jaschek (1985), and Wing 
(1985). Even this list is far from being ex- 
haustive. Some reviews are especially devoted 
to C stars, and others are restricted to M stars, 
each author favoring a particular point. A few 
recent colloquia or workshops are completely 
devoted to such stars: “Physical Processes in 
Red Giants,” Erice in 1980; “Cool Stars with 
Excesses of Heavy Elements,” Strasbourg in 
1984; and “Mass Loss in Red Giants,” Los 
Angeles in 1984. 

In the following, we present the various 
classifications of these red-giant stars, their 
photometry and colors, their space motion and 
space distribution, and their intrinsic proper- 
ties. I must apologize for showing my own 
biases and for omitting many important studies 
and reviews. 

Classifications 

Much more varied than warmer stars, the 
M, S, and C stars can be classified in several 
ways. Classification by spectra is the most ob- 
vious and most common scheme. For the varied 
objects found among the red giants, classifica- 
tion by chemical composition or key abundance 
ratios is important. These objects can also be 
classified by photometry and color. Finally, 
since almost all red giants are variable in light, 
they are classified by type of variability. All of 
these must be considered to form a complete 
picture of a red-giant star. 

Even before a clear understanding of their 
nature was available, red-giant stars were classi- 


fied on the basis of their spectra. The normal 
red giants — the M stars — were classified ac- 
cording to the increasing strength of the TiO 
bands. In contrast, other stars have spectra 
dominated by such carbon-rich molecules as 
CN and C r Quite early, it was realized that 
the tightly bound molecule CO controls the 
spectra because CO forms in cool-star atmo- 
spheres until the less abundant of either carbon 
or oxygen (which are in turn the most abundant 
elements after hydrogen and helium) is used up. 
Whichever of the two elements is more abun- 
dant is then left over to form other molecules, 
which dominate the spectrum. The ratio of 
these two— the C/O ratio — thus becomes a cru- 
cial determinant of red-giant spectra. Other 
abundance parameters such as metallicity and 
enhancement of s-process elements (those heavy 
elements— Sr, Y, Zr, and Ba— produced pri- 
marily by slow neutron capture and beta decay) 
are also key classification criteria easily deter- 
mined from spectra. 

A didactic presentation of the M, S, and C 
stars and related objects and their classification 
through abundance ratios is given in Jaschek 
(1985). The relative position of the different 
groups of stars can be easily visualized in the 
tridimensional representation (Figure 1-1) based 
on the C/O ratio, the metal abundances, and 
the s-process abundances. The C/O abundance 
ratio is the primary criterion differentiating the 
oxygen-rich M (C/O < 1) and the S stars (C/O 
» 0.8 to 1) from the carbon-rich stars (C/O 
> 1). The second gives the cool stars that are 
metal-deficient, such as the CH stars or the 
HdC stars, and the third distinguishes the 
degree of enhancement in the abundance of s- 
process elements. 

Before giving a description of the main types 
of stars mentioned in Figure 1-1, let us recall 
that these stars are located between the Hayashi 
lines for fully convective stars of various solar 
masses (see Figure 40 in de Jager, 1980, adapted 
from published HR diagrams: e.g., from Scalo, 
1976). Wood (1985) found an evolutionary se- 
quence up to the AGB in the M — S — C, but 
the quantitative agreement between the theory 
of shell-flash mixing and observations is still 
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Figure 1-1. Main varieties of red giants, according to the three-dimensional representation described 
in the text (from Jaschek, 1985). 


not satisfactory. Additional insight into the ob- 
servational constraints on theories of mixing 
and nucleosynthesis during the advanced red- 
giant stages of evolution can be gleaned from 
several recent reviews (Scalo, 1981; Iben and 
Renzini, 1983; Iben, 1984, 1985). Although we 
have chosen not to discuss in detail the red 
giants in other galaxies, conclusions drawn 
from these stars in the Magellanic Clouds are 
illuminating, particularly with respect to their 
evolution, and a later section provides an 
overview. 

The M stars , also called oxygen stars in 
parallel with carbon stars because the composi- 
tion is basically oxygen-rich (solar C/O * 0.6), 
are characterized by molecular absorption 
bands of such oxygen-rich molecules as TiO 
and VO in the visible. These bands are very 
weak in MO stars; their strength increases with 


later (cooler) types. In fact, M stars are classi- 
fied from MO to M8 primarily by the increasing 
strength of the bands of the alpha system of 
TiO from 4350 to 4950 A (e.g., Keenan and 
McNeil, 1976; Wing, 1979; Keenan and Pitts, 
1980). A few stars have even been classified M9 
or M10. In the infrared, M stars are character- 
ized by CO, OH, SiO, and H 2 0 absorption 
bands. The 10-/*m region shows the emission 
of silicate grains, which are discussed in Chap- 
ters 2, 4, and 5. Maser lines of CO, OH, H,0, 
and SiO are observed in the coolest M Mira 
stars. 

Among the M variables, some belong to the 
regularly pulsating group of cool stars, the 
Miras (see the section on Variability Types). 
Although the Miras change their spectral sub- 
type and magnitude during their light cycle, 
they occupy a rather distinct zone around the 
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giant and subgiant tracks. However, nonvari- 
able M giants (see the section Nonvariable 
Stars), such as /3 And (MO III), also fall in this 
region of the HR diagram. The OH/IR stars 
(i.e., very red objects with double-peaked 1612- 
MHz OH maser emission) are a continuation 
of M Mira stars to lower temperatures, larger 
masses and luminosities, and longer periods of 
pulsation (Duerbeck and Seitter, 1982). Their 
evolution is discussed by de Jong (1983); the 
most luminous ones (8 to 13 M 0 ) are to be 
identified with the core helium-burning super- 
giant stars, while those of 3 to 8 M Q are AGB 
stars. Some of them are OH objects with unde- 
tected IR counterparts (Jones et al., 1982). They 
pulsate only weakly or not at all, and their 
mass-loss rate is as large as 10" 4 M Q /yr (Baud 
and Habing, 1983); they possibly reach the end 
of their evolution on the AGB (Herman, 1983). 
This is consistent with the very red colors ob- 
served by IRAS (Olnon et al., 1984). If our cur- 
rent ideas of stellar evolution are right, these 
stars soon produce planetary nebulae (Olnon 
et al., 1984). 

The C stars (carbon stars) are distinguished 
by the presence of C 2 Swan bands in their 
spectra (Keenan and Morgan, 1941; Keenan 
and McNeil, 1976; Fujita, 1980). They are fur- 
ther divided into R-type stars (the hotter ones) 
and N-type stars (the cooler ones), the criterion 
being the strength of the violet flux, which is 
very weak in N-type stars (Shane, 1928) — 
weaker than in corresponding M stars and 
weaker than blackbodies at the effective 
temperature. Historically, the R and N types 
were initially used in the HD catalog. The first 
to observe N stars was Secchi (1868). Later, 
Pickering (1898) identified a small group of C 
stars that differed from the N stars in display- 
ing sensible flux short ward of 4500 A, and 
Espin (1889, 1898) remarked that none of the 
stars observed by Secchi were from this group. 
In 1908, Pickering defined the R subtype (i.e., 
the present R type) by the extension of the con- 
tinuum shortward of 4700 A, the earlier stars 
having the larger extension. This deficiency of 
flux in N-type stars continues into the 
ultraviolet. Finally, the R and N types were 


unified by the Cx,y classification of Keenan and 
Morgan (1941). Let us note that, in the Keenan- 
Morgan classification, (Cx,y), the first figure 
(0-9) following the “C” is a temperature in- 
dex obtained from atomic lines; the second 
figure (0-9) is a carbonicity index obtained 
from the strength of the C 2 bands. Underlying 
the Cx,y classification was the assumption that 
x and y (temperature and carbonicity) were in- 
dependent. If so, the Cx classification would 
have arrayed the carbon stars in a temperature 
sequence parallel to that of M stars. Nature is 
not so simple, however, and the ambiguity of 
this classification has been amply demonstrated 
(e.g., Fujita, 1980) both observationally from 
variability studies (Eggen, 1972b), infrared 
broadband colors (Mendoza and Johnson, 
1965), near-infrared narrowband colors 
(Baumert, 1972), and infrared spectrophotome- 
try (Goebel et al., 1978) and theoretically from 
calculations of the effect of carbon enhance- 
ment on opacities (Scalo, 1973) and atmo- 
spheric structure (Johnson, 1982) and on effec- 
tive temperatures (Tsuji, 1981c). 

The early R stars (R0 to R4) are termed “or- 
dinary R stars” because their hydrogen lines 
have strengths similar to those in normal G9 
to K2 giants of the same effective temperature. 
Vandervort (1958) gives the criteria for the 
classification of the R stars between R0 and R8. 

In the carbon stars, the visible spectrum con- 
tains bands of the carbon-rich molecules CN, 
C 2 , C y SiC 2 , and CH instead of the oxide 
bands that characterize the M stars spectra, and 
these bands dominate the spectra of cooler 
stars— the late R (R5 to R9) and the N stars. 
Apart from the CO bands seen in all cool stars 
(though unexpectedly weak in the R0 to R4 
stars), the infrared (1 to 5 ;tm) is increasingly 
chopped up by bands of CN, C 2 , HCN, and 
C 2 H 2 toward later spectral types. Emission at 
1 1 #tm due to graphite and SiC grains is often 
seen. In the radio frequencies, CO, CN, SiO, 
CS, C 2 H, HNC, HCN, C 3 N, and complicated 
molecules such as the cyanopolyyne family 
(HC, ,N) or CH.CN are detected in the ex- 
treme dust-enshrouded C stars. 
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The s-process elements are stronger in late 
carbon stars than in M stars (Utsumi, 1970), 
but the l3 C-rich carbon stars (J stars) appar- 
ently show little or no enhancements of s- 
process elements with respect to Fe (Utsumi, 
1985). 

The spectra of N-type stars show a stronger 
absorption in the blue and the ultraviolet than 
do the R-type spectra. This classic problem, 
termed the “violet opacity,” is discussed fur- 
ther under Photometric Observations. It seems 
that this degree of faintness in the ultraviolet 
region of N-type spectra is not entirely due to 
their lower temperature (de Jager, 1980). The 
absorption by molecular bands (particularly 
C 3 ) and graphite and silicon carbide particles 
could be an important contributor to this 
observed faintness. Stars redder than B-V = 
3 are nearly always carbon stars, of which the 
reddest are the cool N types with B-V » 5 or 
6 (Wing, 1985); for comparison, the solar B- V 
value is 0.62, and MO to M6 giants have B-V 
* 1.55 (Johnson, 1966). 

The early R stars, extensively described by 
Dominy (1982, 1984), have an abnormal oxy- 
gen/carbon ratio relative to the M stars by 0.7 
dex; the nitrogen and oxygen are enhanced by 
only 0.2 dex. The iron abundance is near the 
solar value, except in HD 100764, which has 
an iron abundance 0.6 dex lower than the solar 
one. The s-process elements such as Y, Zr, Mo, 
Ba, La, Ce, Nd, and Sm are not enhanced in 
these hotter carbon stars. A bright member of 
this class, HD 156074, was well analyzed by 
several authors: Greene et al. (1973) for the Fe 
abundance, Wyller (1966), Gordon (1967), and 
Fujita and Tsuji (1976) for the C l2 /C 13 ratio, 
and Yorka (1981) for the N abundance. The ele- 
mental abundances in the early R stars are in- 
terpreted as products of mixing during the 
helium flash (Dominy, 1982; Wood, 1985). On 
the other hand, due to uncertainties both in 
observations (mainly in luminosity) and in evo- 
lution theory, Tsuji (1981c) remarks that it is 
difficult to decide whether N stars (N-type 
Miras are not considered) are produced just 
after the onset of He-shell flash. The carbon- 
star evolutionary “mystery,” paraphrasing 


Iben (1981b), has been studied specifically by 
this author (Iben, 1984) and by Lucy et al. 
(1986). 

Beside these “classical” carbon stars, a class 
of stars that are carbon-rich but have the addi- 
tional peculiarity of hydrogen deficiency — the 
so-called hydrogen-deficient carbon (HdC) 
stars — group together, following the convenient 
classification of Richer (1975): (a) the R Coro- 
nae Borealis (RCB) stars; (b) the helium (He) 
stars; and (c) a few supergiant R stars (e.g., HD 
182040) (see the section Intrinsic Properties) 
that appear to be nonvariable (McKellar and 
Buscombe, 1948; Bidelman, 1953; Warner, 
1963, 1967) and are called the nonvariable HdC 
stars. 

The group of HdC stars is deficient in hy- 
drogen by a factor up to 10 5 or more com- 
pared to the H/Fe ratio in the Sun. The domi- 
nant constituent of their atmosphere is helium, 
although it is not spectroscopically dominant 
in all the HdC stars. Because no HdC star is 
known to be a close binary, the observed abun- 
dances can only be due to the ejection of the 
outer layers during the natural course of evolu- 
tion. The lifetime of the HdC stage appears to 
be - 10 3 years; if so, secular changes should 
be observable. It is proposed that these super- 
giants (see below) quickly evolve to white 
dwarfs (Iben, 1981a, 1985,). 

The RCB stars are characterized, at maxi- 
mum light, by narrow and sharp absorption 
lines of metals and very strong molecular ab- 
sorptions such as CN in R CrB and RY Sgr, 
and C 2 in S Aps and HV 5637. During 
minimum light and the subsequent increase, 
emission lines of helium (e.g., Querci and Quer- 
ci, 1978) and metals usually appear in the spec- 
trum. These stars have a surplus of carbon (by 
a factor of 3 to 10 (Warner, 1967)), as well as 
a surplus of helium, compared to the Sun. 
Earlier, Ludendorff (1906), Berman (1935), and 
Herbig (1949) remarked on these peculiar abun- 
dances in R CrB, and Bidelman (1948) showed 
these peculiarities in XX Cam. They are be- 
lieved to be due to the loss of the hydrogen 
envelope (Sackmann- Juliana et al., 1974). 
Some differences exist between them; in XX 
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Cam, carbon is a little more abundant and hy- 
drogen more deficient than in RY Sgr (Rao et 
al„ 1980) log [(C/H RYSg / (C/H) XXCam ) - 
-1.16. Schonberner (1975) has estimated this 
ratio to be -0.9. The abundance in oxygen is 
also slightly different; it is 3 times more abun- 
dant in R CrB than in RY Sgr. Most RCB stars 
are members of the R type; some of them are 
known to have earlier type spectra, including 
the prototype star itself, which is of spectral 
class cGOep. In the HdC group, they are called 
cool HdC stars. The RCB are eruptive stars 
(and not pulsating stars, such as the “classical” 
C stars; see the section Types of Variability) 
with long periods; moreover, they are the only 
eruptive stars in which the quiescent state is the 
state of maximum brightness. Their minima are 
interpreted by ejections of obscuring matter 
(soot, particles of graphite, etc.). 

Helium stars and RCB stars have several 
common properties, although the He stars have 
a higher temperature (leading to their being also 
called hot hydrogen-deficient stars, like BD 
+ 13 °3224 (Kilkenny and Lynas-Gray, 1982)), 
a larger gravity, and a smaller pulsational pe- 
riod than the RCB stars. MV Sgr appears to 
be a possible evolutionary link between these 
two classes (Hill, 1967) with its light variations 
like those of RCB stars and its spectrum like 
that of the helium star, HD 124448. Tutukov 
and Iben (1985) show theoretical evolutionary 
connections between them. Herbig (1964) noted 
that the strongest lines in the spectrum of MV 
Sgr are from He I and C II. The latter is 
somewhat weaker than He I lines, whereas in 
BD + 10°2179, the lines of the two atoms are 
of the same strength. Warner (1967) analyzed 
some nonvariable helium stars and found abun- 
dance characteristics like those of the non- 
variable HdC stars. Envelopes around the 
helium stars were suspected by Hill (1964), such 
as around HD 168476, which appears like an 
A-type star with many lines of ionized elements 
of the iron group. This star shows that the ex- 
citation temperature for the ionized elements 
is much lower than that for the other elements, 
leading Hill to suggest that they are formed in 
a low-temperature shell. The immediate pro- 


genitors of the helium stars are believed to be 
RCB stars (Schonberner, 1975). However, 
nothing is really known about their common 
origin and by what mechanism they have lost 
their initial hydrogen-rich atmosphere. 

As for the nonvariable HdC stars , the stars 
HD 137613, HD 173409, HD 175893, and HD 
182040 are hydrogen-deficient by factors of 
more than 10 5 compared to the H/Fe ratio in 
the Sun. The deficiency is around a factor of 
50 in HD 148839; this last star shows C 2 and 
C I enhancement, while CH and Balmer lines 
are absent. It represents “an intermediate type 
between normal C stars and the more extreme 
hydrogen-deficient stars” (Warner, 1967). In 
all these stars, the carbon is overabundant by 
factors of 3 to 10, but C 13 is not observed at 
all. All elements heavier than oxygen have a 
roughly solar abundance relative to iron. 

Among carbon-rich stars, we note other par- 
ticular stars such as the CH stars , which are 
G5-K5 giants with strong CH molecular bands 
(G band at 4300 A) and enhanced lines of s- 
process elements in their carbon-rich and metal- 
poor spectrum. They have high radial velocities 
and can be considered Population II carbon 
stars. Stars with a strong G band and low spa- 
tial velocity are also observed; they are called 
CH-like stars. Bidelman and McConnell (1973) 
identified a large number of stars with a very 
weak or absent G band, the so-called weak G- 
band stars ; no correlation was found between 
the CH faintness and the line strength of metal- 
lic lines. The CN-strong stars are carbon stars 
with a metal overabundance and a slight en- 
hancement of carbon (Schmitt, 1971). 

Jaschek (1985) stated that “The S-stars con- 
stitute a natural bridge between M and C stars 
with all kinds of intermediates.” They are late 
(K5-M8) giants. These stars contain metals with 
roughly the same abundances as the M stars, 
but the enhancement of the heavy elements, in- 
cluding Zr, Y, Sr, Ba, and La, is very strong 
and consequently the ZrO bands are much 
more prominent than those of TiO. Oxides of 
yttrium, lanthanum, and vanadium are affected 
by the change in C/O between M and C stars. 
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The C and O atoms are fully trapped in CO; 
this reduction in the free oxygen supply as car- 
bon is enhanced explains the presence of many 
atomic hydrides in the S star spectra. Their B - V 
colors are between -2 and -3. 

The presence of unstable isotopes such as 
99 T c in the stellar spectrum of some S and C 
stars (Little-Marenin and Little, 1979; Smith 
and Lambert, 1985), indicates that the rising to 
the surface of the interior material has occurred 
less than about 10 6 years ago, which corre- 
sponds to a few half-lives of the Tc" isotope, 
the longest lived of the Tc isotopes, under con- 
ditions representative of interiors of AGB stars 
(Iben and Renzini, 1983). 

A very few C and S stars (T Sgr, WZ Cas, 
HR 8714, WX Cyg, T Ara, etc.) present a very 
strong absorption resonance line of lithium at 
X6707 (see the surveys of Torres-Peimbert and 
Wallerstein, 1966, and Boesgaard, 1969). These 
stars are named lithium-rich stars (LRS). The 
derived abundance of lithium is up to 10 5 
times the solar value in WZ Cas. To obtain the 
Li formation or destruction processes in such 
stars, two difficulties must be solved (Scalo, 
1981): to infer the Li/Fe ratio and to infer the 
6 Li/ 7 Li ratio. The first ratio is only tentatively 
estimated (Merchant-Boesgaard, 1970; Waller- 
stein, 1977; de la Reza and Querci, 1978) be- 
cause the X6707 line is largely blended by CN 
red-system lines. The knowledge of the second 
ratio is not easy to obtain through the analysis 
of the lithium resonance atomic line because, 
besides the CN blending, the isotopic shift is 
smaller than the observed Doppler width of the 
stellar line. Some efforts have been made in the 
direction of the molecular isotopic shift by the 
development of computations of molecular 
abundances of LiH, LiOH, and LiO (Mer- 
chant, 1967; Querci and Querci, 1977) and of 
wavelength measurements in the laboratory 
(Plummer et al., 1984). Stellar observations of 
lithium hydride must be made to provide the 
desired ratio. 

Hybrid objects are also observed; for exam- 
ple, the CS and SC stars have characteristics 
of both C and S stars (Stephenson, 1967b; 


Catchpole and Feast, 1971), and the MS stars 
are intermediate between M and S types. Kee- 
nan and Boeshaar (1980) provide reference stars 
in the sequence MS through S and SC (some- 
times called D-line stars (Gordon, 1967)) to C 
stars with weak C 2 bands. 

We conclude this description of the main 
types of cool stars by the Ba II stars (also 
named Ba or barium stars) that constitute 
another family of anomalous giant stars with 
strong Ba II features (especially the X4554 line). 
They share with the early R stars (R0 to R4) 
the fact of being “warm” or “intermediate- 
temperature” stars (see the section on Intrinsic 
Properties ). Most Ba stars have strengthened 
CN absorption (Wing, 1985) relative to normal 
G and K stars, although not as strong as in car- 
bon stars. The s-process elements, especially 
those with 38 < Z < 56, are enhanced (McClure, 
1984; Smith, 1984). An estimate of the strength 
of the barium lines was added to the spectral 
type by Warner (1965). This group has no sharp 
boundaries (Jaschek, 1985). So-called semi- 
barium stars have s-process element abun- 
dances intermediate between those seen in the 
pure Ba II stars and the M stars. The stars for 
which the enhancement is not obvious are 
named “marginal barium stars.” To designate 
the strength of the Ba II characteristics, Keenan 
and Pitts (1980) give a decimal scale that is of- 
ten used. Catchpole et al. (1977) and Yamashita 
and Norimoto (1981) have reobserved the Ba 
II sample of MacConnell et al. (1972) and con- 
clude that about one-third of these stars are 
weak-line objects, such as the subgiant CH stars 
defined by Bond (1974). 

The lack of luminosity determinations is the 
main difficulty in placing the Ba stars on the 
HR diagram, but they are generally believed to 
have luminosities of normal giants. Until the 
recent discovery of their binarity (McClure et 
al., 1980; McClure, 1983), the evolutionary 
status of the barium stars was a puzzle because 
such apparently low-mass stars should not have 
been able to have mixed to the surface nuclearly 
processed material. It now appears fairly cer- 
tain that the abundance anomalies seen in the 
real Ba stars are the result of mass transfer from 
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formerly more massive companions when the 
companions were in the asymptotic giant 
branch stage of evolution (Smith, 1984; 
Lambert, 1985). Whether the same holds true 
for mild or marginal barium stars is, however, 
still an open question (McClure, 1985). Wood 
(1985) suggests that the Ba II stars (and CH 
stars) are probably the result of mass transfer 
from an N-type star to a white-dwarf compan- 
ion in a binary system. 

To understand better the large variety of 
these giants and supergiants, which lie on the 
right side of the HR diagram and cover a large 
range on absolute visual magnitude (see the sec- 
tion on Intrinsic Properties ), several other 
classifications based on criteria other than 
abundance ratios or temperatures are proposed 
(Eggen, 1972b; Feast, 1981; Feast et al., 1982). 
The importance of finding physically meaning- 
ful, significant classifications dividing the ob- 
jects into homogeneous subgroups is stressed. 
Of the various attempts to classify these stars, 
the most famous classification, based on type 
of light variation, is the General Catalogue of 
Variable Stars (GCVS) (Kukarkin et ah, 1958, 
1976), which essentially divides the late-type 
variables into Miras, semiregulars (SR), and ir- 
regulars (L) according to their light-curve am- 
plitude and profile. This point will be discussed 
in the section on Variability Types . As noted 
by Eggen (1972b), the GCVS is a poor guide 
to the type of variation of the N-type red vari- 
ables, maybe because of their extreme ( B-V ) 
colors. This author, from an extensive study of 
red giants, suggests a subdivision of the vari- 
ables on the basis of their place in the lumi- 
nosity/temperature plane (i.e., (M^ o/ , /?-/)- 
plane, adopting Johnson’s (1966) relation be- 
tween R-I and log TJ and from their space 
motions linked to their age. 

Feast (1981) sorts out the variables by a 
spectroscopic argument — the strength of the 
emission lines: (1) the variables with strong 
emission lines, called Xe stars, are Miras (i.e., 
large-amplitude variables); (2) the variables 
with weak emssion lines at some phases, called 
X(e), are semiregulars (i.e., small-amplitude 
variables); and (3) the variables that do not 


show emission lines, called X variables, are 
semiregulars or irregulars, where X denotes the 
spectral type M, S, C, CS, SC, etc. The basic 
idea is that the emission-line strength is linked 
to the extent to which the stellar layers are 
disturbed by the variability (see M. Querci, this 
volume). Feast et al. (1982) relate this division 
to mean infrared colors and define two 
reddening-free parameters b v b r The b l$ b 2 
colors suggest an evolutionary sequence 
M-M(e)-Me that is in fact a sequence of 
increasing light amplitude and decreasing tem- 
perature. The Mira variables that show OH or 
H 2 0 maser emission are at the cold end of the 
b x - b 2 plot, while the carbon stars show a 
large scatter that is also seen in the infrared 
photometric indices (Tsuji, 1985). The SC stars 
form a tight group, whereas the CS stars are 
scattered more widely. The main sources of in- 
frared opacity, different in each stellar type, 
evidently have a strong effect on the infrared 
colors and cause the difference as to the region 
occupied by each type, as demonstrated by 
Feast et al. (1982). 

The main purpose of establishing such em- 
pirical relationships is to provide a norm with 
which red variables in other galaxies (e.g., the 
Magellanic Clouds) or in other regions of our 
Galaxy than in the neighborhood of the Sun 
(e.g., in the galactic bulge) might be compared. 

Photometric Observations 

After the pioneer astronomers who were fas- 
cinated by the red stars — let us note Espin, 
Hale, Kirch, Koch, Maraldi, Pickering, Secchi, 
Schjellerup, Shane, and Sidgraves— the first 
sustained observations of visual light curves of 
red stars began in the late 1930’s at the Harvard 
College Observatory (Campbell, Hughes, 
Payne, Sterne, etc.). For more than 50 years, 
amateur astronomers (working together in such 
groups as AAVSO, AFOEV, GEOS, etc.) have 
joined their observational efforts to those of the 
professionals. More recently, the UBV system 
and its extension into the infrared have been 
extensively used for determining such physical 
parameters of the red giants and supergiants as 
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energy distribution, effective temperature, and 
infrared excesses and for showing correlations 
among them at various phases of the light 
curves (e.g., color/magnitude relation, varia- 
tions of various indices such as (B-V) y ( I-K ), 
and (8 ft m-11 fim) indices with spectral class; 
see below). Astronomers who made important 
investigations in broadband photometry in- 
clude: Bahng, Bergeat, Bessel, Bidelman, 
Blanco, Catchpole, Eggen, Feast, Johnson, 
Landolt, Leighton, Mendoza, Morgan, Neuge- 
bauer, Nicolet, Price, Smak, Stephenson, 
Walker, Westerlund, and Wood. 

Table 1-1 gives broadband magnitudes and 
colors from the ultraviolet to the infrared for 
characteristic stars of different spectral and 
variability types. For x Cyg, an S Mira, it is 
worth noting that the magnitude variations are 
much larger in the ultraviolet than in the infra- 
red region. (A specific example of variations in 
ultraviolet flux for x Cyg is shown in Figure 
1-7, and other examples in the visual are given 
in Figure 1-8.) With the advent of more sensi- 
tive infrared detectors, broadband photometry 
across the entire spectrum became possible, and 
this material has been published and discussed 
in detail. Although readers are probably already 
familiar with these observations and it is im- 
possible to review them here in any case, we 
note a few of those papers containing funda- 
mental broadband photometric data for M 
giants (Johnson, 1966; Eggen, 1967; Johnson 
and Mitchell, 1975), carbon stars (Mendoza and 
Johnson, 1965; Eggen, 1972; Walker, 1979; 
Bergeat and Lunel, 1980; Noguchi et al., 1981), 
and Mira variable stars of all types (Mendoza, 
1967; Catchpole et al., 1979). A useful review 
of photometry of all types of peculiar red giants 
is that of Wing (1985). 

A persistent problem threading its way 
through the astrophysical literature is the vio- 
let flux deficiency of carbon stars. A normal 
M6 giant has a value of ( U-V ) = 2.43 
(Johnson, 1966), while a Ori (M2 lab) has 
(U-V) = 3.96 (Table 1-1). By contrast, a 
warmer N-type carbon star such as TX Psc has 
(U-V) = 6.11 (Table 1-1), and many carbon 
stars have even larger values. This rapid drop 


of flux toward the violet was noticed long ago 
and became the basis for the distinction be- 
tween the carbon stars with sensible flux short- 
ward of 4700 A (R stars) and those with no 
short wavelength flux (N-type stars) (e.g., 
Shane, 1928). On the basis that the flux defi- 
ciency was due to an unknown violet opacity, 
the phenomenon is termed the ‘ ‘violet opaci- 
ty” problem, which is still unsolved. Among 
several suggestions for the agent responsible, 
the leading candidates are the molecule C 3 in 
the photosphere (McKellar and Richardson, 
1955; Bregman and Bregman, 1978) or SiC 
grains in a circumstellar shell (Gilra, 1973; 
Walker, 1976). Observations with IUE 
demonstrate that the violet opacity extends in- 
to the ultraviolet as well (Querci et al., 1982; 
Johnson and O’Brien, 1983). 

The widths of these broadband filters are 
too large ( > 1000 A), however, for the analysis 
of the time behavior of the various bands of 
TiO, VO, ZrO, CO, CN, C 2 , CH, and HCN 
in the visible and the near-infrared, and the ob- 
served variations are difficult to interpret quan- 
titatively in terms of energy dissipated in the 
various stellar layers submitted to acoustic/ 
shock waves. Consequently, intermediate-band 
photometries (with filter bandwidths in the 
range 200 to 500 A), initially developed for in- 
vestigations on A to G stars (i.e., with filters 
adapted to their spectral features), are used for 
the analysis of the warmer stars among M, S, 
and C giants and supergiants; such systems are: 
DDO, Stromgren, Geneva, and Vilnius sys- 
tems. 

In 1966, Wing (1967a) pointed to a spectral 
band of 30 A width near 1 .04 pan that was free 
from serious atomic and molecular blanketing 
in M, S, and C stars. A new epoch was born: 
the variations of the continuum of these late- 
type stars could be followed. Narrowband 
filters adapted to the M star molecular bands 
appeared (Wing, 1971; Lockwood and 
McMillan, 1971). White and Wing (1978) 
classified M supergiants; Baumert (1972) used 
the Wing system to measure the CN absorption 
versus the continuum of C stars; Wing and 
Stock (1973) identified carbon stars in stellar 
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Table 1-1 

Broadband Magnitudes for Characteristic Stars 
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(1) Mendoza (1968), (2) Hackwell (1972), (3) Lee (1970), (4) Mendoza and Johnson (1965), (5) Olson and Richer (1975), (6) Johnson et al. (1966), (7) Walker (1980), (8) Mendoza 
(1967). 



fields, Piccirillo (1977, 1980) investigated the 
S stars with the ZrO and LaO filters, and Yorka 
(1983) published data on R stars. However, let 
us note that some of the filters adapted to 
selected molecules are weakly contaminated by 
other molecules. 

Unfortunately, the observations on the other 
spectral intervals such as X rays and far-UV, 
far-IR, and radio frequencies are scarce. 

The fingers of one hand suffice to count the 
Einstein X-ray satellite observations of red 
giants and supergiants (Ayres et al., 1981). For 
the M supergiants, a Ori (M2 lab) and a Sco 
(Ml lb + B), values of fJf bol are respectively 
<0.03 x 10" 7 and <0.02 x 10" 7 , and for the 
bright giant, 0 Peg (Ml II— III), this ratio is <0.2 
x 10" 7 . This is consistent with the IUE obser- 
vations and demonstrates that there is no stellar 
coronal ( T = 10 6 K) radiation observed from 
such stars. In this frame, perhaps the most im- 
portant result of the observations of soft X rays 
is the discovery, with the Einstein spacecraft, 
of a rather sharp dividing line in the HR 
diagram. For giants, those stars hotter than 
about V-R = 0.80 (K0 III) emit soft X rays in 
detectable amounts, while those cooler do not 
(Ayres et al., 1981). This coronal dividing line 
(the X rays are presumed to arise from a hot 
corona) falls close to the temperature dividing 
line defined by the appearance in the ultraviolet 
spectra of giant stars warmer than approx- 
imately K0 of lines of C IV, indicating the 
presence of a plasma with T = 10 5 K. EX- 
OSAT detections are actually made on Miras 
with a hot companion (R Aqr and CH Cyg), 
individual nonvariable M giants (ft And), Miras 
(RS Aqr), and protoplanetary nebulae (V1016 
Cyg, HM Sge, and RR Tel). It was 
demonstrated by Viotti et al. (1985) that the 
weak X-ray emission at Mira maximum is not 
correlated with the Mira-type variations. Detec- 
tion of X-ray emission in any single M, S, and 
C star would certainly warrant further in- 
vestigation, since it would clearly be an excep- 
tion and so little is presently known about the 
chromospheres of these objects. 

Only disparate attempts on far-UV flux 
measurements exist. The Orbiting Astronomi- 


cal Observatory (OAO-2) UV spectrometer ob- 
served only the M supergiant, a Ori: no varia- 
tions were detected on the seven scans (Code 
and Meade, 1979). On the other hand, Code 
et al. (1980) gave UV magnitude between 1330 
and 4250 A for five M variable and nonvariable 
giants and supergiants. Unfortunately, the large 
bandwidth required for OAO-2 photometry by 
the rapidly decreasing fluxes toward the 
ultraviolet and the complex spectra of these 
stars made the flux determination ambiguous. 
With the Astronomical Netherlands Satellite 
(ANS), Wesselius et al. (1982) have observed 
eight M giants and supergiants, some of them 
many times. The star a Ori, registered six times, 
appears to be variable; however, as the point- 
ing prevision was 1 arc-minute, two or more 
stars could be present in the field of view, and 
the observed variability cannot be conclusive. 
On the basis of three observations, the authors 
conclude that HD 100029 (M0 III) is not vari- 
able. A probable hot UV-bright companion was 
found on HD 216131 (M2 III). 

In their comparison between the current UV 
photometric and spectrophotometric systems 
used for late-type stars, Kjaergaard et al. (1984) 
point out: (a) a strong nonlinearity at low flux 
level for the TD 1 2635 A band; (b) a strong color 
effect on the TD1 2740 A band, as well as on 
the OAO-2 2950 A band, due to red leaks in 
the sensitivity functions; and (c) a strong influ- 
ence of the scattered light of the late-type stars 
on OAO-2 spectrophotometry. Only the ANS 
photometry and the International Ultraviolet 
Explorer (IUE) LWR spectrophotometry are 
found to be without large systematic errors. 
Their Table 2 gives several {UV- V) indices for 
stars of various spectral types and luminosity 
classes. Unfortunately, this table gives no data 
for giants with (B-V) > 1.5 (K5) and for 
supergiants with {B-V) > 1.8; therefore, it is 
applicable only for the hotter R and Barium 
stars, keeping in mind uncertainties in the data 
as commented by the authors. Moreover, the 
presence of emission lines such as Si II XXI 808, 
1817 or the Mg II doublet may cause uncertain- 
ties in these indices. Only the three M giants 
or supergiants, HD 6860 (M0 III), a Ori (Ml 
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Ia), and HD 217906 (M2.5 II-III), are ob- 
served with the satellites TD1, O AO-2, and 
ANS; the accuracy of the measurements does 
not permit quantitative comparisons. 

A new era in research on outer atmospheres 
and chromospheres dawned with the launch of 
IUE, for it made possible fairly routine obser- 
vations of ultraviolet spectra in either a short- 
wavelength (1150 to 2050 A) or a long- 
wavelength (1800 to 3250 A) band. Most studies 
of IUE spectra have focused, naturally enough, 
on the emission lines, and these results are dis- 
cussed in Chapter 2 (observations) and Chapter 
8 (chromospheric modeling). It became clear 
from the earliest ultraviolet observations that, 
in the cooler part of the HR diagram, spectra 
of giants and supergiants were of two distinct 
types, depending on the presence or absence of 
lines of stages of high ionization. All stars show 
emission lines from such easily ionized elements 
as Mg II and Fe II (the archetypal line being 
Mg II X2800), which are representative of a 
warm chromosphere (6000 to 15000 K). Those 
stars to the left of an almost vertical line near 
K0 also show lines from highly ionized elements 
(the archetypal line being C IV XI 550 A) repre- 
sentative of a hot (~ 10 5 K) plasma; those to 
the right of this “temperature” dividing line 
show no high-ionization lines (Linsky and 
Haisch, 1979). The asymmetry of the emission 
lines (Mg II h and k) further indicates an out- 
flow of material in the cooler giants (Stencel 
and Mullan, 1980), and this “mass-loss” divid- 
ing line lies close to the temperature dividing 
line. As described above, observations with the 
Einstein spacecraft of X rays reveal a similar 
dividing line in the HR diagram between stars 
with (to the left) and without (to the right) soft 
X-ray emission (Ayres et al., 1981), which X- 
ray dividing line is also close to the temperature 
dividing line and the mass-loss dividing line. 

These observations accord with and greatly 
expand previous knowledge that violet emission 
lines of Fe II (multiplets 1, 6, and 7) are prac- 
tically ubiquitous in the spectra of M giant stars 
(Boesgaard and Boesgaard, 1976) and are even 
seen in N-type carbon stars (Bidelman and Py- 
per, 1963), but are not seen in warmer stars, 


while the high-excitation 10830 A line of He I 
is almost universally absent from single red- 
giant stars (Zirin, 1982). The common inter- 
pretation of these results, supported by many 
observations and based on a fair amount of 
modeling, is that giant stars warmer than about 
K0 are solar-like in having warm chromo- 
spheres and hot coronae, while cooler giants 
have warm chromospheres and outflowing mat- 
ter (Linsky, 1980, 1982; Brown, 1984; Linsky, 
1985). Except for binary stars, then, we expect 
the red giants of interest here, all of which, 
except perhaps for the warmest R stars, lie 
clearly to the right of the dividing lines in the 
HR diagram, generally to show Mg II in emis- 
sion in some strength and evidence of mass loss 
but no spectral lines from several-times-ionized 
elements. 

Stickland and Sanner (1981) observed the 
far-ultraviolet continua of 15 late K and M 
giant and supergiant stars (as late as M6) with 
IUE and obtained radiation temperatures at 
1850 and 1250 A. Remarkably, the radiation 
temperature at 1850 A of all stars was close to 
3400 K. They also demonstrated that this flux 
originates in the chromosphere. Several studies 
have found a general decrease in the fraction 
of ultraviolet flux in both emission lines and 
continuum, relative to bolometric flux, as one 
goes to later spectral types among the M giants. 
IUE has confirmed and delineated this trend 
down to M6 for stars with angular diameters 
from lunar occultation (Steiman-Cameron et 
al., 1985). 

Along with numerous G and K giants, sev- 
eral M giants and supergiants have been ob- 
served with IUE. Representative spectra of 
three M giants and one supergiant (a Ori) are 
shown in the valuable atlas by Wing et al. 
(1983). When line profiles of Mg II are avail- 
able, the usual procedure is to construct a tem- 
perature profile through the chromosphere. 
Even then, no attempt has yet been made to 
reproduce by any empirical or theoretical model 
the flux of other lines or the continuum. Aver- 
age properties of the chromospheres can be 
deduced from lines of Fe II and C II. A study 
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of 14 noncoronal G8-M3 giants, for example, 
leads to temperatures of 7000 to 9000 K and 
electron densities of — 10 8 cm" 3 , and these ap- 
pear to be independent of spectral type (Car- 
penter et al., 1985). Most of the effort in 
warmer (F, G, and K) stars has been directed 
toward linking the observed line emission to 
other stellar parameters such as mass, luminosi- 
ty, age, or rotational velocity. Although the M, 
S, and C stars of interest here have often been 
neglected, progress can best be judged from 
several recent reviews (Linsky, 1980, 1982, 
1985; Brown, 1984; Bohm-Vitense and Quer- 
ci, 1986). 

More attention has been given to Betelgeuse 
than any other star of interest here; it is often 
taken as the archetype red giant. Although 
spectral lines are discussed more fully in 
Chapter 2, we mention here some of the results 
so far derived. Analyses of the Mg II h and k 
lines, along with the Ca II H and K lines, leads 
to a temperature profile of the upper photo- 
sphere and lower chromosphere (Basri et al., 
1981). Chromospheric densities can be deduced 
from the relative strengths of the lines in the 
C II (UV 0.01) multiplet (Stencel et al., 1981; 
Carpenter, 1984). The extent of the chromo- 
sphere can then also be deduced from the ab- 
solute intensities of the C II lines. For Betel- 
geuse, the chromosphere has a temperature of 
6000 to 8000 K and a thickness on the order of 
the radius of the star. The variation of relative 
velocity within the chromosphere can be ob- 
tained from a study of the line shifts of the 
numerous Fe II lines (Carpenter, 1984). 

Ultraviolet observations of other M, S, and 
C stars we summarize briefly. Barium stars 
show strong emission lines of Mg II similar to 
K giants of the same temperature. Early R stars 
show no ultraviolet chromospheric indicators, 
while later R stars (R8) have weak emission 
lines of Mg II and Fe II (Eaton et al., 1985). 

S stars apparently show a large range in 
strength of chromospheric emission, but very 
little has been published. One special use of 
IUE is to search for hot companions to M, S, 
and C stars. Several types of these peculiar red 
giants have been observed in the short-wave- 


length region for possible white-dwarf compan- 
ions, and a few have been discovered (Dominy 
and Lambert, 1983; Bohm-Vitense and 
Johnson, 1985). 

The first attempt to observe N-type carbon 
stars was negative (Querci et al., 1982): neither 
emission lines nor continuous flux were re- 
corded; the observed flux was lower than 
* 10~ 16 erg/cm 2 A (3a flux limit). Later obser- 
vations with IUE detected weak UV spectra, 
with some continuum down to at least 2850 A, 
for seven bright N-type carbon stars (Johnson 
and O’Brien, 1983; Querci and Querci, 1985; 
Johnson and Luttermoser, 1986). Longward of 
2850 A, the spectrum is basically an absorption 
spectrum — apparently the ultravioletward ex- 
tension of the photospheric spectrum; short- 
ward of 2850 A, only emission lines (C II, Mg 
II, Al II, and Fe II) are detected. The contin- 
uous flux is, however, much weaker than in M 
giants, which implies either that the cool C stars 
have weak chromospheres (low temperatures or 
densities) or that an important opacity source 
is located above the emission-line-forming 
region (see the conclusions obtained in this 
sense in Chapter 2). Yet the strength of the 
chromosphere (e.g., L^/L^) decreases rap- 
idly with spectral type among the M giants. Will 
M7 or M8 giants, which have effective temper- 
atures as low as those of N-type carbon stars, 
have chromospheres as weak? Why are chro- 
mospheres of carbon stars so weak? Is it due 
to lower temperatures or different chemical 
composition or more overlying absorbing ma- 
terial? Observation time on the Space Telescope 
is necessary to observe these stars with enough 
light signal and at various phases to make 
decisive progress on this problem. 

A violet-ultraviolet color-color diagram (Mg 
II-V vs. IUE-V) has been developed and used 
to compare spectrophotometric aspects of IUE- 
LWR observations for K, M, S, and C stars 
(Eaton et al., 1985), as is shown in Figure 1-2 
(see also the section on Intrinsic Properties). As 
anticipated, red giants and supergiants have low 
UV continuous fluxes and display lines from 
only such low ionization stages as Mg II and 
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Figure 1-2. A two-color diagram for stars showing an Mg II index (measuring ratio of the flux in 
a 30 A band at 2800 A to the flux at V) against an IUE color (the ratio of flux between 2585 and 
3200 A (Mg II omitted) to the flux at V). G and K dwarf stars are shown as circled dots, G, K, and 
M giants and supergiants as letters, R stars as asterisks, and N stars as open circles (from Eaton et 
al., 1985). 


Fe II (see Chapter 2), indicating the lack of a 
transition region. 

On the low-energy side, (i.e., th efar-IR), 
the opaqueness of the terrestrial atmosphere 
necessitates the use of satellites. The outstand- 
ing of these— the Infrared Astronomical Sat- 
ellite (IRAS), a joint project of the United 
States, the Netherlands, and the United King- 
dom-produced a wealth of valuable data dur- 
ing 1 year of use (Neugebauer et al., 1984). 
Another one, the Infrared Space Observatory 
(ISO), is planned to be launched by the Euro- 
pean Space Agency (ESA) in 1992 to perform 


imaging, spectroscopic, photometric, and po- 
larimetric investigations in the IR spectral range 
from 3 to 200 /im during an operational lifetime 
of about 18 months. 

The data of IRAS are presently under reduc- 
tion; the first results on cool giants were 
presented at the University of California at Los 
Angeles, meeting, “Mass Loss from Red 
Giants,” held in June 1984, and at the Calgary 
Workshop in 1986. Physical parameters of the 
stellar outer layers are deduced. The low- 
resolution spectrometer (LRS) made spectra 
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between 8 and 22.5 /an with a spectral resolu- 
tion varying between 15 and 30. A continuous 
sequence is found by OInon et al. (1984) from 
classical Miras to the very red Miras like OH 
17.7 - 2.0 or OH 26.5 + 0.6, which are very 
cool objects (Werner et al., 1980); only the 
outermost dust shell is visible at all LRS 
wavelengths, and no absorption features are 
observed. As it is already known from the radio 
observations (Smolinski et al., 1977), Stickland 
(1985) found no contribution to the IR excess 
from free-free radiations on hypergiants. 
Analyzing the IRAS data of many RCB stars, 
Schaefer (1985) finds a bimodal distribution 
among them (see the section Eruptive 
Variables: The RCB Stars). Maps at 60 and 100 
fim around M supergiants are also obtained 
with IRAS. For a Ori, ‘‘shells” appear that are 
not centered on the star. Are these shells ejected 
remnants or parent cloud to a Ori (Wesselius, 
1984)? 

Long-term monitoring of a few oxygen-rich 
Mira variables has already begun to provide in- 
formation about the temporal variation of the 
peak intensity and the shape of some radio 
maser lines (Barcia et al., 1985); their profile 
and their three-peak structure are in agreement 
with theoretical calculations (Zhou and Kaifu, 
1984). Moreover, the mean intensity ratio be- 
tween the v = 1 and v = 2 SiO lines seems to 
agree qualitatively with the theoretical predic- 
tions of Bujarrabal and Nguyen-Q-Rieu (1981). 
Let us also point out some scarce radio flares 
detected on o Cet (Boice et al., 1981), on R Aql 
(Woodsworth and Hughes, 1973; Bowers and 
Kundu, 1981), on V Cyg (Querci et al., 1979), 
and on a Ori (Kellermann and Pauliny-Toth, 
1966; Seaquist, 1967). We stress the importance 
of monitoring radio, visible, and far-UV con- 
tinuum for the detection of such nonthermal 
flares and the analysis of their energy distribu- 
tion. The radio continuum spectrum of a Ori 
obtained with the very large array (VLA) was 
fitted by a power law S v ~ v a with a = 1.32 
(Newell and Hjellming, 1982). A model of 
chromosphere has been deduced with an elec- 
tron temperature of 10000 K and an extension 
of 4 R*. 


In this brief review, we have summarized the 
first photometric attempts made in the far-UV 
and X rays and in the far-IR and radio frequen- 
cies to show the relative lack of monitoring in 
these spectral ranges and the urgency in foster- 
ing such observations. The reasons for the 
paucity of observations are mainly that: (a) the 
stellar flux is very low, and it can be registered 
with current instruments for only a few bright 
stars; and (b) the flux may be strongly variable 
(why not?), and the phase of the star’s max- 
imum has not always been considered. 

Space Motion and Space Distribution 

Although detection of M, S, and C stars has 
mainly been due to objective-prism surveys, 
new techniques involving photoelectric spectral 
scanners or photography of star fields through 
narrowband filters are adopted to identify cool 
stars in crowded fields such as the fields of rich 
star clusters (Palmer and Wing, 1982). These 
filters are selected from Wing’s (1971) eight- 
color photometry previously used to identify 
carbon stars (e.g., Wing and Stock, 1973) or 
to classify M supergiants (White and Wing, 
1978). 

In the solar neighborhood, only 1 percent 
of all the stars are giants, and the G, K, M, and 
correlated types of giants make up at most 5 
percent of this 1 percent (Jaschek, 1985). The 
apparent and space distributions in the solar 
neighborhood are well known through early re- 
views on the question. M stars are found in 
large numbers in the galactic center direction 
(2000 per square degree according to Blanco, 
1965), while the C stars avoid this direction. On 
the other hand, the C stars are strongly concen- 
trated toward the galactic equator, compared 
to the late M stars. The majority of carbon stars 
appear to be spiral-arm objects. Among the car- 
bon stars, the R stars have smaller concentra- 
tion toward the galactic plane than the N stars; 
moreover, the R stars are found in great num- 
bers where no N stars are present. The disper- 
sion of the early R stars above the galactic plane 
indicates that they are representative of a pop- 
ulation that is older and less massive than the 
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N stars (Stephenson, 1973). Alksne and Ikau- 
nieks (1971) note an apparent concentration of 
later (R5-R9) stars toward the galactic equator 
similar to that for Mira and semiregular C stars, 
while type Lb (irregular variations of light) 
variables are concentrated more to the equator. 
The radial distribution and hydrogen deficien- 
cy appear to be correlated: R stars such as HD 
100764 have low radial velocity (5 km/s) and 
do not suffer the severe hydrogen deficiency 
seen in supergiant R stars such as the high- 
velocity star, HD 182040, in which the Balmer 
lines are absent (Sanford, 1944). 

The nonvariable HdC stars have strong 
similarities with the RCB (variable HdC stars): 
chemical composition, space distribution (|/ij 
from 0.008 to 0.024 and |^| from 0.010 to 
0.035). These values of radial velocities and 
proper motions are the characteristics of high- 
velocity stars. With space velocities in excess of 
400 km/s, these stars are members of the halo 
population. The He stars (hotter HdC stars) 
seem to be spread out farther from the galac- 
tic plane than the cooler ones, as is also in- 
dicated by their somewhat larger radial veloci- 
ties (Warner, 1967). In summary, all the cool 
nonvariable HdC stars appear to be closely con- 
centrated to the galactic center, whereas the 
RCB stars are not noticeably concentrated to- 
ward the galactic plane and galactic center. The 
He stars are widely distributed. 

As for the M stars, if the late types (M7- 
M 10) are concentrated in directions close to the 
galactic center, as quoted above, M5-M6 types 
are not as concentrated, and M2-M4 types are 
nearly equal in number in the galactic center 
and the anticenter directions. Moreover, the 
early M types present local groupings along the 
galactic equator. It is concluded by several 
authors (quoted by Alksne and Ikaunieks, 
1971) that half or more of early M-type stars 
are members of the galactic spiral arms, while 
late M types belong to the disk population. An 
illustration is provided by the careful study of 
722 M giant stars (with a few S and C stars) 
in the Perseus spiral arm (Ichikawa, 1981). As 
pointed out by Gehrz et al. (1981), OH/IR stars 
are strongly confined to the galactic plane and 


are usually found at such large distances (6 to 8 
kpc) (Baud et al., 1981; de Jong, 1983) that in- 
terstellar extinction at infrared wavelengths 
cannot be neglected. Space densities ( n ^ and 
effective heights of the space distribution ( z e ) 
from the galactic plane of late-type M giants 
in the solar neighborhood observed by various 
authors are listed by Mikami and Ishida (1981; 
see their Table 3). These authors compare them 
with their statistical estimates based on the 
Two-Micron Sky Survey data (Neugebauer and 
Leighton, 1969) and find them to be in good 
agreement; moreover, they determine the space 
density and distribution of M supergiants 
(Table 1-2). 

S stars are as numerous in the direction of 
the galactic center (like the late M stars) as they 
are concentrated in the galactic spiral arms. 

Discussing the kinematics and stellar popu- 
lation of the Ba II stars, Eggen (1972a) con- 
cludes that their space motions resemble those 
of an old disk population of 1.0 to 1.5 M© 
stars. In their sample of Ba II stars, Catchpole 
et al. (1977) found a difference in velocity 
dispersion from Ba-strong versus Ba-weak 
stars, which is explained by Rocca-Volmerange 
and Audouze (1979) as a reflection of an age- 
versus-abundance correlation. 

Other investigations studied the distribution 
of late-type red giants in the direction of the 
central regions of our Galaxy and of other near- 
by galaxies, the Large and Small Magellanic 
Clouds, and other galaxies such as M31. Blan- 
co et al. (1978, 1984) survey Baade’s Window 
(that surrounds the globular cluster, NGC 
6522), which constitutes a uniquely complete 
sample of the red-giant population of the galac- 
tic nuclear bulge. The striking feature is a very 
large number of M giants, but no C or S stars 
were found, even though they should have been 
detected had they existed. Such a result con- 
firms a previous survey in the Sagittarius I 
region (Baade, 1951), in a direction closer to 
the galactic center than Baade’s Window. 
Therefore, the ratio of the number density of 
carbon stars to the number density of M 
giants— the C/M ratio— is negligibly small in 
the galactic bulge. 
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Table 1-2 

Space Densities and Distributions of Nearby M Giants and Supergiants* 


Spectral 

Type 

Member Density 
( 1 0” 6 stars/pc 3 ) 

Effective Height 
(pc) 

K0-K4 

40.0 

170 

K5-M1 

3.5 

200 

M2-M4 

3.1 

310 

M5-M6 

0.9 

390 

>M7 

0.6 

380 

M supergiants 

0.054 

210 


•Adapted from Mikami and Ishida (1981). 


The distribution of carbon and M-type 
giants in the Magellanic Clouds (Aaronson, 
1984) and nearby galaxies is important for in- 
terpreting AGB evolution (e.g., Iben, 1984; 
Wood, 1985) and for understanding fundamen- 
tal galactic properties, such as the metal abun- 
dance (C/M is strongly correlated with this 
quantity (e.g., Scalo, 1981)), the distance to the 
galaxy, and the full AGB luminosity function 
for the galactic field studied. Outstanding is the 
lack of late-type M stars in the central regions 
of the Small Magellanic Cloud (SMC), remind- 
ing us of the almost complete lack of C stars 
in the bulge of our Galaxy (Blanco et al., 1978). 
Also, differences in the C/M ratio between the 
SMC and the Large Magellanic Cloud (LMC) 
are noticeable. A survey of high degree of com- 
pleteness for carbon stars and M-type giants in 
the Magellanic Clouds (Blanco and McCarthy, 
1983) indicates that the ratio of the surface den- 
sity of carbon stars and that of giants of type 
M6 or later in the SMC varies from 19.2 ± 0.8 
at the center to 4.7 ± 0.4 in the periphery, while 
the C/M ratio is 2.2 ± 0.1 throughout the 
LMC. It is suggested that the mixture of stellar 
populations is roughly uniform across the 
LMC, but not in the SMC. Recalling the neg- 
ligible value of the C/M ratio in the galactic 


nuclear regions, Blanco and McCarthy (1983) 
note that their population includes older and 
less massive red giants than those observed in 
the Magellanic Clouds. A recent review of the 
distribution and motions of red giants in the 
Magellanic Clouds is by Catchpole and Feast 
(1985). 

A field in M31 has been investigated by 
Richer, Crabtree, and Pritchet as quoted by 
Richer (1985) in which C/M is about 0.12, im- 
plying a metal abundance higher than that in 
either the LMC or NGC 300, but lower than 
that in the galactic center. In addition, there ap- 
pears to be a deficiency of luminous AGB stars 
as in the LMC. 

To conclude on the spatial distribution, we 
note that carbon stars have been discovered in 
all dwarf spheroidal galaxies (Aaronson et al., 
1983). 

Intrinsic Properties 

Red giants and supergiants (M, S, C, Ba II, 
MS, CS, CH, and CN stars, etc.) occur in a 
large range of stellar masses and ages because 
they are found in halo, old disk, and young disk 
populations. The observational data indicate 
that they have mixed nuclearly processed ma- 
terial from the core to the atmosphere for 
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luminosities such as -10 = M bol = + 3, ages 
from - 10 7 to 10 10 yr, and a broad range of 
masses from ~1 M© to = 20 Mq, with a 
mean mass of 1.2 Mq (Scalo, 1976, 1984). A 
summary of the physical properties (luminosi- 
ties, temperatures, masses, abundance of ele- 
ments, population types, etc.) of red giants of 
the disk population, as well as their character- 
istic spectral features, is given in Table 1-3 
(from Scalo, 1981). 

Absolute Magnitudes and Colors. The N-type 
stars and the various S-, MS-, and CS-type stars 
belong to luminosities -4 > M bol > -6 (Richer, 
1981; Bessel et al., 1983; Lloyd-Evans, 1983). 
They are asymptotic giant branch (AGB) stars 
where the helium shell flashing occurs in the in- 
termediate mass stars. Ishida (1960) gives an 
M of -1.5 to -2.0 for the N-type stars. The 
R5-R8 stars are also very bright but probably 
not well classified (Eggen, 1972a). The hotter 
R stars (R0 to R4 types), the Ba II and CH 
stars, fall in a fainter range of luminosity ( + 3 
£ M bo/ S -3); they are not so high on the 
AGB as the N-, S-, MS-, and CS-type stars 
quoted above. They are probably older, lower 
mass stars. With the well-known relation 
among the emission core of Ca II H and K lines 
and the luminosity of late-type stars (Stratton, 
1925; Wilson, 1959; Wilson and Bappu, 1957; 
Warner, 1969), Richer (1975) obtains the range 
in M between -0. 1 and +1.1 for the hot R 
stars. Bright members of this class are HD 
182040 with V = 6.99, B-V = 1.07, U-B = 
0.65, R-I = 0.55, and I-L = 0.44 and HD 
156074 with F= 7.61, B- V= 1.1, R-I = 0.3, 
and I-L = 2.2. 

Although the great majority of these stars 
are variable, we must keep in mind that some 
are considered nonvariable, such as the early 
R stars, which are predominantly nonvariable 
giants (see the section Nonvariable Stars). 
Distinctions in absolute magnitudes or in col- 
ors are noted between variables and non- 
variables of the same type or same class. The 
variable R stars are redder than the latest R 
nonvariables (Vandervort, 1958). The non- 
variable R and N stars both have a continuum 


1104 magnitude M(104) 1 -1.2 (Baumert, 
1974); they are approximately 1.5 to 5 mag 
fainter than the variables, while the early non- 
variable R stars are 2 mag fainter than the late 
R nonvariable stars. Baumert (1974) notes that 
the brightness increases in the sequence non- 
variables, irregular Lb giants, semiregular SRb 
giants, Miras, and semiregular SRa giants. (See 
the section Variability Types of Giants and 
Supergiants for the definition of these various 
variables.) No clear correlations were observed 
between absolute magnitude and CN strength 
or temperature. 

To conclude on the “classical” carbon stars, 
we stress that the early R stars are physically 
distinct from the later carbon stars. Such a 
dichotomy is not only striking from visible in- 
dices, as mentioned above, but also from the 
(R-I) index (Eggen, 1972b) and from a far-UV 
index (Eaton et al., 1985). From the latter, 
the carbon stars seem to be separated into 
three groups in the color diagram (Mg 
II-V)/(IUE-V) (Figure 1-2). Early R stars fall 
into two groups among the late G and early K 
giants; late R stars (R5 to R8) distinctly overlap 
with the reddest K and M giants, joining 
smoothly with the N stars at their bluer end. 
The N sequence extends 2.5 mag farther to the 
red. On the other side, Figure 1-2 illustrates the 
relation between color and chromospheric emis- 
sion for C stars and other giants: the fraction 
of stellar luminosity radiated in the 
chromospheric Mg II lines decreases with ef- 
fective temperature. Low chromospheric emis- 
sion characterizes the N and R stars (also the 
G dwarfs) and a few K stars, whereas most of 
the K and M giants are located above the car- 
bon sequence by up to 2 mag. 

All the properties of the HdC stars are well 
summarized by Warner (1967), Wallerstein 
(1973), and Hunger (1975). 

The atmospheric parameters of the RCB 
stars indicate high luminosity; they are as lu- 
minous as other HdC stars (Warner, 1967). 
These stars are supergiants (Orlov and 
Rodriguez, 1977). Herbig (1958) gives M bgl ~ 
-4.0 to -6.0 for the prototype star, R CrB. In 
his Table 20, Glasby (1968) gives 25 RCB stars 
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Table 1-3 

Physical Properties and Characteristic Spectral Features of Disk Population Red Giants*, 
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Symbols: N = K giant normal abundance, A = Absent, + or - = somewhat enhanced or deficient with respect to normal K giants, + + or — = strongly 
enhanced or deficient, ? = unknown or uncertain, < ) = only a small fraction of the spectral type possess the anomaly. 

The straight line across the table delineates the "hot" stars, early R and Ba II stars (above the line) for which -4000 < T co)or <5500 K, and the cooler 
ones (below the line) for which - 2000 < T co|or < 3800 K. ° 


with maximum and minimum magnitudes and 
spectral types. Alksne and Ikaunieks (1971) 
show the spectral class, Howarth (1976) shows 
the maximum magnitude and spectral class, and 
Sherwood (1976) adds some stars to Glasby’s 
list. Among them, four are observed in the 
Large Magellanic Cloud: W Men, SY Hyi, HV 
5637, and HV 12842, with the absolute mag- 
nitude at maximum being -4.8, -6.6, -3.2, and 
-4.9, respectively. Duerbeck and Seitter (1982) 
give a list of RCB stars with their spectral type, 
their V, B- V, and U-B values, and their pulsa- 
tional period. 

Various estimates have been made concern- 
ing the helium stars; they are based on the 
strength of interstellar lines and proper mo- 
tions. For MV Sgr, Klemola (1961) gives M y 
* -2.0 to -4.0, and consequently, M bo( ^ 
-4.0 to -6.0; Herbig (1964) gives a value of -5.0 
for M bgr 

The analysis of the interstellar D lines of the 
typically hydrogen deficient C (HdC) non- 
variable star, HD 182040, shows that this R star 
is highly luminous and confirms Warner’s 
(1967) conclusion that the HdC stars are super- 
giants (Utsumi and Yamashita, 1971). From 
observations of the Ca II K emission cores and 
the Wilson-Bappu effect, Richer (1975) gives 
M » -4.1 for the three nonvariable HdC 
stars: HD 137613, HD 173409, and HD 182040. 
Because their atmospheric parameters indicate 
high luminosity, they are as luminous as their 
variable analogs (Warner, 1967). The hydro- 
gen-deficient (supergiant) R stars are physical- 
ly distinct from the N and “classical” giant R 
stars. The He, RCB, and nonvariable HdC 
stars have a common range in values of mean 
M value -2.0 to -5.0. 

V 

MacConnell et al. (1972) estimate that 1 per- 
cent of the G and K giants are Ba II stars; 
therefore, the Ba II are Population I giants 
(classical and very weak line objects), and their 
absolute magnitudes are similar to those of G 
and K giants (Baumert, 1974). The classical Ba 
II stars are all binaries, but this may not be true 
for the mild Ba stars. Bond and Neff (1969) and 
Gow (1976) showed the excess absorption in the 
violet region of the Ba II stars: a controversy 


developed about the reality and the identity of 
this feature (McClure, 1984, and references 
therein). The difficulty in placing the Ba II stars 
in the HR diagram is mainly one of determin- 
ing their luminosities and Morgan-Keenan spec- 
tral classification. From his analysis of many 
studies, McClure (1984) concludes that most Ba 
II stars have “normal or slightly above normal 
(<Af> - 0.0 mag) G and K giant luminosi- 
ties.” However, the classical bright Ba II star, 
£ Cap, has a very high luminosity: A/ y « -3. 
This implies that the star has a large mass or 
that it could be in a helium shell flashing stage 
of evolution on the second AGB (McClure, 
1984). 

For the CH stars, we have M(104) = -3.9 
as derived by Baumert (1974), who believes that 
the CH stars measured by Wing and Stock 
(1973) are at the extreme tip of the red-giant 
branch. 

Among the brightest stars are early M-type 
supergiants (M0 to M4 types) with -8 = M y § 
-5. St others and Leung (1971), Humphreys 
(1978, 1979), and Cowley and Hutchings (1978) 
collected such bright objects in clusters, nearby 
galaxies, and the Large Magellanic Cloud. 
Their M m would be -7 to -10 (Glass, 1979), 
like the M supergiant near 30 Doradus, the car- 
bon stars in NGC 1783 or NGC 2477, or q Cas 
and VX Sgr (Sargent, 1961). The most lumi- 
nous of the nearby supergiants is \i Cep. With 
its M - -8 and its M, « -10, it is about one 
million times brighter than the Sun. All these 
stars are also considered to be super- 
supergiants, a term introduced by Feast and 
Thackaray (1956), otherwise called hypergiants 
(Van Genderen, 1979). 

For the M Miras, Clayton and Feast (1969) 
conclude that the absolute magnitudes vary 
smoothly within the period: M m = -3.0 to 
-1.0 at maximum light and = -1.5 to 0.1 
at mean light intensity, with periods varying 
from - 180 to -500 days. The shorter period 
Miras (-130 days) deviate from this period/ 
luminosity relation. They have M m - -1.6 
and Mj = 0.1. Robertson and Feast (1981) 
show that the stellar luminosity decreases with 
period for the galactic Miras, and Glass and 
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Lloyd-Evans (1981) reach the same conclusion 
for the Miras in the Large Magellanic Cloud. 

Information on colors and luminosities are 
best displayed as an observational HR diagram, 
and several attempts to construct such diagrams 
for galactic red giants have been made (e.g., 
Scalo, 1976; Tsuji, 1981c). A better result 
comes from the Magellanic Clouds, where the 
known distance permits absolute luminosities 
to be determined (Wood et al., 1983). (See the 
section Aspects of Evolution of Long-Period 
Variables (LPV’s) in the Magellanic Clouds.) 

Effective Temperatures. Effective temperatures 
can be measured or estimated for at least a few 
stars from the relation F = , where F is 

the integrated energy flux at the stellar surface. 
This flux can in turn be derived by any of 
several means, all depending on some variation 
of the relation between observed energy flux 
(/), emitted flux (F), and angular diameter (0): 
f = F(0 2 / 4), where the flux is either mono- 
chromatic or integrated. That is, the effective 
temperature can be inferred for any star for 
which the angular diameter and total integrated 
(bolometric) flux (outside the Earth’s atmo- 
sphere) can be measured. This latter quantity 
is often difficult to obtain for late-type stars 
since the bolometric correction is so large. In 
fact, the desired information — complete spec- 
trophotometry — is seldom available for the M, 
S, and C stars (Strecker et al., 1979, for M 
stars; Goebel et al., 1980, for a carbon star; and 
Augason et al., 1986, for S stars), and fluxes 
obtained from broadband colors are used in- 
stead (e.g., Ridgway et al., 1980a). 

Angular diameters are measured or inferred 
by several means: 

1 . They can be directly measured by lunar 
occultation for stars near the ecliptic, and 
this is the principal source of our current 
data (e.g., Ridgway et al., 1980; Beavers 
et al., 1982). 

2. For a few stars, values of angular 
diameters are available from Michelson 
interferometry or from speckle in- 
terferometry (e.g., Balega et al., 1982; 
Bonneau et al., 1982; see also Chapter 2). 


3. Angular diameters can be indirectly in- 
ferred from the relation between visual 
surface brightness and unreddened ( V-R ) 
color (Barnes and Evans, 1976; Barnes 
et al., 1978; Eaton and Poe, 1984) al- 
though the validity of this relation for the 
coolest giant stars has not been estab- 
lished beyond doubt. 

4. Substitution of the monochromatic flux 
from an appropriate model for the flux 
at the stellar surface allows one to infer 
the angular diameter of any selected star 
for which the monochromatic flux can 
be observed (see above formula). 

If the integrated flux can be observed, the 
effective temperature can also be calculated. 
Since the infrared region is often least influ- 
enced by molecular bands, it offers the best 
choice for application, and the method is gen- 
erally called the method of infrared photometry 
(Blackwell and Shallis, 1977; Blackwell et al., 
1980). For the coolest red giants generally, this 
method provides the best present information 
and has been used to define a temperature scale 
for M giants (Tsuji, 1981a), for S stars (Auga- 
son et al., 1986), and for N-type carbon stars 
(Tsuji, 1981b). 

A final method of determination of effec- 
tive temperature, that of model fitting, does not 
depend on the angular diameter. Instead, one 
directly compares computed and observed flux 
curves and chooses the effective temperature of 
the best-fitting model. Widely used in hot and 
intermediate stars, the method has not been 
much used in cool stars (e.g., Querci et al., 
1974; Bouchet et al., 1983; Steiman-Cameron 
and Johnson, 1986) because of the lack of com- 
plete spectrophotometry and the deficiencies of 
the models, although it has been applied to 
calibrate a scale of effective temperatures 
for M giants (Tsuji, 1978). As these deficien- 
cies are corrected, this method will assume 
greater importance. In this scope, the com- 
parison between observed and synthetic spec- 
tra should give an effective temperature (e.g., 
Querci and Querci, 1976). However, it remains 
the problem of uniqueness of the solution. 
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Table 1-4 

Effective Temperatures [T J for M Giant Stars* 

•ff 


Spectral 

Type 

T 

c 

(K) 

V-K 

Color 

T « (K)1 

T « < K > 2 

r, f (K) 3 

MO 

3750 

3.78 

3895 

3900 

3800 

Ml 

3640 

4.02 

3810 

3800 

3865 

M2 

3530 

4.30 

3730 

3700 

3800 

M3 

3400 

4.64 

3640 

3600 

3640 

M4 

3250 

5.10 

3560 

3500 

3460 

M5 

3000 

5.96 

3420 

3300 

3310 

M6 

2600 

6.84 

3250 

3200 

3280 

M8 

— 

— 

— 

2300 

— 


1. Effective temperature from broadband colors and lunar occultation angular diameters (Ridgway et al., 
1 980a, b). 

2. Effective temperature from model fitting (Tsuji, 1978). 

3. Effective temperature by method of infrared photometry (Mira variables are not considered) (Tsuji, 
1981a). 

*Color temperature is on the Wing system (Ridgway et al., 1980). 


In the past, a rough method was to use a 
spectroscopic temperature or excitation temper- 
ature for the effective temperature. 

Results for M giant stars are collected in 
Table 1-4. Although these are the best current 
values, there is still considerable uncertainty, 
and much more work is needed. A graphic 
comparison of the results obtained for M stars 
by lunar occultation (Ridgway et al., 1980a, b), 
by infrared photometry (Tsuji, 1981a), and by 
the method of model fitting (Steiman-Cameron 
and Johnson, 1986) is shown in Chapter 7 
(Figure 7-14). As yet, there is no reliable tem- 
perature scale for M giants or supergiants later 
than M6, where most stars are Mira variables, 
the very red ones being OH-1R objects. 

Several obstacles hinder the decisive deter- 
mination of effective temperatures for the R 


stars. No angular diameters are available, mol- 
ecules distort all broadband colors, the method 
of infrared photometry fails because of an ex- 
cess of flux in the L band compared to G-K 
giants, blackbody colors are inconsistent, and 
insufficient spectrophotometry is available for 
model fitting. These and other problems have 
been exhaustively discussed by Dominy (1982, 
1984), who finds that effective temperatures for 
three R0 stars lie in the range 4500 5 T rf g 
4850 K. However, some R0 stars may be even 
hotter. Dominy finds that R5 stars are nearly 
as hot as R0 stars. The temperature scale for 
the later R stars is presently quite uncertain, and 
these stars may not be closely related to the ear- 
ly R stars. 

Effective temperatures for N-type carbon 
stars have been notoriously difficult to obtain, 
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and even now they are uncertain except for a 
few stars. Angular diameters have been mea- 
sured for six non-Mira N-type carbon stars, of 
which five (TX Psc, X Cnc, AQ Sgr, SZ Sgr, 
and TW Oph) are irregular variables and one 
(Y Tau) is an SRa variable (Ridgway et al., 
1977; Walker et al., 1979; Ridgway et al., 
1980b; Ridgway et al., 1982). Although one star 
(TW Oph) is so heavily reddened that the ef- 
fective temperature is unreliable, values for the 
other stars form the cornerstone of any further 
effort. 

Tsuji (1981b) has obtained values of effec- 
tive temperatures of 31 N-type carbon stars 
(SRb and Lb variables) by the method of in- 
frared photometry (Blackwell et al., 1980). 
Angular diameters were deduced from 
previously published models (Querci et al., 
1974; Querci and Querci, 1975) and L band 
photometry (Noguchi et al., 1977), corrected 
for the effects of molecular bands (however, see 
discussion by Bouchet, 1984a). Combined with 
complete broadband photometry, these yield 
effective temperatures. For the stars in com- 
mon, the results of Tsuji agree with those ob- 
tained from measured angular diameters, and 
we believe that these are the best values current- 
ly available. The effective temperatures are all 
confined within the range 2400 = T ff = 3200 
K, and most fall between 2600 and 3100 K. It 
is well known that the “temperature classes” 
of the C star classification are not well cor- 
related with effective temperatures (e.g., Tsu- 
ji, 1981b); Tsuji (1985a) presents arguments to 
explain this situation. 

A scale of effective temperature for non- 
Mira S stars has been deduced (Augason et al., 
1986) both by model fitting and by the method 
of infrared photometry, and the complicating 
effect of the unknown C/O ratio is revealed. 
There is also no temperature scale for S or C 
Mira variables. 

Once fundamental effective temperatures 
have been determined for a sufficient number 
of calibrating stars, it will be useful to calibrate 
a color temperature derived from a carefully 
chosen set of filters (Wing, 1985). When this 
is possible for all types of peculiar red giants, 


values of effective temperature can be routine- 
ly obtained. Many effective temperature in- 
dicators are presently taken into account among 
the cool stars. Since the bolometric flux is very 
sensitive to temperature variations and the in- 
frared monochromatic flux is not, the color in- 
dex, m bol ~L , can be used as a temperature in- 
dicator (proportional to 7^ ff ) (Tsuji, 1981a, 
1981b, 1985a). Another is the I(104)-L(400) in- 
dex by Wing and Rinsland (1981) because the 
two spectral regions, 104 (1.04 pan) and 400 
(4.00 jrni), are free from molecular absorptions 
in the cooler stars. 

For cool carbon stars, Tsuji (1981c) shows 
that all broadband photometric colors correlate 
fairly well with his effective temperature scale. 
Interestingly, although the intrinsic ( R-I ) col- 
or is nearly the same for the N-type irregular 
stars, these have a range of temperature, and 
(R-I) is therefore an indicator, not of 
temperature, but of interstellar reddening. The 
(I-L) color is suggested by Tsuji as the best 
temperature discriminant for both carbon stars 
and M giant stars. Wing (1985) also suggests 
the use of the I-K index. 

Although the intrinsic infrared colors, (J-K) 
and ( H-K ), for cool carbon stars correlate with 
effective temperature and have been used to in- 
fer effective temperature, the range in both col- 
ors is relatively small (-0.6), and the scatter 
in (J-K) is large. These colors do not appear 
to be reliable temperature indicators. One 
reason for the scatter might be the imprecision 
of the effective temperatures. On the other 
hand, the claim has been made, on the basis 
of comparative studies of carbon stars in the 
Magellanic Clouds and the Galaxy, that the in- 
frared colors of carbon stars are primarily 
determined by strengths of molecular bands 
and only secondarily by effective temperatures 
(Cohen et al., 1981). Perhaps the answer lies 
somewhere between; infrared colors may well 
be influenced both by effective temperature and 
by chemical composition. 

The temperature range covered by the RCB 
stars as a group is very large. For example, RS 
Tel has a spectral class as late as R8 with a sur- 
face temperature of 2500 K. Two very hot 
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stars — the B star, MV Sgr, and the 09 star, 
V348 Sgr — are also considered as RCB stars be- 
cause of their variability behavior. 

Radii. The giants, and chiefly the supergiants, 
are enormously distended. Most of them have 
sizes approaching that of our solar system. 
Some examples follow. If we adopt an effec- 
tive temperature of 3000 K for fi Cep, one of 
the most luminous supergiant M stars, we find 
a radius of some 1590 times the solar value; this 
means that the photospheric layers of the star 
fill in our planetary system out to the orbit of 
Saturn! The extent of the Cep outer envelope 
layers are not yet measured nowadays. The 
photosphere of the supergiant Betelgeuse (a 
Ori: M2 lab) is roughly the size of Jupiter’s or- 
bit ( R ~ 650 R 0 ). Polarimetric observations 
give a radius of the a Ori outer envelope of 
roughly 3700 R Q . (Our planetary system would 
be filled out to the Uranus orbit!) Observed 
from the Earth, the outer envelope of a Ori has 
an apparent diameter equal to 1/10 of that of 
the Sun. When we consider the HR diagram, 
we have to bear in mind that two supergiants 
of the same luminosity but with different 
temperatures do not have the same diameter. 
For example, the M, S, and C stars have radii 
100 times that of a main-sequence B star. The 
cool hypergiants are still more extended stars. 
Finally, the S Mira, x Cyg, has a photospheric 
radius R ~ 240 R 0 (Hinkle et al., 1982). 

The direct measurement of angular diameter 
of single stars uses the lunar occultation tech- 
nique (e.g., Ridgway et al., 1980a, 1980b; Bea- 
vers et al., 1982) and numerous interferometric 
techniques such as speckle interferometry (e.g., 
see the review of Dainty, 1981, and references 
therein), as noted in the previous section. A 
photometric approach may also furnish the stel- 
lar radius; if one knows the absolute bolometric 
magnitude, M bo[ (or the luminosity), and the 
effective temperature, 7 ff , the radius of the 

star comes from: L = 4 n R 2 a T 

eff 

Description of the geometrical shell exten- 
sion is given by M. Querci (this volume). As 
an example, let us add that, using the very long 


baseline interferometer (VLBI), Bowers et al. 
(1980) found some masing regions in OH/IR 
stars from > 0.5 " to > 0.04 ". The extent of the 
1612-MHz region can be quite large, from 1000 
to 10000 AU. The apparent diameter is consis- 
tent with the general predictions of the maser 
models (Goldreich and Scoville, 1976). 

From the knowledge of the stellar radius and 
mass, the acceleration of gravity in a stellar at- 
mosphere may be determined. In these stars, 
the g eff values are rather badly known since the 
masses are linked to an uncertain evolution 
theory (Wood, 1985); they are considered to be 
0.1 < i* < 10. Determination of the surface 

°eff 

gravity for these peculiar red-giant stars re- 
mains a pressing unsolved problem. Much in- 
sight can be gained by examining the exhaus- 
tively discussed parallel problem for the well- 
observed K2 giant, Arcturus (Trimble and Bell, 
1981). (Examples of indirect methods of esti- 
mating the surface gravity of late-type giants 
are given in Chapter 7.) 

Aspects of Evolution of Long-Period 
Variables (LPV’s) in the Magellanic Clouds 

Studying long-period variables (i.e., the 
most extreme — in luminosity— red giants and 
supergiants) in the Magellanic Clouds releases 
us from the problem of distance determinations 
(hence, absolute magnitude determinations) in 
comparison to LPV’s in the galactic disc for 
which existing distance determinations are gen- 
erally based on statistical parallaxes (since they 
are not known to be related to open clusters). 
Also, the LPV’s in the Magellanic Clouds cover 
a wide range in mass and luminosity and are 
generally free from significant interstellar red- 
dening (Wood, 1982). 

Wood et al. (1983) obtain IR photometry 
(JHK) and low-dispersion red spectra of 90 
LPV’s in the Small Magellanic Cloud (SMC) 
and the Large Magellanic Cloud (LMC), com- 
pleting the IR photometry data of Glass and 
Feast (1982 and references herein) on Miras in 
the LMC. Their results, which are summarized 
below, largely contribute to progress in the 
knowledge of the red-giant or supergiant phase 
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of evolution that ends by drastic events such 
as planetary nebula ejection or supernova ex- 
plosions. One must also recall that these obser- 
vations include the brightest and, hence, most 
massive stars, some considerably more massive 
than the intermediate-mass stars that form the 
bulk of the M, S, and C stars. 

Evolutionary calculations show that the 
LPV’s fall into two distinct groups: (1) non- 
degenerate core helium (or carbon) burning 
supergiants with mass M = 9 M 0 on their first 
appearance as red giants; (2) asymptotic giant 
branch (AGB) stars (among them Miras) with 
degenerate carbon/oxygen cores and hydrogen- 
and helium-burning shells, climbing up the 
giant branch for the second time. Relevant 
properties of the AGB stars must be noted: (1) 
theoretical maximum possible luminosity for an 
AGB star is M bol - -7.1; (2) their main- 
sequence masses are =9 M Q ; and (3) they are 
able to dredge up carbon and s-process elements 
to the surface through undergoing helium shell 
flashes, so that they participate in the enrich- 
ment of the interstellar medium in these ele- 
ments at the onset of current loss of their 
envelopes by stellar winds. 

Observational data of LPV’s in the Magel- 
lanic Clouds confirm the division of these stars 
in supergiants and AGB’s: (1) An absolute K- 
magnitude (M K ) against period (P) plot (see 
Figure 2 in Wood et al., 1983) shows that the 
stars fall into two sequences separated by the 
theoretical AGB limit in M K : below this limit, 
stars are AGB stars, whereas above the limit, 
they are supergiants. (2) In this plot, carbon (N) 
stars and S stars (rich in s-process element Zr) 
are located in AGB sequence only, as expected 
from theoretical evolution calculations for 
which C and s-process elements are dredged up 
to the surface during helium shell flashes suf- 
fered by AGB stars. (3) K light curves indicate 
that the AGB stars have pulsation amplitudes 
of 0.5 to 1 .0 mag, whereas the supergiants have 
smaller K amplitudes (<0.25 mag). (4) Finally, 
the mean J-K colors are rather similar for the 
supergiants whatever the period, while they 
become redder with period for the AGB stars. 


Wood et al. (1983) discuss the evolution and 
physical properties of the LPV’s in the Magel- 
lanic Clouds through the (< M bol >, ^-dia- 
gram (Figure 1-3) on which the two regions oc- 
cupied by the supergiants and the AGB stars 
are again clearly delimited (by dotted lines on 
the figure). For each bolometric luminosity, the 
core mass, M , of the AGB stars is indicated 
(theoretically, a linear relation exists between 
these quantities), as well as the main-sequence 
mass, M m , for the supergiants. (These quan- 
tities are related regardless of subsequent mass 
loss.) Theoretical lines of constant mass (the 
pulsation mass) are also plotted assuming the 
LPV’s are first-overtone pulsators (see details 
in the section Modes of Pulsation of the Long- 
Period Variables). From Figure 1-3, Wood et 
al. note that the pulsation masses: (1) for the 
supergiants, range from ~7 to -30 M Q , with 
a strong concentration from - 10 to -25 M Q ; 
(2) for the AGB oxygen-rich stars, tend to range 
from -0.7 to 7 M Q ; and (3) for the carbon 
stars in their sample (mainly in the field of NGC 
371 and, hence, perhaps coeval), tend to con- 
centrate from 1.0 to 1.5 M Q . The evolutionary 
implications of Figure 1-3 are as follows. 

The supergiants (including some galactic 
supergiants) describe a continuous band from 
(P, M bo/ ) - (100 days, -6.5) to (850 days, 
-8.2), with all but three of the Magellanic 
Cloud objects having P > 400 days. On the 
other hand, the initial main-sequence masses in 
this range of luminosities and P > 400 days go 
from —16 to 23 M Q . The pulsation masses, 
which are present-day masses, indicate that the 
supergiants have lost up to half their mass since 
the main-sequence phase: the more luminous 
and the initially more massive the supergiants, 
the more mass they have lost in previous evo- 
lutionary phases. The supergiant evolution re- 
sults from: (1) the decrease in mass, causing P 
to increase (P a M" 1/2 ), and (2) the change in 
radius with evolution (P a R 3/2 ). 

As AGB stars evolve, they increase their 
luminosity with time, and correspondingly, they 
increase their periods, P. It is worth noting that, 
for the first time, stars on the upper AGB (i.e., 
with luminosity right up to the AGB limit 
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(M bol ~ 7.1)) have been identified. The arrows 
on the constant mass lines for AGB’s in Figure 
1-3 indicate that these lines are approximate 
evolutionary tracks. (Mass-loss effects are 
shown to be qualitatively unimportant.) To 
each luminosity corresponds an observed maxi- 
mum period, the approximate position of which 
is given on the dotted line. As discussed by 
Wood et al. (1983), the normal evolution of an 
AGB variable evolving on an evolutionary track 
with M < 3.5 M Q stops when its track crosses 
the dotted line (maximum period reached); at 
this moment, the total mass of the star is still 
greater than the core mass, showing that the 
envelope hydrogen and helium has not yet been 
fully destroyed by nuclear-burning shells. The 
authors suggest that the evolution is terminated 
by the loss of the stellar envelope, followed by 
the formation of a planetary nebula, with 
nebula mass they estimate to be from ~0.1 
M Q (for the oldest evolved stars) to 2.1 M Q 
(for stars of initial mass 5 M Q ). The mecha- 
nism they suggest for ejection of the envelope 
is the switch in mode of pulsation from first 
overtone to fundamental. This suggestion is 
supported by theoretical pulsation models, in- 
dicating that, as an AGB star increases its 
luminosity, it pulsates in modes of lower and 
lower order until the fundamental one- 
thought we have to keep in mind that the mode 
of pulsation remains a disputed question (see 
the section Modes of Pulsation of Long-Period 
Variables ). 

The galactic OH/IR stars (see the section 
Classification) with P ~ 1000 days are a good 
example of these stars pulsating in the funda- 
mental mode and being in a transition phase 
between LPV and planetary nebula. When the 
evolution track of a more massive AGB star 
with M > 3.5 M Q reaches the AGB maximum 
luminosity limit, such a massive star will pro- 
duce a supernova by the ignition of carbon in 
the degenerate core. 

Finally, Wood et al. (1983) discuss further 
important points related to nucleosynthesis on 
the AGB. During helium shell flashes, 12 C and 
s-process elements are dredged up to the stellar 
surface. Carbon stars have dredged up suffi- 


cient carbon to be characterized by C/O > 1 . 
The S stars, located on the AGB between M and 
C stars, have not yet dredged up sufficient car- 
bon to have C/O > 1 , but they have dredged 
up sufficient Zr, an s-process element, so that 
ZrO bands are observed to be enhanced in these 
stars. Carbon stars in the Magellanic Clouds 
were already known to have their higher lu- 
minosity at M bol 6. This is confirmed by 

Figure 1-3, in which all the more luminous 
AGB stars are noncarbon stars, among them 
S stars. It is clear that s-process elements are 
dredged up in these upper AGB stars, but C/O 
remains less than unity. An explanation favored 
by Wood et al. to the fact that some AGB stars 
of type S are found above the most luminous 
C stars is that the 12 C dredged up is converted 
to 14 N during quiescent evolution between 
shell flashes by CNO cycling at the base of the 
envelope convection zone; therefore, the AGB 
stars are the source of primary nitrogen (since 
14 N is synthesized from the hydrogen and he- 
lium in the star at its birth), as well as s-process 
and carbon elements. 

VARIABILITY TYPES OF 
GIANTS AND SUPERGIANTS: 

A SURVEY OF OBSERVATIONS 

AND INTERPRETATIONS 

The brightness fluctuations of these intrinsic 
or physical variable stars are caused by geomet- 
rical and physical factors: fluctuations in di- 
ameter, temperature, pressure, molecular and 
dust opacity strengths, etc. Apparently, com- 
plex phenomena make up their entire unstable 
atmosphere: shock waves, strong coupling of 
convection and pulsation, and certainly mo- 
tions of clouds of matter from the photosphere 
to the circumstellar layers (if any) and infall- 
ing of this matter at other phases of the varia- 
tion, perturbation in the layers due to compan- 
ion tidal effect. 

The visual light curves tend to be regular, 
with large amplitudes for the Miras; they tend 
to be less regular and/or erratic, with much 
smaller amplitudes for the semiregular and ir- 
regular variables. Among all these variables, we 
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P (days) 

Figure 1-3 . Bolometric luminosity (M bol ) plotted against period (P) for long-period variables (LPV y s) 
in the Magellanic Clouds (SMC and LMC ). Also included are some galactic super giant variables. Stars 
labeled t( LMC Miras” are LPV’s studied by Glass and Feast (1982). Dotted lines delineate the regions 
occupied by supergiants and AG B stars. Continuous lines are lines of constant mass assuming that 
LPV’s are first-overtone pulsators. Also indicated are theoretical values of the core masses (M J of 
AGB stars and the initial main-sequence masses (M m ) of supergiants. See text for comments (from 
Wood et al.y 1983). 


find stars with all different values of the ratio 
C/O ( < 1 for M-type stars, <1 for S-type stars, 
and >1 for C stars). They may be either giants 
or supergiants (Kukarkin et al., 1958). All these 
stars can also be identified as pulsating or erup- 
tive variable stars (Strohmeier, 1972). 

In the pulsating group , we recognize three 
types — Miras, semiregulars (SR), and irregulars 
(L) — as follows: 

1 . The Mira Ceti-type stars are long-period 
stars with visible amplitude variation 
over 2 magnitudes (and for some cases 
up to 5 and larger), with well-expressed 
periodicity, good regularity, and periods 
from 2 months to 1 year or more. (They 
belong to the AGB long-period variables 
(LPV); see Figure 1-3.) Abound the max- 


imum, characteristic emission spectra are 
observed, mainly the hydrogen Balmer 
series. They occupy the tip of the AGB. 
Representative stars are o Cet, R Lep, 
X Cyg, etc. 

2. The semiregular variables are further 
classified SRa, SRb, and SRc. Their visi- 
ble amplitude variations are smaller than 
those of the Miras. The SRa light curves 
are not regular and have strong varia- 
tions from one period to another. Dif- 
ferent durations of cycles with ir- 
regularities and some constancy of 
brightness are observed in the SRb light 
curves. The SRa and SRb variables are 
giants; the SRc’s are supergiants with 
SRb behavior. Some of the giants are 
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Figure 1-4. 2. 7-fim light curves: Top: S Ori is a Mira-type variable with a period of 419 days; bottom: 
RTHya is an SRa-type variable; center: the irregular (Lb) variable V352 Ori shows a small flux varia- 
tion on a time of 2 years (from Maran et al. f 1977). 


probably on the AGB, others on the first 
ascent of the red-giant branch (RGB). 
Typical representatives are WZ Cas, q 
Per, S Per, etc. 

3. The true irregular variables are called Lb 
or Lc. The former are giants; among 
them, 0 Peg, CO Cyg, and V Aql are the 
typical representatives. The latter are su- 
pergiants; the best known is TZ Cas. 

Figure 1-4 shows the three characteristic 
curves of Mira, SRa, and Lb stars (Maran et 
aL, 1977). 

R CrB is the typical representative of the 
RCB stars, the eruptive group . These stars are 
of spectral classes, F, G, K, and R (hot carbon 
stars) and have a very high luminosity, as 
described in the previous section. They are 
characterized by slow nonperiodic drops in 
brightness of different amplitudes (1 to 9 mag- 
nitudes) and of different duration (up to several 


hundred days for some). The increase in bright- 
ness is generally rapid (a few days). Besides R 
CrB, other representative stars are XX Cam, 
RY Sgr, and SU Tau. The star, q Cas, is one 
of these stars which appear in the literature as 
pulsating (SRc) or eruptive (RCB) stars. 

Let us remark that we do not intend to study 
here certain peculiar but rare types of variables, 
mainly among the carbon spectral-type stars 
such as Isb (V374 Aql, UX Cas, and DY Aur), 
RVa (AC Her, etc.), or the cepheids with car- 
bon molecular bands (RU Cam, etc.). These 
types of variables are described by various 
authors (e.g., Payne-Gaposchkin, 1954a). 

In the late-type star family, the carbon stars 
are not only different by the C/O ratio, which 
is greater than one, but also by the behavior 
of their light curves (i.e., perhaps by specific 
reactions of their atmospheres to mechanical 
energy transport). First of all, the distribution 
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Table 1-5 

Percentage of Variability Types 
in M, S, and C Stars* 


Spectral Type 

Type of Variability! % ) 
Mira SR Lb 

M 

46 

29 

25 

S 

58 

23 

19 

C 

20 

38 

42 


*From Alksne and Ikaunieks (1971). 


of the different types of variability is very dif- 
ferent for the three spectral classes, M, S, and 
C. Table 1-5 (Alksne and Ikaunieks, 1971) 
shows that the number of semiregular and ir- 
regular stars increases sharply as one goes to 
carbon stars. This might be partially due to a 
selection effect, since these stars are faint and 
their light variations are not so easily followed 
as in the M stars. Half of the variable stars of 
M and S classes are Mira variables. For carbon 
stars, only 20 percent are Mira variables, 40 
percent are SR, and 30 percent are Lb. The 
Mira and semiregular carbon stars have the 
largest mean periods; carbon Miras have the 
smallest photographic amplitude variations 
among all the Mira variables, whereas the ir- 
regular carbon stars have the largest ones 
among all the Lb variables. 

The study of the cool star variability is at 
present based on visible and IR photometric 
data : Far-UV and far-IR data are too rare (see 
the section on Photometric Observations) to 
study the photometric variability. The present 
data are borrowed from a large number of 
papers, and consequently, they suffer from too 
many inhomogeneities like observations at dif- 
ferent phases, use of different photoelectric 
systems ( UBV , DDO y etc.), and data reduction 
method not always being adapted to cool stars. 
To achieve homogeneous data covering several 
periods, Bouchet (1984b) carried out a 
photometric analysis (from the blue to far- 
infrared ranges) of a sample of galactic carbon 


stars at ESO-La Silla during 4 years. This work 
is an example of a long set of data for Mira, 
SR, and L carbon stars made under the same 
observational conditions. It will be a good sam- 
ple for testing the evolution and/or dynamical 
pulsation theories. Another extensive example 
concerning infrared photometry on M, S, and 
C stars comes from Catchpole et al. (1979). 

Before describing the general characteristics 
of each group of variables, we focus on two im- 
portant remarks: 

1. The Mira variables, as well as the SRb, 
Lb, or HdC stars, have microvariations 
with a duration much smaller than the 
length of the period or the pseudoperiod 
generally known (some hours or days 
compared to several months or years). 

2 . In the HR diagram, some stages of the 
evolution of the variables can be followed 
during a human life. For example, some 
of these stars become apparently planetary 
nebulae in less than 10 years (e.g., HD 
59643, HM Sge, or V1016 Cyg). (See the 
section Irreversible Changes in SR and L 
Stars.) 

The Miras 

By “Miras,” we distinguish a type of light 
curve, not a spectral type of stars. Many papers 
have confused the term and sometimes speak 
about Miras without any definition of the spec- 
tral type of the stars; this can occasionally be 
confusing for a nonspecialist. Therefore, we 
have to be clear: we point out that the Mira 
variable-type stars exist among the three spec- 
tral types: M (R Aqr, R Aql, R Cae, R Cen, 
o Cet, R Hya, U Ori, etc.), S (R And, R Cyg, 
X Cyg, etc.), and C (U Cyg, V Cyg, T Dra, R 
Lep, SS Vir, etc.). Moreover, the chemical 
equilibrium of the elements inside the atmo- 
sphere of these later stars shows that it is im- 
possible to find OH in carbon stars. Conse- 
quently, the papers on OH masers usually refer 
to the spectral-type M stars. This is meaningful 
when looking for eventual correlations among 
the observed features. 
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The Miras are the most observed and best 
known among the intrinsic variable stars. The 
earlier curves of light were recorded during the 
18th century: 

During the autumn of 1596, the shepherd 
David Fabricius discovered a variable star 
in the constellation of the Whale, a star of 
the 3rd magnitude. In 1603, the astronomer 
Bayer added this star, which was of the 4th 
magnitude, onto his famous atlas and la- 
beled it with the letter o. In 1635, Dutch 
astronomers Holwarda and Fullenius like- 
wise showed that this star is variable, and 
in 1638, Holwarda gave the amplitude of 
variation to be from the 2nd to the 10th 
magnitude. Hevelius observed this star dur- 
ing 50 consecutive days and determined its 
period: 331 days. At that time, this star was 
the sole variable known in the sky. It was 
therefore named “Mira Ceti” or the “magic 
star of the Baleen constellation” (Jacchia, 
1933). In 1839, Argelander, the German 
astronomer who is often called the father of 
variable star astronomy, observed a max- 
imum of 2.2. 

X Cyg was discovered by G. Kirch in 
1681, R Hya by G.P. Maraldi in 1704, and 
R Leo by Koch in 1782. In the following 
century, Baxendell discovered, among many 
stars, S Aql (1863) and U Boo (1880); 
Pogson discovered R Cyg in 1852 and R 
Cas, R UMa, and S UMa in the following 
year. R Car was discovered at Cordoba in 
1871. 

Let us quickly describe the famous work 
of Padre Angelo Secchi. In 1868, he observed 
nearly 100 red stars with the 9-inch Merz 
refractor of the Collegio Romano and a 
small direct-vision spectroscope. Among 
them, he included 15 long-period variables 
from Schjellerup’s catalog (1866). The first 
slit spectrogram of LPV stars could indeed 
be this of o Cet taken by J. Wilsing with the 
Potsdam refractor in February 1896. Vogel 
(1896) remarked on the great intensity of the 
bright hydrogen lines and the absence of He. 
In 1897, the Reverend Walter Sidgreaves 
(1897) pointed out the rapid variation of o 


Cet during a week, “but a marked change 
in the relative intensities of the yellow-green 
and the blue radiation appears to have taken 
place during the cloudy week between 
December 2 and December 11.” He 
confirmed the absence of He noted by Vogel 
the previous year: “Of the hydrogen lines 
He is still absent, lost, or much weakened 
in the calcium absorption.” The notion of 
stellar envelope is already contained here, 
100 years ago.... 

Nearly all photoelectric V magnitudes’ of o 
Cet and T Cas at bright phases lie above the 
curves given by the American Association of 
Variable Star Observers (AAVSO) and may in- 
dicate a systematic difference between the 
AAVSO and the V magnitude (Lockwood and 
Wing, 1971). Stanton (1983) gives an experi- 
mental relation between the visual magnitude 
(m ) and the photometric measurements V and 
B-V. 

The Mira light variations with phase are 
sometimes very large in the visible. For exam- 
ple, the Mira, x Cyg (spectral type S), has a 
visual amplitude variation that is about 9 mag- 
nitudes (4000 times brighter at maximum light 
than at minimum). Consequently, UV and vis- 
ual spectrograms can hardly be obtained, espe- 
cially near the light minimum of the Mira stars 
(Maehara, 1971). The Mira stars were analyzed 
mainly by the Japanese group (Maehara, 1968; 
Yamashita and Maehara, 1977, 1978; Yama- 
shita et al., 1978; and others). 

Most of the Mira changes are cyclic, with 
a period equal to that of the visual light curve. 
Some obvious inconsistencies appear in the ob- 
servational data, as well as phase lags between 
light curves in different colors, which are 
caused by the actions of the cyclical tempera- 
ture variations on the various layers of their 
very extended atmospheres. Joshi et al. (1980) 
give the evolution of the effective temperature 
during 12 days for o Cet: on October 6, 1979, 
Mira was an M5.5 III star with 2470 K, and on 
October 18, 1979, it reached M4.5 III with 2760 
K. Because the period of individual light cycles 
sometimes differs from the star’s mean period 
by several percent, the observable properties are 
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not repeated exactly from cycle-to-cycle 
(Barnes, 1973; Wing, 1980). 

Fluxes of some M Miras show a significant 
deviation from a stellar blackbody beyond 8 pm 
due to the presence of a cold thin dust shell 
(\0-pm excess radiation), whereas some have a 
great excess of infrared (see M. Querci, this 
volume) due to the presence of a thick cir- 
cumstellar dust envelope (Epchtein et al., 1980). 
Among the M Miras with the largest infrared 
excess at 10 pm y we know: IRC + 10011, IRC 
+ 50137, and IRC -10529. Some of them have 
masers; the OH satellite line at 1612 MHz is 
seen in many sources, representatives of which 
are IRC + 10011, WX Sgr, R Aql, NML Tau, 
IRC + 50137, VY CMa, and NML Cyg. U Her 
and U Ori are known to be 1665-MHz emitters. 
The M star, RV Hya, has no OH microwave 
emission (Zuckerman, 1980). For the first time 
in a carbon-star envelope, Henkel et al. (1983) 
found a maser line in IRC + 10216. Like OH, 
the observed SiS maser has two horn features. 

Some changes from OH type I to OH type 
II Miras were observed during recent years, 
such as on R Leo, for which Lesqueren (1983) 
reports an unusual intensity variation: the flux 
at 1667 MHz has steadily decreased from 
L5E-22 W/m 2 , reached during the 1980 cycle, 
to a smaller value 0.1E-22 W/m 2 at the max- 
imum of the 1982 cycle. Another striking ex- 
ample is U Ori. Pataki and Kolena (1974a, 
1974b) found a bright 1612-MHz flash up to 
22 Jy and a disappearance of the 1667-MHz line 
below the 0.4-Jy detection level on May 28, 
1974 (the 1612-MHz line is characteristic of the 
OH type II maser star), whereas before 1974, 
particularly on July 28 and November 26, 1973, 
only the main lines (at 1667 and 1665 MHz, 
characteristic lines of OH type I maser stars) 
are present in emission. Moreover, on May 28, 
1974, no radical change is registered in the out- 
put of the 1665-MHz main line. Coincidental- 
ly, the 43.122-GHz SiO maser line was not 
detected on June 2, 1974, by Snyder and Buhl 
(1975), but it was detected exactly two U Ori 
periods later by Balister et al. (1977) with an 
energy level of 89 Jy. On the basis of these 
observations, Cimerman (1979) develops a 


monitoring that shows no variation on the 
1667-MHz main line until May 1977. On 
January 25, 1978, this line appears with 0.8 Jy 
and -42 km/s. The 1665-MHz line was not 
observed as often as the 1612-MHz line because 
of its weakness. On March 1, 1978, the 1612- 
MHz line has a flux peak of 2.2 Jy, which is 
about one-half of that observed by Pataki and 
Kolena (1974b) in May 1974. On July 15, 1975, 
the spectrum of 1612 MHz consists of a sharp 
feature at -42 km/s and a group of blended 
lines peaking at -46 km/s. This large blend 
evolves to a sharp feature at -47 km/s, and the 
-42-km/s peak remains constant (spectrum of 
January 25, 1978). Besides the variation of the 
line profile, a variation on the intensity of the 
1612-MHz line occurs: the average level of 
45.E-22 W/m 2 detected by Pataki and Kolena 
(1974a, 1974b) remains constant until the begin- 
ning of 1976, when it drops to one-half its 
previous level in 100 days and continues to 
decline (see correlations with visible variations 
below). This decline is confirmed by Jewell et 
al. (1979; see their Figure 15). The main-line 
emission has exhibited highly circular polarized 
components and large semiperiodic variations 
(Fix, 1979). Cimerman (1979) attributes the 
origin of these events to the central star; this 
assumption is corroborated by the coincidence 
in velocity between one of the 1612-MHz lines 
and the optical emission lines. All these lines 
are interpreted as the outcome of a shock or 
magnetohydrodynamic (MHD) wave originat- 
ing under the photosphere of U Ori (Merrill, 
1960; Wallerstein, 1973). The model of Elitzur 
et al. (1976), which explains the 1612-MHz 
maser in sources with a thick shell like IRC 
+ 10011, seems to be inadequate for U Ori , 
where the dust shell is very thin. In this case, 
another pumping process for the 1612-MHz line 
may exist in addition to the one that uses the 
dust grains (Cimerman, 1979). A mapping of 
the U Ori neighborhood was made by the VLBI 
on the principal and satellite frequencies. Un- 
fortunately, the variations on 1665-MHz lines 
(Fix et al., 1980) were made 3 years after the 
1612-MHz ones (Reid et al., 1977), and com- 
parisons are therefore hazardous. 
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Phase lag. Since their discovery, the Miras 
have attracted generations of astronomers until 
the present day; the pioneers are Herschel, Sec- 
chi, Pettit and Nicholson, Hetzler, Hoffleit, 
Campbell, Stern, Joy, Merrill, and others. Oc- 
casionally, astronomers of more recent times 
have made the same discoveries as those of their 
forerunners for want of looking at their old 
publications. The study of the phase delay be- 
tween the visible, photographic, and near- 
infrared light curves illustrates this fact to 
perfection. 

Traditionally, the epoch of zero phase refers 
to the visual maximum. The phase lag (or shift 
or delay) of the maximum of an infrared or ra- 
dio curve is the epoch of this maximum, minus 
the epoch of the visible maximum coming im- 
mediately before it. The first correlations be- 
tween the variations at different wavelengths 
were made by Pettit and Nicholson (1933). As 
an example, they have shown that o Cet, R 
LMi, R Leo, R Hya, R Aql, and \ Cyg have 
a real maximum of energy shifted by a phase 
lag of about 0.14 after the visual maximum. 
Hetzler (1936) has monitored 30 Mira variables 
with photographic plates at an effective wave- 
length of 8500 A (AX = 100 A). His light curves 
have a much smaller amplitude than the cor- 
responding visual light curves because of the 
longer wavelength effect and the decreased 
blanketing by TiO bands. Hetzler showed that 
the maxima of the infrared and the visual 
curves are reached at about the same time, but 
the infrared maximum persists longer. R Lyn 
and R Cam are two characteristic examples of 
simultaneity. R Hya is an exception; the in- 
frared maximum precedes the visual one. In this 
case, the IR phase lag is negative. This 
phenomenon is to be confirmed by careful 
monitoring of many Mira-type stars at several 
wavelengths. Generally speaking, the infrared 
light maximum typically lags the visual by 0.1 
to 0.2 period, agreeing with all the previous 
quoted observations and with those of Men- 
doza (1967), Lockwood and Wing (1971), and 
Frogel (1971). 

An analysis of the spectrum of the spectral 
M class Miras allowed the astronomers of the 


late 1960’s to distinguish the behavior of the 
continuum from that of the absorption lines. 
With the 1104 filter, Wing (1967a) showed that 
the I(1.04-/*m) magnitude of x Cyg was brighter 
at phase + 49 days than at 0 day by 0. 1 magni- 
tude; this is very similar to the phase lag noted 
in radiometric curves obtained with the vacuum 
thermocouple by Pettit and Nicholson (1933). 
He explains this phase lag by a two-layer model. 
The infrared magnitude measured at a contin- 
uum point refers to the photosphere, while the 
heavily blanketed visual magnitudes presuma- 
bly refer to the upper layers of the atmosphere. 
He concluded that the temperature variations 
of the two regions appear to be out of phase. 
At present, such an explanation must be viewed 
in the context of shock-wave interpretation (see 
M. Querci, this volume). 

Fillit et al. (1977) have analyzed the time 
variation of some type I OH Mira sources and 
the correlation with their visual magnitude. 
More than 2-year monitoring of seven OH/IR 
stars was done on six Miras (R LMi, R Cas, Y 
Cas, RS Vir, U Her, and S CrB) and on an SRc 
variable (S Per). Even allowing for experi- 
mental uncertainties, a systematic phase lag of 
about 20 days between OH lines and visible 
light is observed on S CrB (Figure 1-5). These 
authors implicitly assume that the shape of the 
radio and visible curves are the same. This as- 
sumption must be proved. 

The phase delay of the [V = 1, J = 2-1] 
SiO curve on R Leo, is greater than that of the 
2.1-fim curve, whereas the 2.7-/4m and SiO 
curves have the same lag as on o Cet. The SiO 
masers observed by Spencer et al. (1981) are 
interpreted as collisional pumping by Bujar- 
rabal and Nguyen-Q-Rieu (1981). This requires 
high temperatures and densities in the SiO 
layers (T K ~ 2000 K, n H ~ 10 10 cm -3 ), and 
consequently these layers are close to the star 
(R ~ 1.5 /? + ). Hinkle et al. (1976) indicate that 
the SiO strength in Miras varies markedly with 
phase (see also Nguyen-Q-Rieu, this volume). 

The behavior of another maser line, H 2 0 at 
1.35 cm, and of lines at 2.2 /un and in the visible 
were compared (Schwartz et al., 1974) for some 
Mira stars (U Her, S CrB, U Ori, R Aql, and 
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Figure 1-5. The flux and the magnitude for each of the OH maser lines: 1612 MHz, 1665 MHz, and 
1667 MHz lines of S CrB were fitted by least squares using the visual light curve (from Fillit et al., 1977). 


W Hya; see their Figures 1 through 5). The 
mean curves show the phase lags of radio and 
infrared curves relative to visual ones. These 
lags are characteristic on U Her and S CrB. 
Some stars do not have the usual phase lags 
(e.g., the H 2 0 maser curve of U Ori has a 
large lag compared to the infrared curve at 2.2 
tim). Although the IR curve of R Aql has a 
phase lag of 0.2 compared to the visual curve, 
the 1.35-cm curve has the same phase as the 
visual one. W Hya shows an analogous be- 
havior. This nonconventional behavior requires 
comment. Is it real? Yes, if the observed values 
at the different wavelengths are obtained simul- 
taneously; no, if there is a time lag between 
them: in this case, we find the presence of large 
amplitude variations during small intervals of 
time (a small percentage of the period). If the 
quoted values are from different cycles, another 


interpretation is possible: the behavior of the 
maser lines with the phase is not exactly the 
same in all the cycles. This was already dem- 
onstrated by Wing (1967b). Cox and Parker 
(1979) demonstrate that the H 2 0 was often not 
stable for more than a few cycles of stellar 
brightening. (A more comprehensive descrip- 
tion of the nature of the different maser pumps 
is made by Nguyen-Q-Rieu in this volume.) 

From 1972 onward, observations of some 
M Miras were obtained by using the satellites 
of the U.S. Air Force (Maran et al., 1977). Be- 
cause the spectral interval on the water-vapor 
bands around 2.7 ^tm was chosen, correlations 
with the H 2 0 maser should be possible with- 
out the perturbations introduced by the variable 
telluric water vapor. Figure 1-6 shows the 
2.7-/im infrared and the visual curves of two 
Miras, S CMi and o Cet. 
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Figure 1-6. 2. 7-/t infrared curve (top) and AA VSO visual curve (bottom) for two Mira variables, S 
CMi and o Cet. The left ordinates are the infrared fluxes and the right ones are for the visual magnitude 
(from Maran et at., 1977). 


Amplitude of Light Variations at Different 
Wavelengths. The S Mira, x Cyg, is a very good 
example for this purpose. Maehara (1971) ob- 
served a visual amplitude variation of about 9 
magnitudes, whereas the photographic infra- 
red variation was only about 2 magnitudes. 
Although similar data for the ultraviolet region 
is incomplete, some information on the 
variability of the energy distribution of x Cyg 
from 2500 to 3350 A at two different phases 
is shown on Figure 1-7. The flux longward of 
about 3050 A becomes considerably fainter at 
phase 0.18 than at the light maximum (phase 
0.04). The total energy emitted in the range 
from 2000 to 3240 A is 4.5 E-ll erg/cm 2 s at 
phase 0.04 and is lower by a factor of 3.53 at 
phase 0.18. The stellar continuum is even more 
depressed at phase 0.22 on a high-resolution 
spectrum. This remarkable change is mainly 
due to a change in the temperature of the region 
of continuum-formation (Cassatella et al., 
1980). 

Lockwood and Wing (1971) show examples 
of some M Miras for which the curves of the 
filters, V, 78, 87, and 88, are depressed by TiO, 


and the filter- 105 curve by VO. These curves 
are strongly correlated with the differential 
temperature effect, the local continuum, and 
the variation of diameter of the forming layers, 
so that these curves are the best reflection of 
the physical parameter variations of the at- 
mosphere along the period. Figure 1-8 describes 
the curves of two Miras, o Cet and T Cas, in 
which the minima appear to occur at the same 
phase. The behavior of R Aql is different from 
that of o Cet and T Cas because the visual 
minima of these two stars have a phase lag 
relative to the infrared minima. This is shown 
on the 2.7-^m curve obtained with the U.S. Air 
Force satellite (Maran et al., 1977, 1980). When 
one considers the infrared spectral range, one 
sees that the amplitudes of variation 
[F /F . 1 decrease with increasing 

1 max min J 0 

wavelengths from 1 .2 to 3.5 jtm and are roughly 
constant between 3.5 and 10 /im. The 
amplitudes of the 1612-MHz OH masers which 
are linked with the IR fluxes are significantly 
less than, or at most equal to, the infrared 
amplitudes (Harvey et al., 1974). Herman 
(1983) followed OH/IR Miras during a period 
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Figure 1-7. Ultraviolet energy distribution of the S Mira variable x Cyg from 2500 to 3350 A at phases 
0.04 (thick line) and 0.18 (thin line). The total energy flux in this spectral region is less at the latter 
phase by a factor 3.53 (see text) (from Cassatella et al., 1980). 


of 3 to 5 years with the 25-m Dwingeloo ra- 
diotelescope in order to study their varia- 
bility. 

After considering the amplitude variations 
in visible, infrared, radio, and ultraviolet wave- 
lengths, we now follow the amplitude variation 
at one wavelength and point out the character- 
istic behavior of the Miras. Maran et al. (1980) 
show five well-observed minima of S CMi (with 
the U.S. Air Force satellite), where we note the 
agreement of flux density level at all five 
minima. The flux of S CMi is remarkably repre- 


sentative of all minima. The authors emphasize 
that the minimum is the “normal state of a 
long-variable star.” However, there has been 
a good deal of controversy about the reproduc- 
tion of Mira minima. Glasby (1968) said that 
o Cet has variable minima: “on occasion [the 
minimum] has only been as faint as magnitude 
8.0 while at other times it has sunk as low as 
10.2 magnitudes.” Campbell and Payne (1930) 
show variation of magnitude of the minima of 
S Pav. On the other hand, we sometimes ob- 
serve a strong contrast between consecutive 
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Figure 1-8. Light curves of o Cet (Mira) at left and T Cas at right, in V and at five in- 
frared wavelengths. The wavelength designations to the right of the curves are defined 
in Lockwood and Wing’s (1971) Table 1. The observed cycles of the stars are numbered 
serially; the IR observations start with cycle 1 in 1965 and are coded as follows: 

-1 = + ,0 = 0,1= • , 2 = 0,3 = m,4= D,5 = A,6= t\, 7 = T , 8 = V , 

9 =+ , and 10 = 0 (from Lockwood and Wing, 1971). 

maxima of Mira stars, such as the two maxima ima tend to alternate, and bright maxima tend 

of R Aql in 1972-1973 and those of S CMi in to occur before the predicted date (Harrington, 

1973 and 1974 (Maran et al., 1977). In the late 1965; Keenan, 1966). This could be correlated 

18th century, the irregularity of the maximum with the intensity of the pulse (i.e., the size of 

light was already noticed by William Hershel the energy reservoirs in the cycle). RS Cen, R 

on o Cet, where it is extremely well pronounced Cet, V Cen, T Gru, R Vir, RT Cen, and SS Her 

in the visual curve. In 1779, it attained the first are Miras with alternative bright and faint 

magnitude, being equal to Aldebaran in bright- maxima. 

ness. For many Miras, the cycle-to-cycle dif- Different Profiles of Light Curves. The 

ferences are not erratic; bright and faint max- AAVSO report 38 (1983) gives many visual light 
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curves of Mira variables. We warmly congratu- 
late all the AAVSO members for their arduous 
and comprehensive work; we are convinced of 
its usefulness to professional astronomers. 

The light curves of carbon stars are more 
gradual and more symmetrical than those of M 
class stars (Figure 37 in Alksne and Ikaunieks, 
1971). Campbell (1955) shows that the S star 
light curves are very similar to those of M stars. 
As a rule in M and S spectral types, the rise to 
maximum is somewhat more rapid than the de- 
cline to minimum. In general, D - 0.44, with 
D = (epoch of maximum-epoch of minimum)/ 
period. However, some exceptions exist, such 
as: o Cet (D = 0.62) and R Cyg (D = 0.7). 
(Because D is always larger than zero, it rep- 
resents the fraction of the period needed to go 
from the minimum to the maximum.) 

All the Mira light curves do not have the 
same shape (Figure 1-9). They can present: 

1 . A rapid rise to the maximum and a slow 
decline to the minimum: the majority of 
the Mira light curves have this shape, 
with D - 0.45 (R Gem, U Her, and U 
Ori). 

2. A slow and long rise to the maximum 
with a rapid decline with D — 0.7 (S CMi 
and U CMi). 

3. A narrow and sharp maximum with a 
large minimum (R Aqr, R Cas, Z Cas, 
U Cet, T Col, R Crv, R For, R Pav, T 
Sgr, RU Sgr, and R UMa). 

4. A large and round maximum and a nar- 
row and deep minimum (Z Cap, S Car, 
RT Cen, X Cet, T Hor, and S UMa). 
These stars have a period P < 225 days 
(Celis, 1977) so that the magnitude varies 
rapidly around the minimum. 

5. A sawtooth shape with sharp maxima 
and minima with D ~ 0.5 (RT Cyg and 
V Oph). The variations between their ex- 
tremes are rectilinear segments. These 
stars have periods shorter than 150 days 
(SS Her, W Pup, and R Vir; Celis, 1977). 


6. A smoothed sinusoidal shape (SS Vir, R 
Dra, W Cas, and U UMi). 

7. Two maxima of the same intensity and 
two minima with different depths per cy- 
cle (R Cen and R Nor, both discovered 
by Gould at Cordoba in 1871). The main 
maximum is followed by the main mini- 
mum. This curious light curve was repro- 
duced theoretically by Wood (1979) in an 
attempt to produce a larger post-shock 
velocity maximum in an isothermal pul- 
sation model. It appears that this situa- 
tion is apparently stable (over 75 cycles) 
where shocks with different velocities al- 
ternate. This visual curve also shows vari- 
ation from cycle to cycle in the rate of 
rise and fall, in the amplitude and in the 
height of the maxima, and in the depth 
of the minima (Marino et al., 1979). CR 
Mus shows a scatter of approximately 
one magnitude in the level of the V light 
curve; the shape of this maximum varies 
extensively from cycle to cycle. 

8. A hump on the rising branch, around the 
phase 0.7, so that the minimum is moved 
at earlier phase (R Lep, T Tuc, and S 
Vol). Such a hump is also observed by 
Hetzler (1936) on R Tri, R Vir, R Hya, 
T Cep, and x Cyg. For R Tri (Figure 
1-1 0a) and R Vir (Figure l-10b), the 
shoulder (or hump) is much more impor- 
tant in the photographic infrared (X = 
8500 A and AX = 100 A) than in the visi- 
ble, whereas for T Cep (Figure l-10c), the 
shoulder rises in proportion to the bright- 
ness variation in the 8500 A region and 
in the visible. To explain the 8500 A be- 
havior of R Tri and R Vir, Wing (1967b) 
suggests that the effective wavelength of 
Hetzler’s photographic infrared agrees 
well with his curves 87 and 88, which are 
likewise affected by moderately strong 
TiO absorption. For these stars, it 
therefore seems that the visible region is 
less affected by the TiO opacity. The 
humps on the rising branch of the light 
curve are generally centered at <t> = 0.7 
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Figure 7-70. Mira light curves illustrating 
humps on the rising branch (from Hetzler , 
1936), 


(e.g., T Cas in Figure 1-8). The majority 
of the stars studied by Lockwood and 
Wing (1971) show evidence of humps in 
at least a few cycles. In some of these 
stars, the humps have reoccurred (R Cam 
and T Cep), but in some stars, the humps 
are spread out during the phase interval 
0.6 to 0.8. Glasby (1968) indicates the 
presence of humps in the descending 
branch of the light curves. Vernon- 
Robinson (1929) observed secondary 
oscillation on the descending branch of 
XY Cas. The hump phenomenon is pre- 
sent in some Miras and not in others, and 
it comes and goes in some of them. The 
visual curve of x Cyg is a good example 
of a hump on the ascending branch 
(Figure 1-9). T Cas has a hump on the 
ascending branch of the light curve; its 
maxima and minima and the size of the 
hump vary from one cycle to the follow- 
ing one (Leung, 1971). The humps are 
considered as a rapid drop of the radius 


of the star (Lockwood and Wing, 1971). 

9. Inflection points in the increasing and de- 
creasing parts of the light curves (R Cha, 
T Hya, R Ind, V Mon, T Nor, and RT 
Sgr); their periods range from 200 to 350 
days (Celis, 1977). 

Nearly all the short-period stars (less than 
150 days) do not have very constant elements 
p and q with respect to the length of a given 
period P (Celis, 1977). As an example, T Gru 
(P = 136.6 days, Ml Iae-M2 Ibe) has P periods 
between 121 and 138 days, Q periods from 124 
to 145 days, and partial period p from 54 to 
83 days. The amplitudes go from 1 .80 to 3.05, 
and the light minima differ up to 1.5 mag. 

All these light curves might reflect the inter- 
action of the opacities and density variations 
of the different atmospheric layers and the 
propagation of shock waves (see M. Querci, 
this volume). 

Are light curves at other wavelengths similar 
to the visual ones? The two stars in Figure 1-8 
have V light variations of a strikingly different 
nature, which is preserved at the different ob- 
served wavelengths. The hump on the rising 
branch of the AAVSO curve for T Cas also ap- 
pears on the infrared curves, including the con- 
tinuum curve, 1104. 

Variation in Colors During Mira Light Cy- 
cles. All the observed Miras (for M, S, and C 
spectral types) show variations in color during 
their cycles of light, which are first indicated 
by the analysis of color indices obtained from 
visual observations (generally from AAVSO or 
Association Frangaise des Observateurs 
d’Etoiles Variables (AFOEV)) and broadband 
photometry U, B , K, 7, K , etc. Narrowband 
photometry with adapted filters for characteris- 
tic atomic (He I, Ha, and Na I) and molecular 
(TiO, ZrO, VO, etc.) features is proposed, at 
first by Wing (1967b). Spectrophotometry is 
also a good technique for investigating the vari- 
ation in colors during the Mira phases. We will 
not elaborate here on the variations of line in- 
tensity or profile; they are dealt with by M. 
Querci, this volume. 
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Figure 1-11. Mean color curves of seven Mira 
variables: S Pav, R Hor, T Gru, RS Sco, W 
Cen, R Oct, and S Oct. Ordinates and abscissas 
are color index and phase, respectively, all 
periods being reduced to the same horizontal 
scale. Maxima and minima are marked with 
short vertical lines and labeled “M” and “m” 
(from Campbell and Payne, 1930). 

As early as 1930, Cannon (1930) showed 
that the color index (difference between the 
photographic and photometric magnitudes) is 
at its greatest about 40 days before maximum 
and at its minimum for 100 to 150 days after 
maximum. Campbell and Payne (1930) ob- 
tained photographic light curves and used the 
AAVSO visual ones to give the color indices 
from observations assembled into 10-day 
means; these are obtained by subtracting the 
visual magnitude from the photographic one 
(Figure 1-1 1). Two opposite effects contribute 
to the photographic minus visual index: (1) the 
great strengthening of TiO absorption, which 
cuts down the red part of the spectrum and di- 
minishes the visual brightness more than it re- 
duces the photographic light; and (2) during the 
fall, the decrease of the temperature of the pho- 
tographic and visual forming layers which 
should strengthen the red end of the spectrum. 

During the Mira light variations, the B-V 
values vary and the spectral classes usually 



Figure 1-12. Color-color curves for luminosity 
class III nonvariable stars and Mira stars of 
spectral type M and S. Ordinates, in (I'-K) 
magnitudes; abscissas, in (V-I ') magnitudes. 
Heavy dots (plotted three times on shifted scales 
of ordinates) = nonvariable stars from K0 to 
M6; smaller dots = loops described during the 
cycle of variation for stars of average period 
239, 350, and 456 days, class Me; broken lines 
= Se stars, mean period 380 days; dotted lines 
extend the relation for nonvariable stars 
through points corresponding to the fall from 
maximum. On the lowest curve, points corre- 
sponding to K0 III, M4 III, M5 III, and M6 
III are identified. Circles show the relation be- 
tween maximal colors for Mira stars of classes 
M4e, M6e, and M7e (from Payne-Gaposchkin, 
1975). 

change, the earliest spectral class appearing 
near the visual maximum light. 

It is interesting to compare the color indices 
of Miras during their cycles with the indices of 
the nonvariable stars of the same spectral type 
and luminosity class. Payne-Gaposchkin (1975; 
Figure 1-12) shows color/color curves for Mira 
stars and nonvariable stars of luminosity class 
III. During the decline from the visual maxi- 
mum and a little later, the representative points 
of the Miras follow the relation shown by the 
nonvariable stars of the same spectral type. The 
I ' filter used in Figure 1-12 is related to the 
Johnson system by (I-K) = 0.745(1 '-K) - 0.13. 
Two remarks will be made about this figure: 
the length and the maximum width of the loop 
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(I -K), and consequently, its area versus (V-I ') 
seem to be correlated with the period/ampli- 
tude relation of the M Mira stars; secondly, 
during the maximum phase, the stars are bluer 
than during the minimum. This was also ob- 
tained by Mendoza (1967), using his multicolor 
photometry. The analyzed set of stars indicates 
that, with the light variation in V, the B- V in- 
dex and the other color indices change. Typ- 
ically, as the star becomes fainter in V (it is run- 
ning to the minimum), it becomes redder; how- 


ever, the C star, LW Cyg, is an exception. 

Spinrad and Wing (1969) demonstrate that 
the variation of the molecular (TiO + VO) in- 
dex around 1 fim is correlated to the observed 
temperature (Figure 1-13). On the other hand, 
the nonvariables and small amplitude variables 
are represented by their mean line running from 
near the origin to the box labeled RX Boo in 
Figure 1-13. (This box contains all five scans 
of the semiregular variable, RX Boo.) Each 
Mira moves on the color index/temperature 



Figure 1-13 . The relations between band strength and temperature are obtained from 
observations in the 1-fim region . The ordinate is the sum of TiO and VO indices (unit 
= 0.0/ mag) and is calibrated in terms of spectral type . The abscissa is 0 = 5040/T, 
where temperature T is obtained from a blackbody fit to the continuum . Observa- 
tions are connected in chronological order (from Spinrad and Wing , 1969). 


i 


43 


i 



plane along the period: 

1 . R Tri (A) is a relatively early-type Mira; 
its earliest spectral type was M3 , obtained 
at its 1965 maximum, but its correspond- 
ing temperature was that of a normal K4 
giant. 

2. x Cyg ( x ) was observed during a large 
part of two successive cycles, and its 
minimum approach paths were seen to 
be different in the two cycles. 

3 . o Cet ( • ) is shown from minimum to the 
post-maximum in the 1965 cycle. Its 
molecular index is much stronger than in 
non-Miras of the same temperature; con- 
sequently, this star is much hotter than 
the temperature obtained from its spec- 
tral type. 

Another extensive project is from Maehara 
and Yamashita (1978), who use a photoelectric 
scanning spectrophotometer. The wavelength 
region covered is 3700 to 5500 A with 13 A res- 
olution. They give the energy distribution mea- 
sured at flux peaks for a suitable sample of stars 
of luminosity class lab to III and spectral range 
from MO to M8 stars. Three S stars (V Cnc, R 
And, and W And) and one carbon star (SS Vir) 
are also included. The comparison of these 
curves with those of nonvariable M stars gives 
the spectral type variation of each Mira star 
during the period: o Cet ranges from M5 to M8, 
R Aur from M6 to M8, R Cnc from M6 to M8, 
etc. The S stars have an M-star behavior in the 
spectral region of this investigation because 
they suffer from TiO absorptions. The C-type 
Mira, SS Vir (C6,3e), shows an entirely dif- 
ferent distribution. The flux is too small to be 
measured shortward of X4250. The variations 
of the color indices in the spectral range 4250 
to 5000 A are very small compared to those of 
M and S stars; the variation of the C 2 band 
head at X4737 is also relatively small. 

Krempec (1975) shows changes of the cen- 
tral depth in spectral features, such as in C 2 
and CN molecular bands and in atomic D-lines, 
with the light phase. 


A comparison of the color changes observed 
on two prototype stars with two maxima per 
period (R Cen and BH Cru) shows differences: 
(1) the two U-B curves have reversed shapes, 
and (2) on the R Cen curve, the color index de- 
creases at the second maxima, while in BH Cru, 
a sharp increase is observed. The B-V curves 
of these two stars show only small differences. 

To conclude this study on the variations of 
color indices, we point out that there is a large 
variation in the color index for stars of the same 
spectral subclass; the strength of the hypotheti- 
cal shock wave that causes the variability should 
not be correlated with the spectral subclass. The 
same remark applies to the SRa variables ac- 
cording to Alksne and Ikaunieks (1971). In 
their review of carbon stars, these authors em- 
phasize that the mean values of B- V and V-I 
indices decrease somewhat in the sequence 
Mira, SR, L. 

Period Changes. Little is known about the 
frequency of period changes in Mira-type stars 
because a long-time base of observations is re- 
quired and the different shape and amplitude 
of successive cycles make the inferences more 
doubtful. Nevertheless, three kinds of period 
changes are listed: 

1 . A secular evolution by smooth increasing 
or decreasing of the length of the period 
or by a sine or cosine variation of the 
latter. 

2. Some shape changes of the light curve be- 
tween two or several consecutive periods. 

3. An abrupt change of period and an 
abrupt shift of epoch. 

From a 20-year survey, Hoffleit (1976) 
shows that log (A/VP) is spread around a linear 
function of log P, where A P is the maximum 
difference between the various periods of the 
same star and P is the mean period. This 
observed relationship holds for Miras, SR 
variables, W Vir stars, classical cepheids, and 
RR Lyr stars. 

Hoffleit (1979) has drawn a catalog of 356 
stars that have, at one time or another, been 
recorded as having changing periods; these stars 
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were selected in the General Catalogue of Vari- 
able Stars (Kukarkin, 1976), in the Geschicht 
und Literarur des Lichtwechsels (Prager and 
Schneller, 1934, 1963), and in some Informa- 
tion Bulletins of Variable Stars (IBVS). Among 
332 Miras, 273 individual increases in period 
(e.g., for V Cas and U Her) and 290 decreases 
(e.g., for T Cep) were noted. This catalog is 
available on request to its author. To have some 
idea as to the number of variable stars with dif- 
ferent kinds of changing period, we must refer 
to the 108 Sagittarius stars examined in suffi- 
cient detail by Hoffleit (1979). The observations 
could not be suitably represented by a constant 
period over 50 years, but the O-C plots for a 
constant period are best represented by one 
abrupt change. Three stars have the O-C curve 
represented by a parabola, indicating that the 
period is progressively increasing or decreasing. 
One star has its O-C curve represented by a 
sine-term, which suggests that the period is 
alternatively increasing and decreasing. 

The four stars, V462 Cyg, V734 Cyg, HO 
Lyr, and MV Sgr, have cyclical changes in pe- 
riod. The O-C curves are represented by sine 
or cosine terms (Hoffleit, 1979). U Boo and S 
Ser have sinusoidally changing periods (Sterne 
and Campbell, 1937), as well as six others (R 
And, RS Cen, Z Cyg, TU Cyg, W Her, and W 
Hya) given by Prager and Schneller (1934). 
However, this kind of ephemeris is a very 
simple mathematical representation of very 
complex physics, and also it breaks down dur- 
ing 50 years. Thus, several later catalogs 
dropped these representations in favor of listing 
discrete periods for successive time intervals. 

The conclusions of Nudjenko (1974) about 
the period changes are noticeably different 
from those of Hoffleit (1979). Forty-three 
(O-C) curves of the former demonstrate that 
the changes appear to be smoothly continuous 
rather than abrupt like these observed by Hof- 
fleit. Investigations of all of the data are re- 
quired. Table 1-6 summarizes the number of 
stars that have changed their period n times 
(Hoffleit, 1979). 

Sterne and Campbell (1937) and Hoffleit 
(1979) come to the same conclusion that the 


Table 1-6 

Number of Stars Having Changed Period* n Times 


n 

No. Stars 

n 

No. Stars 

2 

146 

7 

1 

3 

110 

8 

0 

4 

46 

9 

1 

5 

15 

10 

0 

6 

4 

IV 

o 

2 


* From Hoffleit (1979). 

stars, R Cnc, R Cen, and X Lib, have changed 
periods. 

Some examples of stars that have a pro- 
gressively decreasing period length are: 

1 . HS Aql, with its period ranging from 267 
days in the 1930’s, through 263.8 days in 
the late 1950’s, to approximately 260 
days in 1970’s (Thompson, 1981). This 
star might be climbing the asymptotic 
branch. 

2. W Tau, with its period ranging from 273 
days in 1887 to 253 days in 1956 (Schnel- 
ler, 1965). 

3. R Aql, with its regularly decreasing pe- 
riod (Hoffleit, 1979; Schneller, 1965) in 
a sinusoidal steady manner (Payne and 
Campbell, 1930). 

4. R Hya, with its period of about 500 days 
in the 18th century to about 390 days as 
noted by Merrill (1946). These period 
changes are used by Wood and Zarro 
(1981) to provide direct observational 
confirmation of the theory of helium 
flash shell (see their Figure 3). 

A list of such stars may be found in Hof- 
fleit’s catalog (1979). 

On the other hand, some stars, like V Cas 
and U Her, have increasing period length 
(Lockwood and Wing, 1971). Again, Hoffleit’s 
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catalog (1979) lists stars with likely smooth 
varying periods. 

It is useful to know that the prototype star, 
o Cet,was noted in Prager and Schneller (1934) 
to have 12 discrete periods ranging from 322.5 
to 335.4 days. For a total of 335 epochs of max- 
ima, the best constant period is 331.7 days, giv- 
ing a nonsystematic spread in O-C of 82 days 
and A(0-C)/P = 0.247. Consequently, Hof- 
fleit (1979) concludes that “this star is no longer 
considered as having a significantly changing 
period.” 

In 1936, Hetzler (1936) shows some small 
variations on consecutive light curves of Miras. 
The perturbation appears in the visible AAVSO 
light curves and in the infrared ones at X8500. 
The latter wavelength is affected moderately by 
strong TiO absorption (Wing, 1971). Figure 
l-10b shows two light curves (IR on the top and 
visual on the bottom) for two periods of R Vir 
separated by a period which is not indicated. 
The hump shown in the first period is not pres- 
ent in the second one (or is much less notice- 
able). 

R Car is an astonishing star; it shows a wide 
maximum cycle followed by a narrow maxi- 
mum one (Celis, 1977). 

Abrupt changes of period (rapid fluctua- 
tions of the length of period) appear in some 
Miras; they are visible from the differences 
(O-C) between the observed epoch of maxi- 
mum and the computed dates derived from for- 
mulas (JF> Max = P x E). This phenomenon 
was observed on T Phe and U Tuc by Campbell 
(1926) and on RT Eri (Payne and Campbell, 
1930); it is not correlated with the length of the 
next cycle, such as the shift of maximum de- 
scribed in the next paragraph. 

Abrupt shifts of epoch (fluctuations of 
phase) are also pointed out by Ludendorff 
(1928) on R Lup and by Shapley (1929) on Z 
Aqr. For the latter, the (O-C) curve 
demonstrates that the period remains constant 
over the interval covered by observations (1896 
to 1929), although there appears to be a jump 
in epoch of about -55 days near JD 2420900. 
This change seems to take place within about 
two periods, but the observations are not 


spaced closely enough to determine whether the 
change is due to an abnormal lengthening of 
a maximum or a minimum. New photometric 
and spectroscopic observations are needed on 
such stars. 

Multiple Periods. The light curves can also 
be analyzed by other means, such as the super- 
position of two or more distinct periodic light 
curves. Payne-Gaposchkin (1954b) notes that 
the ratio P /P, is approximately the same for 
the stars of a given spectral type, but it differs 
from one spectral type to another. She gives 
P 2 /P { ~ 9.4 for the M stars and 12.2 for N 
stars. The frequency distribution of P 1 /P x for 
M-type variables has a rather flat maximum be- 
tween 7.5 and 10, with a slight peak between 
9 and 10 (Houck, 1963). With a small sample 
of N-type stars, Houck shows a maximum fre- 
quency of P 2 /Pj between 12 and 13. These 
values agree well with Gaposchkin’s earlier re- 
sults. Houck (1963) gives six Mira variables 
with secondary periods: SV And, U Per, V 
Hya, V545 Cen, Y Per, and V1280 Sgr. Among 
the variables with M-type spectra, there are 
five with P /Pj ratio in the range 22.5 to 51.5 
and one N star with this ratio about 45.4. In 
these quoted papers, P 2 should be interpreted 
in terms of a beat frequency phenomenon 
resulting from interference between two slightly 
different periods within the atmosphere of the 
stars. 

Leung (1971) describes a method of analysis 
of superposition of two or more nonsinusoidal 
period components. Leung (1980) applies this 
method to 50-year observations on six selected 
M and S Miras: o Cet, R And, T Cep, T Cas, 
R Hya, and S CrB. Two components were de- 
termined for the majority of the observed Miras 
(Figure 1-14). For five variables, the period 
ratios, P 2 /P ]5 of the two individual nonsinu- 
soidal period components are around the value, 
1.04 ± 0.01. In general, the shorter period P { 
is associated with the larger amplitude of varia- 
tion, and the amplitude ratios, A 2 /A ]9 range 
from 0.15 to 0.28. The amplitude and period 
ratios of period components of S CrB do not 
fit into the general behavior pattern described 
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Figure 1-14. AA VSO observations of T Cas (10-day mean). The computed light curve (smooth curve) 
is based on two period components. The observations not included in the determination of period 
components are denoted as crosses. The abscissa is Julian Day (from Leung, 1980). 


here. The o Cet light curve needs the superposi- 
tion of three periodic components. 

Houck (1963) gives two periods for the C 
star, V Hya: P J = 533 days and P 2 = 6500 
days (beat frequency). These values are of the 
same order as those of Campbell (1943), who 
finds P { = 539 days and P 2 = 18 years. The 
latter period gives extra-deep minima reached 
in 1889, 1908, 1925, and 1943. Analyzing the 
data gathered by AAVSO from 1860 to 1968 
and using the method of analysis of Leung 
(1971), Leung (1973) points out one period of 
6670 days with an amplitude of about 3.5 mag- 
nitude. The other component has a period of 
530.7 days, an amplitude variation of 1.1 mag- 
nitude, and a shape close to sinusoidal, so that 
the larger period component has the larger am- 
plitude. Note that the opposite holds for the red 
supergiant variables (Leung, 1973). He con- 
cludes that V Hya may be a late Mira of lb 


type, for which the long-period component is 
not due to a normal mode of radial pulsation, 
and that the dip preceding the minimum is 
probably due to nonlinear effects or cold spots 
on the star surface. Reanalyzing the data of V 
Hya, Leung (1980) concludes that this carbon 
star has an extremely peculiar light variation 
and proposes four distinct period components. 

To conclude, we recall Leung (1980), who 
suggests that the primary period, P {9 is the 
pulsational period and that the modulating pe- 
riod, P , may be derived from a nonradial 
mode. The multiple periodicity may reveal in- 
formation on the modes of oscillation of these 
stars. 

Short-Term Fluctuations. In addition to the 
rapid evolution in the o Cet effective tempera- 
ture on 12 days reported by Joshi et al. (1980), 
we want to detail some scarce observations 
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which seem to indicate short-term fluctuations 
in the Mira variables. 

The J = 1—0, v = 1 , and v = 2 SiO maser 
emission lines of the Mira variable, R Leo, 
show phase-dependence effects in both line 
shape and intensity. At the premaximum 
phases, they are strong narrow features of de- 
creasing intensity. Around the visual maximum 
(<p = 0.98, 0.02, 0.04), they appear to be respec- 
tively about 30 and 15 percent broader without 
any trace of narrow features (Clark et al., 
1982). The effective duration of these changes 
continues up to <p « 0.40. 

With the Griffin-type radial-velocity spec- 
trometer and masks especially designed for 
Mira spectra, Pierce et al. (1979) find nightly 
variations of radial velocity of emission and ab- 
sorption lines in the spectra of Miras. More pre- 
cisely, the range in the measurements for each 
night is larger than the average uncertainty for 
ajset of observations covering 7 nights between 
^ = 0.9 and <p = 0.2. For example: on T Cep, 
the nights at <p = 0.02 and <p = 0.12 show a 
variable behavior of the absorption lines (the 
emission lines are too faint to be measured); on 
X Cyg, the nights at <p = 0.07, <p = 0.11, and 
ip =0.13 show a variable behavior of the emis- 
sion lines; and on o Cet, the radial velocity is 
variable in the absorption lines for the nights 
around ip = 0.76 and ip - 0.21 and in the emis- 
sion lines for the nights at ip = 0.12. The 
range of these nightly radial-velocity variations 
is between -2 and -13 km/s. These short-term 
variations should be interpreted by the interac- 
tion of shock waves with nonuniform atmo- 
spheric layers. Also, they should be linked to 
eventual flares such as those observed in R Aql 
(Woodworth and Hughes, 1973, 1977) or V Cyg 
(Querci et al., 1979); however, simultaneous 
observations are not available. Finally, it might 
be that these variations in radial velocity are a 
consequence of matter swept along magnetic 
tubes (see M. Querci, this volume). 

Correlations Between Different Observed 
Quantities. Numerous correlations between the 
observed quantities such as period length, light- 
curve shape, IR excess, visible and IR color di- 
agram, and OH radio-line separation exist in 


the literature. Apparently, some of them are 
surprising, because we do not understand for 
the moment why they happen. However, we 
shall note in detail here those of the greatest 
importance and shall endeavor to draw conclu- 
sions on the relationships between the derived 
physical parameters. 

The first correlations between the observed 
quantities were made at the Harvard Observa- 
tory. Campbell (1925) first demonstrated that 
a statistical relationship exists between the 
shape of the light curve and the period. The 
asymmetry becomes gradually more pro- 
nounced as the period increases. In all three 
spectral classes (M, S, or C), Mira light ampli- 
tudes are smaller for shorter periods and greater 
for longer periods. This was again explained by 
J.G. Garcia (1980). 

Using the data of Kukarkin et al. (1958), 
Merrill (1960) found that the period of M, S, 
and C stars is distributed in a broad maximum 
around 300 days for M stars, 360 days for S 
stars, and 400 days for C stars. Consequently, 
the C Miras have, on the average, the longest 
periods. 

Using two groups of Miras, one in which 
emission lines have low radial velocities (<20 
km/s) and another in which they have high ra- 
dial velocities (> 80 km/s; from Merrill's cata- 
log, 1941), Ahnert (1969) finds that the larger 
the scattering of cycles in percentage of periods, 
the earlier the subgroup of Me stars. A corre- 
sponding relation exists between the scattering 
and the period length, a relation which seems 
trivial because it is a well-known fact that the 
longer the period, the later the Me subgroup 
(Glasby, 1968). 

Harrington (1965) investigated 165 Miras 
from among the data of Campbell (1955) and 
found a correlation between the difference in 
magnitude of two successive maxima and the 
time interval separating the two maxima. When 
the interval is under the average value, the sec- 
ond maximum is brighter (45 percent of the 
considered stars show this correlation with a 
significance of at least 95 percent). These light 
curves may be interpreted by the outward prop- 
agating disturbances taking more or less time 
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to reach the star surface, this time being a func- 
tion of the energy carried. 

The kinematic studies of Feast (1962) imply 
that the shorter period Mira variables are Pop- 
ulation II stars, while those of longer periods 
belong to more intermediate populations. These 
conclusions suggest that the period amplitude 
of the Mira stars may be a function of their 
metal abundance. 

To determine the difference in metal abun- 
dances between short and long-period Mira var- 
iables, DeGioia-Eastwood et al. (1981) used the 
relationship derived between |W|, the velocity 
of the star perpendicular to the galactic plane, 
and d(U-B), the ultraviolet excess, given by Eg- 
gen et al. (1962). To calibrate the relationship 
between [Fe/H] and the kinematic properties, 
they use Carney’s (1979) relation between 
d(U-B) and [Fe/H]. They have examined two 
groups of Mira variables: the stars with periods 
149 < P < 200 days were found to be deficient 
in [Fe/H] by more than 1 dex, compared to the 
group with 350 < P < 410 days. The metal 
deficiencies in the short-period Mira variables 
are substantial, and consequently, the mass lost 
is proportional not only to the outflow velocity 
and the mass of the shell, but also to the abun- 
dance of heavy elements contained in the 
ejected matter. 

Eggen (1973) shows that the periods and col- 
ors of small-amplitude red variables are not re- 
lated. Eggen (1975) deduces a linear relation in 
the plane (log P,(R-I) o ) for the Miras, the in- 
dex (R-I) 0 being measured at phase 0.25 (ap- 
proximately medium light). This linear relation 
is given by: 

( R-I) 0 = -0.45 mag + 0.90 log P 

It is probable that this relation becomes non- 
linear for stars with periods longer than 500 
days that are redder than R-I =1.8 mag (Eg- 
gen, 1975); although the number of such stars 
is not important, observations are needed to de- 
termine the (( R-I) 0 , log P) relationship. The 
halo Mira variables have no correlation be- 
tween the period and the color R-I at phase 
0.25 (medium light). Feast (1980) points out 
that Miras show a slow rise in brightness with 


period, from M bo{ ~ -4.1 at 200 days to 
M boJ ~ -4.6 at 400 days. The short-period 
(~ 135 days) Miras do not follow this law and 
are quite fainter. 

The extension of the visible to the red and 
radio spectral ranges, and consequently, the de- 
tection of the IR excess and of the maser lines, 
leads to many correlations for an improved 
knowledge of the Mira behavior. Moreover, the 
analysis of M, S, and C stellar spectra using 
model atmosphere synthetic spectra is very use- 
ful because it shows that the stellar continuum 
could be reached in small intervals around 1.04, 
3.6, and 8 fim (see H.R. Johnson, this volume). 

Lockwood and Wing (1971) and Wing and 
Lockwood (1973) find a direct correlation be- 
tween the period length and the amplitude vari- 
ation at 1 .04 fim for a large sample of M stars, 
whereas Keenan et al. (1974) find a direct cor- 
relation between the mean spectral type at max- 
imum light and the period for the Me stars. 
Ukita (1982) shows (Figure 1-15) this relation- 
ship for the Mira stars with strong OH emis- 
sion, weak OH emission, non-OH Miras, and 
Miras which have not been observed in OH. 
Frogel (1971) established a correlation between 
the variation of the [1 .2 /im]-[3.5 fim] color in- 
dex and the amplitude variation at 2.2 fim for 
stars with little or large excess 10-jun radiation. 
A tendency for larger period stars to have 
longer amplitudes and larger [0.8 fim]-[2.2 *mi] 
or [I-K\ color indices was noted by Hyland et 
al. (1972). The continuum amplitude variations 
can be explained by variations in the effective 
temperature, whereas both the effective tem- 
perature and molecular opacities act together 
on other spectral ranges. Using the 2.2-fim spec- 
tral range, Harvey et al. (1974) again found the 
period-amplitude relation previously obtained 
in the visible and at 1 .04 fim; they also estab- 
lished a relation between the variation of the 
[3.5 fim]-[\0.0 fim] color index and the period 
amplitude. DeGioia-Eastwood et al. (1981) 
found that the [8.7 ^m]-[11.4 fim] color aug- 
ments monotonically with increasing period for 
a sample of 41 Mira variables of M spectral 
type. However, since the photometric phase at 
the time of each measurement had not been 
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Figure 1-15 . Plot of the 1. 04-yjn intensity variations against stellar period for oxygen-rich Mira variables: 
filled circles represent OH Miras in which OH emission was detected [OH luminosity =SD 2 >0.1 
Jy (kpcf]; filled triangles are OH Miras with weak OH emission [SD 2 < 0.1 Jy(kpc) 2 ]; open circles 
are non-OH Miras; dots are Miras not having been observed in OH or for which upper limits are 
poorer than the above criterion (from Ukita, 1982). 


taken into account (of course, how could it?), 
some of the scatter may be due to intrinsic 
variability. We conclude that the height of the 
silicate feature increases with the period ampli- 
tude, and consequently, the ratio of mass loss 
and dust formation, as well as the temperature 
of the silicate emitting layers, are directly linked 
to the stellar parameters which dictate the 
period amplitude. The height of the silicate 
feature does not seem to be due to a stochastic 
process in any Mira variable. 

Like DeGioia-Eastwood et al. (1981), we 
point out that the IR excess is not indicative of 
the total amount of mass in the shell, which is 
governed by the mass loss, but only of the mass 
of dust (and not gas) contained in the ther- 
mosphere , which extends to that radius in the 
shell where the dust is sufficiently warm to 
radiate significantly at 10 fim. Therefore, the 
[8.7 fxm]-[HA pm] versus period relation can 
be interpreted as an indication that the observed 
mass of dust in a given thermosphere is an 
almost linear function of the period. 

Bearing in mind that the Population II stars 


have a short period and that they are metal- 
deficient, we assume that two stars are equal 
in mass and luminosity and differ only in 
period; the previous correlation permits us to 
conclude that the Population II stars are also 
infrared-excess-deficient, and consequently, 
that dust formation may be a function of the 
metal abundance. 

In some stars, Wilson (1970) noted that the 
OH 1612-MHz masers and the IR variations are 
clearly correlated. Harvey et al. (1974) devel- 
oped a monitoring program to measure, at 
monthly intervals, the maser at 1612 MHz and 
the broadband infrared fluxes at six wave- 
lengths for OH/IR sources. They essentially 
confirm the observed correlation between the 
infrared and the 1612-MHz variability of the 
OH/IR stars and the OH clouds. This coupl- 
ing mechanism is consistent with a radiative 
pumping of the masers, possibly at 2.8, 35, 53, 
80, and 120 ^m (Litvak, 1969; Litvak and 
Dickinson, 1972; Elitzur, 1978; Elitzur et al., 
1976; Bujarrabal et al., 1980; Epchtein et al., 
1980). 
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A monitoring of OH maser emission from 
the M-type stars, IRC + 10011, NML Tau, U 
Ori, S CrB, WX Ser, U Her, R Aql, and R Cas 
is made by Jewell et al. (1979). The maser emis- 
sion is observed in the 1612-, 1665-, and 
1667-MHz lines of OH and 22-GHz line of 
H 2 0. They conclude that, in most cases (their 
Figures 7 through 13), for all main and satellite 
lines, the maser pumping mechanism must be 
radiative in nature. For the satellite line, this 
has already been established with some certain- 
ty by Harvey et al. (1974). The data presented 
by Fillit et al. (1977) also suggest a radiative 
pumping for the satellite line and, with less cer- 
tainty, for the main lines. For these latter lines, 
the work of Jewell et al. (1979) establishes con- 
clusively that the main lines are also related to 
the optical flux of the star (see later for con- 
troversy). Schwartz et al. (1974) suggest two 
possible correlations between the 1.35-cm H 2 0 
maser and the 2.2-/un infrared flux, either an 
exponential rate which fits the intensity of the 
maser versus the IR flux at 2.2 fim, or a linear 
relation between them after a threshold of the 
infrared flux value (which is probably corre- 
lated with some processes quenching the maser 
at low level of infrared flux or a possible colli- 
sional deexcitation of H 2 0 molecules). The 
pump may be linearly proportional to the in- 
frared flux above the threshold if the maser is 
saturated. 

Werner et al. (1980) simultaneously moni- 
tored far-IR, near-IR, and radio observations 
on five very red objects with 1612-MHz emis- 
sions. The direct comparison of these data 
strongly supports the hypothesis that the maser 
emission is pumped by 35-jun photons (see also 
Nguyen-Q-Rieu, this volume). 

In the case of U Ori, the comparison of dif- 
ferent line velocities produces three velocity 
groups (Cimerman, 1979). The most important 
is the group at -42 km/s, which includes the 
zero-volt optical absorption lines, one line of 
the SiO transitions, and an OH line of each fre- 
quency seen in all observations. The second 
group includes the optical emission lines and 
one line of the 1612-MHz features at -47 km/s. 
Lines scattered around the velocity of -36 km/s 


form the last group. Garrigue (1980) attempts 
to correlate the OH maser emission with the 
observed visible light curve. In his U Ori 
monitoring, he finds an entire period of light 
shifted by 2 mag above its usual mean value for 
each phase, whereas the light curve of this Mira 
is usually very regular. Garrigue and Men- 
nessier (1980) point out that the beginning of 
this epoch of higher brightness began approx- 
imately 650 days before the 1612-MHz flare 
discovered by Pataki and Kolena (1974b) in 
November 1974 and that the lifetimes of both 
visual brightness and 1612-MHz satellite in- 
creases have 550-day durations. This confirms 
that the OH emission in the satellite line is due 
to a perturbation which propagates throughout 
the stellar atmosphere. The proposed scenario 
is as follows: 

1 . An overbrightening perturbation occurs 

in the photosphere, disturbing the visual 
light curve during 550 days. 

2 . It propagates thoughout the stellar at- 

mosphere with a mean velocity of 20 
km/s (Slutz, 1976). 

3 . After about 650 days, the perturbation 

reaches the dust shell and disturbs the 
IR radiation emitted by it (at about 10 14 
cm; from Gehrz and Woolf, 1971). 

The shell IR radiation reaches the U Ori 
OH-emitting region (10 15 cm; from Reid et al., 
1977) in 1 day, and “by OH radiation coupling, 
is responsible to the maser effect during 550 
days” (Garrigue and Mennessier, 1980). This 
interpretation is a little different from the one 
proposed by Cimerman (1979), who does not 
take account of the 650-day delay between the 
beginning of the visible and the radio 
variations. 

Four papers (Dickinson et al., 1973; Harvey 
et al., 1974; Dickinson et al., 1975; Morris et 
al., 1979) argue in favor of a direct correlation 
between the expansion velocity given by the OH 
maser circumstellar lines and the period. Lepine 
et al. (1976) show that this correlation is doubt- 
ful for their Mira sample. Using the relation- 
ship of Figure 1-15 and the hypothetical cor- 
relation expansion velocity versus period 
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amplitude, one might conclude that a direct 
correlation exists between the OH expansion 
velocity and the amplitude variation on the 1104 
filter. Research in this way was attempted by 
Ukita (1982), who confirms the uncertainty 
stressed by Lepine et al. (1976). With data of 
type I and II OH Miras (Figure 1*16), he finds 
that the expansion velocity has a direct correla- 
tion with the amplitude variation at 1 .04 ^m 
and an inverse correlation with the period 
length. Knowing that the type II OH/IR Mira 
stars have larger color indices than the OH 
Miras and, consequently, have thicker 
envelopes (Harvey et al., 1974; Dickinson et al., 
1975; Olnon, 1977), Ukita (1982) assumes that 
the mass-loss rates increase in the order: non- 
OH Miras, OH Miras, and type II OH/IR 
stars. Ukita interprets his result in terms of a 
dust-driven wind enhanced by pulsation, a 
hybrid model of mass loss (Wood, 1979; 
Deguchi, 1980). 

Looking for the SiO masers in a large sam- 
ple of Miras, SRb, and supergiants, Spencer et 


al. (1981) conclude that no global properties of 
these stars correlate with the SiO luminosity; 
this indicates that the SiO masers are located 
very close to the stars, where the local condi- 
tions affect the maser intensity much more than 
the global properties of the photosphere, such 
as the stellar temperature, affect it. The local 
conditions could be governed by clumps of gas, 
convective cells, or turbulent eddies. The varia- 
tion of the SiO maser intensity and profile 
pointed out by Clark et al. (1982a, 1982b) in 
the Mira, R Cas, is a good example (see their 
Figures 1 through 5). Moreover, in Miras, the 
SiO maser flux is correlated with the bolometric 
flux (Cahn and Elitzur, 1979; Cahn, 1981) and 
with the infrared flux (Hjalmarson and 
Olofsson, 1979). (Let us remark that, when the 
lines are weak, broad, and roughly parabolic, 
they are due to nonmaser emission formed in 
a thick expanding envelope (Robinson and Van 
BlerKom, 1981).) 

Using VLBI observations, Bowers et al. 
(1980) correlate the mass-loss rates with the size 



Figure 1-16. Plot of expansion velocity versus the quantity A I x 400/P , for oxygen-rich Mira variables. 
(A I and marks are defined in Figure 1-15.) The crosses represent type II OH/IR stars (from Ukita , 1982). 
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of the masing regions; stars with mass-loss rates 
less than 10' 5 M© yr’ 1 have masing regions 
<1000 AU in extent, while stars with stronger 
mass-loss rates have larger masing regions. 

Wolff and Carlson (1982) point out a 
“general tendency of both line width and in- 
tegrated intensity of SiO thermal emission lines 
to increase with infrared color excess.” 

Scharlach and Woolf (1979) look for a cor- 
relation between the maser activity and the 
amount of ejected matter, which is correlated 
with the subphotospheric activity (Thomas, 
1973) as indicated by the hydrogen emission. 
They therefore searched for a correlation be- 
tween the maser activity and the H<5 emission 
line (which is the only Balmer line not disturbed 
by TiO absorptions); no correlation was found. 

Modes of Pulsation of Long-Period Variables. 

Miras are radial pulsators and are driven mainly 
by hydrogen ionization effects (e.g., Cox, 
1984). The real germane problem is to know 
which radial mode is the primary mode of 
pulsation, the fundamental or the first over- 
tone. Both modes have been proposed by dif- 
ferent groups through comparisons between 
observational values of the pulsation constant 
Q and theoretical Q-values. 

The pulsation constant Q (in days) relates 
the period P (in days) and the density q by: Q 
= P(q/q © ) l/2 , otherwise: 

Q = P(M/Mo) m {R/RoY v2 , 

with M , the mass, and R , the radius (i.e., a 
period-mass-radius (PMR) relation exists). 

A primordial observational constraint to the 
Q determination concerns the knowledge of 
radii. Their evaluation by direct methods (giv- 
ing angular diameters) such as occultation 
technics or speckle interferometry, or by 
photometry is debated by Wood (1981) and 
Willson (1982). (See also the section on Radii.) 
In the direct method, it is shown that radii vary 
from one spectral region to another; also 
distances to individual objects must be known. 
The photometric approach involves the 
luminosity definition: 

l = 4tt r 2 o rL 


The radius may be expressed in the function of 
L, hence the absolute bolometric magnitude 
Af bol , and of F ff . Therefore, knowing obser- 
vationally the A/ bol /log P relation and the log 
F ff /P relation for LPV’s, we get observa- 
tional Q M~ U1 /P relations. However, the 
uncertainty in the effective temperature deter- 
mination (hence, in the radius determination) 
for the Miras remains the major problem in 
estimating Q-values from observations. (See the 
extensive discussions in Willson, 1982, and in 
Wood, 1981, 1982.) This is shown in the sum- 
mary of observed Q-values given by Fox and 
Wood (1982, and references herein): (1) for the 
Population II small-amplitude variables (e.g., 
in globular clusters 47 Tuc, M4, M22, and co 
Cen), using the Ridgway et al. (1980) 
temperature scale (an extension to Miras of a 
scale based on nonvariable M giants), the pulsa- 
tion mode is the first overtone for stars in the 
two first-quoted clusters, and the fundamen- 
tal for stars in the two latter, without clear ex- 
planation; (2) for Miras in the solar 
neighborhood, using a blackbody temperature 
scale based on lunar occultation angular 
diameters, the mode is the first overtone, 
whereas the use of the Ridgway et al. F ff scale 
might lead to the fundamental mode; (3) in the 
Magellanic Clouds, for the two groups of long- 
period variables (see the section Aspects of 
Evolution of Long-Period Variables in the 
Magellanic Clouds ), the pulsation mode ap- 
pears to be first overtone if blackbody temper- 
atures are used, or fundamental mode if the 
Ridgway et al. F cff are favored. In fact, one 
can have confidence in the blackbody color 
temperatures for LPV’s of large amplitude that 
are linked to infrared lunar occultation 
diameters of Mira variables (as described in 
Glass and Feast, 1982); therefore, first-overtone 
pulsations are favored. 

Pulsation models, as well as linear 
nonadiabatic models (Fox and Wood, 1982), 
have been calculated as nonlinear nonadiabatic 
models (Wood, 1974; Tuchman et al., 1979; 
Ostlie et al., 1982). In the extensive grid of 
linear nonadiabatic pulsation models of Fox 
and Wood (1982), theoretical values of the 
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pulsation constant are derived for the first three 
radial pulsation modes in red giants (galactic 
disc and Population II) and supergiants. Main 
results are as follows: 

• In the standard models, for Miras in the 
galactic disc, Q is not constant for either 
mode. Q-values for the first overtone 
asymptotically approaches Q x ~ 0.04 at 
low luminosity or high mass, and short 
P, but Q { may take values up to about 
0.075 (see Figure 3 in Fox and Wood, 
1982); no simple period-mass-radius 
relation exists. For the fundamental 
mode, Q-values increase with P 9 and the 
relation P M* R ~ b = Q with a - 1.0 
and b - 2.0 holds. 

• The ratio/yPj can take values from 
-2 to 3 for values of P { < 630 days; 
P Q /P { ~ 7 occurs when P Q ~ 2000- 
5000 days in some massive luminous 
stars. Beat periods exist between and 
P 2 of ~3 P { to -20 P x (Wood, 1982); 
some observed beat periods are given in 
the section Multiple Periods. 

Note that, for the LPV’s in Magellanic 
Clouds (models with slightly decreased metal 
abundances), the Q-values for the first overtone 
tend to cluster around Q (days) = 0.38. 

Therefore, although the ratio of PfP { 
seems to be well determined theoretically, deter- 
mination of the mode of pulsation itself of the 
Miras requires more than the Q-value deter- 
mination. Wood (1981) suggests that the 
strongest evidence in favor of the first-overtone 
pulsation mode is the presence of secondary 
periodicities caused by resonant coupling be- 
tween the first-overtone (corresponding to the 
primary mode of pulsation) and the fundamen- 
tal mode (the secondary mode of pulsation) 
when the latter has a period almost twice that 
of the first overtone. In fact, Wood (1982), 
among various examples, quotes massive 
supergiants (see Figure 1-3) which present 
secondary periods with P — 2000-7000 days 
(i.e., about seven times the length of the 
primary pulsation period) in agreement with the 


P /P - 7 found in supergiant models of Fox 
and Wood (1982). Also, carbon stars among 
AGB variables are known to have such an 
observed P Q /P l ratio. Another support to the 
first-overtone primary mode of pulsation comes 
from Figure 1-3. As commented in the section 
Aspects of Evolution of Long-Period Variables 
in Magellanic Clouds , an AGB LPV switches 
from the first-overtone mode to the fundamen- 
tal mode when crossing the dotted line cor- 
responding to the long-period edge of the AGB 
region, to subsequently form a planetary 
nebula. If the primary mode is the fundamen- 
tal one, the event that could suddenly stop the 
AGB evolution of an LPV still having substan- 
tial envelope mass is not clear. Evidence for the 
fundamental mode of pulsation in Miras has 
been given by Willson (1981, 1982, and 
references herein). Again, the controversy 
comes mainly from the used radii. Also, 
Willson analyzes shock propagation in Mira at- 
mospheres to determine Q-values indicating the 
fundamental mode of pulsation. Wood (1981, 
1982) extensively discussed her arguments and 
did not find them convincing. 

The SRa-SRb-SRc Variables 

Sir Williams Hershel discovered the first two 
variable stars of this type: fi Cep in 1782 and 
a Her in 1796, with their magnitude varying in 
a “random” manner between 4-6 and 3-4, 
respectively. John Hershel attributed these ir- 
regularities to some periodic veiling effects 
caused by interposing dark nebulous matter. 

The amplitudes of SR variables are smaller 
than those of Miras, and the lengths of their 
cycles are generally smaller than the period of 
the Miras, except for certain supergiants in 
which periods of many years are observed (e.g., 
a Ori). 

It is very difficult to define a period for stars 
in which the random brightness fluctations are 
more marked than the regular period. The 
duration of the observed cycle often varies by 
more than 1 month, and the individual max- 
ima by more than a magnitude. However, the 
classification into the three following groups 
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(more or less arbitrary) is generally accepted: 

1 . The semiregulars of late spectral class 
(M, S, and C), denoted SRa, are giants. 
Many of these stars differ from the 
Mira types only in the smaller amplitude 
of light variations. Frequently, their 
curves have strong variations from one 
cycle to another. Typical representative 
stars are RU And, S Aql, Z Aqr, T Cnc, 
WZ Cas, T Cen, RS Cyg, RS Gem, and 
R UMi. 

2 . The SRb-type semiregulars have a poor- 
ly expressed periodicity. Different 
periods of the individual cycles prevent 
the prediction of the epochs of 
minimum and maximum brightness. 
These stars sometimes temporarily 
replace periodic changes by slow ir- 
regular variations or even by the con- 
stancy of the brightness. Typical stars 
are V Aqr, V Boo, RX Boo, UU Aur, 
X Cnc, Y CVn, TT Cyg, R Dor, RY 
Dra, UX Dra, T Ind, R Lyr, VY Leo, 
W Ori, L 2 Pup, q Per, SW Vir, and S 
Set. The SRb stars are giants. 

3 . The SRc-type semiregulars are super- 
giants with an SRb behavior. The visual 
light changes are generally of the order 
of a single magnitude or less (Maeder, 
1980). Typical representatives of the SRc 
type are VY CMa, Cep, RS Cnc, RW 
Cyg, TV Gem, Y Lyn, a Ori, S Per, and 
VX Sgr. 

It appears that we find a larger fraction of 
semiregular variables among the carbon stars 
(C spectral class) than among the M and S stars. 
Almost half of the known semiregulars are N 
stars, the remainder being of M class; very few 
of S and R classes are known. Moreover, it is 
of interest to note that the N-type stars have 
larger primary and secondary cycles than the 
M class stars (for example, see Table 9 in 
Glasby, 1968). 

The frequency-distribution curve of class C 
Mira and SR stars are very different. The 
former has a large maximum between 350 and 


450 days, while the latter has two maxima 
around 150 and 400 days, clearly separated by 
a minimum around 270 days (Alksne and 
Ikaunieks, 1971). Analyzing this minimum in 
the distribution, Guzeva Yakimova (1960) 
noted the absence of carbon stars with periods 
(Miras) or cycles (SRa) between 250 and 285 
days. With observations based on 316 stars 
(mainly spectral class M), Glasby (1968) uses 
the shortest cycle for SR stars (cycle intrinsically 
related to a possible pulsation of the star itself) 
and constructs a frequency-distribution curve 
for the cycle length. Two maxima are detected, 
around 85 and 135 days, on a period range 
from about 30 to 1000 days or more. 

Different Time-Scale Light Variations of 
Semiregular Variables. The first approach to 
these stars could be light curves obtained by a 
visual estimate of the brightness; such tasks are 
well suited to amateur observers. Visual curves 
of semiregular stars do not have the strong 
similarity among them as do the Miras, but are 
quite individual. Moreover, instead of the mean 
periods of the Miras, the SR variables may be 
characterized by a form of periodicity which is 
hidden by irregular variations in brightness 
(Glasby, 1968). 

Among the semiregular variable stars, the 
SRa are most similar to the Miras, with which 
we find many common characteristics such as 
Balmer emission around some light maxima as 
in WZ Cas, changing periods (Hoffleit, 1979), 
and other characteristics. 

As for the SRb stars, there is often an er- 
roneous classification between them and the ir- 
regular variables, Lb, because the observations 
are not continued for a sufficient duration. 
Maran et al. (1980) give a good example with 
VY Leo (56 Leo), which clearly shows the slow 
irregular variations that temporarily replace the 
periodic changes. Two groups of Brazilian 
astronomers have separately detected variations 
with a time scale of about 1 hour on the M8 II 
star, R Crt: Gomes Balboa et al. (1982) in the 
22-GHz H 2 0 line and Livi and Bergmann 
(1982) in the DDO magnitudes and colors (over 
TiO bands and Ca I lines). R Crt is the first SRb 
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Figure 1-17. Photometric light curves of the SRb variable, TW Hor (from Bouchet, 1984a). 


variable to show such rapid variations. 

Bouchet (1984a) analyzes the C star, TW 
Hor, by broadband photometry from the U 
band to 30 fit n (i.e., using U, B, V, R, I, J, K, 
L, M, N, Q, and P filters; Figure 1-17) carried 
out over a 4-year period and by spec- 
trophotometry from 1 to 13 jtm. The variations 
of brightness are similar in the J, H, L, and V 
bands. They are less pronounced in the K and 
M bands. The V-R and K-/ indices vary slightly 
in phase with the V-magnitude. The B- V index 
remains nearly constant over the period, but the 
U-B varies oppositely to the other colors and 
magnitudes. Strong and rapid oscillations on 
U-B are pointed out between = 0.80 and <p 
= 0.90 (see his Figure 2); they are probably 
connected with the X3280 Fe II VI emission-line 
variations (Bouchet et al., 1983; see also M. 
Querci, this volume). These oscillations should 
duplicate rapid temperature or opacity varia- 
tions due to nonthermal chromospheric excita- 


tion (Querci and Querci, 1983, 1985a). Surely, 
important scientific information is contained in 
these variations, but many opacity contributors 
such as C 3 , violet CN, SiC, or graphite are in- 
cluded in the flux of each filter, and the restora- 
tion of each one is impossible. A new 
photometric system for carbon stars is needed 
(Querci and Querci, 1985b). A feature noted in 
the Bouchet’s (1984a) paper is the variation on 
the N flux where SiC particles may form, which 
is either due to the variation of the amount of 
emitting SiC or to the fading of the exciting 
emission. Moreover, the time-scale variations 
over a week on N and P filters are not quoted, 
but they are not negligible. At phase 0.25, the 
N flux varies from 8.1 to 6.3 x 10' 25 Wm" 2 
Hz' 1 , while Q varies from 2.8 to 5.7 x 10" 25 
Wm' 2 Hz' 1 . 

Lastly, we describe the supergiants (SRc). 
Because the time interval between two con- 
secutive maxima for the supergiants is between 
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several hundred and several thousand days with 
amplitude variations which are generally small 
(<1.5 mag), several years of observations are 
needed to decide whether a supergiant such as 
q Cas belongs to the SRc class or to the Lc class 
(Smith, 1976). Maeder (1980) observed on the 
Geneva photometry file that the cyclical 
amplitude variations of the SRc stars decrease 
by several tenths of magnitude over some years; 
he concludes that “over half a century or more 
they may exhibit changes up to 2 magnitudes.” 
Moreover, he points out that the limit in the 
HR diagram where the light variation 
amplitudes of the G-M supergiants become im- 
portant coincides very closely with the limit of 
appearance of very deep convective zones in 
these stars. Using these two observations, we 
may conclude that the SRc supergiants stay no 
more than 100 years near the Hayashi limit 
because, in this area of the HR diagram, the 
amplitude variations are at their maximum and 
the observations show that the light maxima 
decrease about 2 magnitudes in one century. 

The hypergiant, q Cas, shows large varia- 
tions in luminosity and spectral class. Before 
1930, it had a K class spectrum; in 1943, 
Morgan et al. (1943) classified it as F8Ia,b, and 
during 1946-1947, the star underwent a deep 
light minimum (A m v ~ 2.0), probably because 
of the large amount of ejected matter, and con- 
sequently, the spectrum looked like that of an 
M star. In 1948, Greenstein (1948) pointed out 
signatures of matter falling back on the star a 
day after a sudden expansion. Joshi and Rau- 
tela (1978) found that the color temperature of 
q Cas was higher than that of 6 CMa (a non- 
variable star of the same spectral class) by about 
650 K in 1920 and 100 K in 1974. As quoted 
by Glasby (1968): “One marked peculiarity of 
the light curve of q Cas is the apparent existence 
of three fairly distinct modes with quite abrupt 
changes occurring from one to the next. Such 
changes have been observed in 191 1 , 1922, and 
1946 when a very deep minimum occurred 
which was well followed visually and photoelec- 
trically.” 

Some other observational proofs of long- 
term variations have been made by the analysis 


of the variable polarization. For example, 
Tinbergen et al. (1981) interpret the yearly 
variations of a Ori and a Sco by very large 
eddies (size about 1.5 x 10 8 km) with a 
lifetime of 1 year. One, or at most two or three, 
of these large moving elements are seen at the 
same time; this agrees with Schwarzschild’s 
(1975) estimate that a cell moving up and down 
in the stellar atmosphere with a sonic speed 
(supposedly constant and equal to 5 km/s) 
crosses the whole stellar diameter in 1 year. In 
a Ori, Goldberg (1979) finds random fluctua- 
tions on the time scales of 1 year or less noted 
above. He also finds a cycle of about 5 years 
with both brightness and photospheric radial- 
velocity variations, which demonstrates the cor- 
rectness of the period of 5.781 years of Jones 
(1928), Sanford (1933), and Pettit (1945). 
(However, see the discussion by M. Querci, this 
volume.) 

In addition to these long-term disturbances, 
night-to-night variations are observed in some 
super or hypergiants. In q Cas (Figure 1-18), 
Joshi and Rautela (1978) show variations in the 
slope of the Paschen continua and the Balmer 
jump and suggest that they are brought about 
by the circumstellar cloud surrounding the star. 
The daily disturbances are probably due to local 
interaction between the photosphere and sur- 
rounding matter in small puffs or to local mo- 
tions in the atmosphere itself, leading to ir- 
regular variations of brightness and polari- 
zation. Sargent (1961) suggested that the weak 
Balmer absorption lines of q Cas during the late 
1950’s are probably filled in by emission lines, 
and he argued that there is no deficiency of 
hydrogen in the irregular variables. We con- 
clude that a monitoring with high-resolution 
spectroscopy could clarify the problem. 

Looking for the middle-term color changes 
shown by q Cas (Beardsley, 1961), Landolt 
(1968) concludes that “such a very long-term 
project as this color change problem is one to 
which serious amateur-astronomers using stan- 
dard UBV filters, a photoelectric cell, and a 
small reflector could usefully contribute.” 

The above variabilities remind us of those 
observed in the hotter supergiants, from A to 
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Figure 1-18 . Energy distribution curves of 
q Cas and S CMa: filled circles = q Cas; open 
circles = S CMa (from Joshi and Rautela, 
1978). 

G spectral classes, with the three different time 
scales: (1) the nightly variations, (2) the cyclic 
variations of several hundred days, and (3) the 
long-term variations of several years which are 
sometimes quasi-catastrophic for the star. 

Correlations Between Observed Quantities. 

Definite correlations between periods, absolute 
magnitudes, spectral classes, and phase lags are 
less obvious in the SR variables than in Miras. 
Generally, the clouds of points in the figures 
(Feast et al., 1980) are scattered around the 
mathematical relation deduced from these 
figures. However, efforts have been made to 
rectify the situation: Glasby (1968) reports 
studies from Palmer which conclude that, sta- 
tistically, the amplitude variations may be ac- 
cepted as constant. Or rather, he points out that 
the behavior of the maximum and minimum 


apparent brightness is virtually the same in the 
SRa stars, whereas it is very different in the 
Miras, in which two successive maxima could 
vary appreciably while the minima show very 
little variation. Moreover, he notes that the 
longer the period, the redder the spectrum — 
an association also noted by Strohmeier (1972). 

Although the monitoring of H 2 0 micro- 
wave (Schwartz et al., 1974) on the SR vari- 
ables, RX Boo (SRa), NML Cyg, and VX Sgr 
(SRc), was of such low signal-to-noise ratio that 
no correlation between microwave variations 
and visual luminosity was possible, Cox and 
Parker (1979) found a steady increase of the 
H 2 0 line with no sign of cyclic variation seen 
at optical wavelength in RX Boo (SRb). More- 
over, sharp maxima appear in the maser 
strength with no apparent correlation with the 
visible. The monitoring of the two H 2 0 peaks 
of the SRb star, R Crt, gives a different 
behavior of these two features against the visual 
brightness. Spencer et al. (1981) conclude that 
the SiO maser intensities are correlated not with 
the stellar properties but with local conditions 
in the stellar atmosphere of the SRa, SRb, and 
SRc stars (like the Miras). When the visual 
brightness of VX Sgr (SRc) monotonically in- 
creases by about 3 magnitudes (Schwartz et al., 
1974), the 2.2-fim flux increases by a factor 2. 
Then, and then only, does the H 2 0 microwave 
emission appear. It seems that the masering 
regions and the photosphere of the star are not 
usually coupled mechanically because the stellar 
photosphere suffers some periodic variations, 
but the masers have quite nonvariable radial 
velocities. 

In his investigation on the SRa carbon star, 
Y CVn, with photoelectric B, V, and R filters 
and spectroscopic observations, Vetesnik 
(1982b, 1982c) gives a C 2 Swan radial-velocity 
oscillation of ±2 km/s. Although a gap exists 
in the photometric observations, he suggests a 
likely relationship between the light variations 
and the radial-velocity oscillations. Moreover, 
he gives the light elements of the minima by 
JD.. = 2436097.3 + 251.6E. However, this 

Mm 

new period disagrees with the mean period 
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given in the GCVS (P = 158 days); the chang- 
ing period is quite evident unless there is a 
mistake in the GCVS. 

Another SRa carbon star, UX Dra, shows 
a light curve with a marked regularity, alter- 
nating one deep and one shallow minimum, 
which suggests the light curve of a photometric 
binary (Vetesnik, 1982a). On the other hand, 
the radial velocity deduced from the C 2 Swan 
has a half-amplitude of 2.2 km/s and a period 
of 340 days, which is double the mean period 
obtained with the light curves. Confirmation 
by continuous observations is needed. 

Referring to 20-year observations in the 
Geneva photometry, Maeder (1980) gives the 
V amplitude variations for several SRc 
supergiants of type la, lab, and II. For a given 
spectral type, the larger the brightness ampli- 
tude, the more luminous is the supergiant; how- 
ever, Strohmeier (1972) argues that the cycle 
length has no apparent correlation with the 
luminosity. This apparent contradiction in- 
dicates that a very high quality of measurement 
is necessary for investigating SR stars. (The 
Geneva photometry is a good example.) 
Moreover, Maeder (1980) points out a new 
observed correlation for the G-M supergiants; 
for the same amplitude variation, the higher the 
luminosity, the earlier the spectral type is. 
Observations of numerous supergiants are 
necessary to confirm Maeder’s conclusions 
because he has used only three la supergiants, 
nine lab and three II supergiants of M0-M2 
class, and three lab and two lb supergiants for 
spectral classes later than M2. Because the ir- 
regular variable stars are generally poorly 
monitored by ground-based observers, 
numerous observations are needed to con- 
tribute to the understanding of these stars. 

Nevertheless, from 24 red supergiant vari- 
ables in the LMC, Feast et al. (1980) deduce 
a period-luminosity relation. The measurements 
of Glass (1979) come close to a theoretical P-L 
relationship derived from the work of Stothers 
and Leung (1971). 

The satellites hail a new epoch. Mar an et al. 
(1980) demonstrate that VY Leo is clearly an 
SRa variable with a period of about 1 year, 


although it was classified Lb in the GCVS. 
Thus, this star can no longer be used as a nor- 
mal giant (Wing, 1980). Consequently, it is 
urgent to use satellites for long series of uninter- 
rupted observations of selected variable stars; 
such sustained observations will help to point 
out the various time-scale phenomena involved 
in these star atmospheres. 

Superposition of Periods, Chaos, or Random- 
ness in Light Curves. A primal question is: are 
the semiregular stars multiperiodic, chaotic, or 
truly random? 

In the last decades, astronomers have mainly 
developed multiperiodic analysis, and many ap- 
parent light-curve irregularities may be explain- 
ed by a complex mixing of two or more oscilla- 
tions, each of which is more or less regular and 
varies independently. Harmonic analysis of the 
observations must be employed to identify the 
individual oscillations. RS Gem (SRb, M II) is 
a representative star, with a secondary varia- 
tion which is observed as a second maximum 
progressing along the light curve relative to the 
primary cycle (Figure 1-19). UZ Per is one of 
the SRa variables which has been studied in 
depth over several years; it apparently has three 
cycles of different length: a short period of 90.8 
days, a longer period of 922 days, and a third 
one which is estimated to be more than 5000 
days long (Glasby, 1968). 

Variable stars with a double period are more 
often observed among C stars than among M 
and S stars. A list of those stars are given in 
Table 26 of Alksne and Ikaunieks (1971). 
Values of the ratio of the secondary to the 
primary periods (P 2 /P,) are quoted in the 
works of Payne-Gaposchkin (1954b), Houck 
(1963), and Alksne and Ikaunieks (1971), in 
which the majority of the stars are SRa or SRb 
of both M and C classes. We will not discuss 
it again. Leung (1980) found that this period 
ratio (secondary to primary) is about 10 for the 
M supergiants (SRc). Stothers and Leung (1971) 
suggest two possible explanations for these two 
periods. First, the primary period is connected 
to the radial fundamental mode of pulsation, 
and the long secondary period is tentatively in- 
terpreted as the convective turnover time of 
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giant cells in the stellar envelope. A second 
possible interpretation is that the primary 
period may represent the first radial overtone, 
while the secondary period may represent the 
fundamental mode itself. Stothers (1972) 
prefers the first explanation, while Fox and 
Wood (1982) argue in favor of the second one 
and give a period ratio of about 7 (see the sec- 
tion Modes of Pulsation of Long-Period 
Variables and de la Reza, this volume). Figure 
1-20 shows the S Per light curve and the com- 
puted two-component light curve from Leung 
(1980). 

Nowadays, two other types of physical ex- 
planations for the behavior of the semiregular 
and irregular stars are proposed: the ran- 
domness and the chaos. In the randomness pro- 
cess, the events are largely independent of each 
other, and it is impossible to make accurate 
predictions of any range of variations with cer- 
tainty. However, in a chaotic process, which 
has been mainly analyzed in mathematics and 
physics, the observed variations reflect the 
unstable motions which are produced by a col- 
lective and cooperative behavior of the matter 
when some acting forces are largely amplified 
(Whitney, 1984). Nonlinear equations are used 
to investigate such behavior. 

Here, we make two observational remarks: 

1 . The analysis of semiregular and irregular 
variables requires a much longer series 
of uninterrupted observations than have 
been done until now; only very few stars 
(Houck, 1963; Glasby, 1968; Hoffleit, 
1980; Leung, 1980) have been monitored 
to this extent. 

2 . Although observations of new periods in 
cepheids or binaries with improved ac- 
curacy permit us to obtain more precise 
quantitative parameters of these stars, 
the accuracy of the observations on the 
semiregular stars may determine the 
nature of the variations amidst the pro- 
cesses described above and may allow us 
to choose between two competing phys- 
ical theories (Whitney, 1984). Observa- 
tions with improved accuracy are crucial 
for semiregular and irregular stars. 


As we know, there are many more irregular 
variables among the carbon stars than among 
the M stars. Let us give some examples of 
semiregular C stars with apparent behavior 
changes in their light curve, illustrating chaotic 
variations. First, we quote the SRa star, S Cam. 
This star was observed from the beginning of 
the 20th century. Ludendorff (1923) derived a 
period of 325.5 days and a light amplitude of 
2.8 mag from his 1911-1912 observations. 
Later, Campbell (1941, 1947) and Payne- 
Gaposchkin (1944) indicate a period of 326.5 
days and a light amplitude of 2.5 mag. Nielsen 
(1952) obtained a period of 324 days. All these 
results seem homogeneous. However, May all 
(1960a, 1964, 1966) suggests that S Cam light 
amplitudes decreased to 2.0 mag, and Romano 
(1950) points out time intervals with constant 
brightness and other ones with additional waves 
on the light curve. From her 1966-1971 obser- 
vations, Krempec (1973) deduces a period of 
326.4 days and a variable behavior in the light 
amplitude. During some time intervals, the light 
amplitude is above 1.3 mag, and during some 
other ones (presenting very flat minima), the 
light amplitude is only about 0.8 mag. During 
the Krempec’ s observational period (1966- 
1971), the SRb star, RV Aur, presented regular 
time variations with a period of 229 days, but 
the behavior of this star was sometimes quite 
irregular. UV Cam is another SR star with an 
apparent period (294 days) and considerable ir- 
regularities (Krempec, 1973). U Cam is one of 
the SRb stars with multiple periods. Payne- 
Gaposchkin (1944) finds three periods for this 
star (223, 435, and 3000 days). Kukarkin (1949) 
confirms the two periods of Jacchia (1933): 419 
days (as the period of fundamental oscillations) 
and a superposed wave of 3000 days. Krempec 
(1973) shows that the period of 400 days is 
probable. The light variations were often 
irregular. 

As quoted by Krempec (1973), there is no 
regularity in the light variations of the SRb star, 
X Cnc, and the observations give considerably 
different periods, such as 970 or 365 days. The 
light amplitude does not exceed a few tenths of 
magnitude, but Mendoza (1967) and Eggen 
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Figure 1-21. Light curve of V Aql (from Glasby, 1968). 


(1967) stated that the photometric color index, 
B-V , ranges from 2.9 to 3.4 and 3.0, respec- 
tively. Hagen (1925) notes that X Cnc is im- 
bedded in a dark symmetric cloud of size 0°8. 
This cloud is confirmed by a strong continuous 
ultraviolet absorption (Mendoza, 1967) and a 
strong continuous infrared emission (Gillett et 
al., 1971). TT Cyg, also an SRb carbon star, 
is an astonishing object; it shows very rapid ir- 
regular light variations during the 1923-1924 
and 1928-1929 observations of Parenago 
(1938). During other time intervals, it has a 
period of 118 ± 18 days (Payne-Gaposchkin, 
1944; Krempec, 1973). 

To conclude this analysis of the semiregular 
stars, we emphasize some occasional striking 
behavior like the one observed by Garcia et al. 
(1977) on the SRb star, R Dor. On the basis of 
455 visual observations, these authors show the 
occasional appearance of an 8-day harmonic 
period besides the 332-day fundamental period. 
At the present time, it is difficult to give a 
physical explanation to the excitation of the 
8-day mode. 

The Lb and Lc Stars 

Generalities. The Lb and Lc stars, giants and 
supergiants, respectively, are slowly varying 
without any trace of periodicity and with small 


amplitude (<1.5 mag). Deep minima appear 
at irregular intervals. The light curves seem to 
be composed by slow waves of varying 
amplitude. Strohmeier (1972) notes that an ex- 
tremely weak periodicity occasionally appears. 

Some variables are initially classified Lb or 
Lc and are later shifted to other variability 
classes as new observations reveal unsuspected 
periodicity. Therefore, the red irregular stars 
with unknown spectral class and luminosity are 
first assigned to the irregular star class while 
awaiting new observations. 

Few such stars are regularly observed by 
spectroscopy. Even today the presence of 
Balmer emission lines has not been successfully 
demonstrated in these variable stars. 

Although the Miras of spectral class M are 
scattered over a large spectral range— M2 to 
M 10— the irregular stars, as well as the semi- 
regulars, are concentrated at M6 and M6.5, 
with none later than M7 (Cameron and Nassau, 
1956); the latest should be RX Boo (M7-M8). 
Joy (1942) showed that the concentration of ir- 
regular M stars occurs at M5 and M6. This dis- 
crepancy occurs because a small difference ap- 
pears between the Mt. Wilson (Joy) and Case 
(Cameron and Nassau) systems of classifica- 
tion. We must keep in mind that these dif- 
ferences are partially due to the variation of 
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Figure 1-22 . (a) Portions of the visual light curves of two supergiant variables of small amplitude . 
The ordinates are in magnitudes, and the abscissas are in thousands of days from an arbitrary zero- 
point. Note that the scales for the two variables are not the same . The curve for fi Cep is from Hassen- 
stein (1938), and that of VV Cep is from Fredrick (1960). (b) Portions of the light curves of three 
red supergiant variables of large amplitude. The ordinates of S Persei and VX Sagittarii are in units 
of visual magnitude, while photographic magnitudes are given for VY Canis Majoris. The abscissas 
are in thousands of days from an arbitrary zero-point . The visual light curves are based on AA VSO 
observations, and the photographic curve was determined by Robinson (1970) from Harvard plates 
(from Smith, 1976). 


spectral class of the star during its light varia- 
tions. Wilson and Merrill (1942) give a long list 
of irregular variables with M spectral types. 
They add two S spectral-type stars: SU Mon 
and AD Cyg. The stars n Cep, a Ori, a Her, 
and a Sco are also classified irregular by 
Strohmeier (1972). Many M class irregular 
variables are connected with young O-associ- 
ations of Population I (Crawford et al., 1955). 

The carbon Mira variables are found in the 
low-temperature C6-C7 subclasses, while the 
semiregular and irregular variables are chiefly 
found in the earlier subclasses (Alksne and 
Ikaunieks, 1971). Typical representative Lb car- 
bon stars are V Aql (Figure 1-21), W CMa, 
V460 Cyg, U Hya, V Hya, T Lyr, BL Ori, TX 
Psc, and VY UMa. 

In any case, it is generally difficult to decide 
whether a supergiant belongs to class SRc or 
class Lc because the time interval between suc- 
cessive maxima is very long. This time varies 
from a few hundreds to a few thousands of 
days, and extensive observations are needed 


before a conclusion can be reached. This raises 
the question that another kind of distinction 
between the supergiants must be investigated. 
Smith (1976) suggests a much more physical 
splitting among the supergiants: those with 
small amplitude variations and those with large 
amplitude variations (Figure l-22a, b). In the 
first group, we find \i Cep and VV Cep, the 
latter being an eclipsing binary. Their mean 
maximum amplitude is 1 . 1 mag with a disper- 
sion of about 0.2 to 0.3 mag. The second group 
includes stars like S Per, VY CMa, and VX Sgr. 
Using the published line curves of Robinson 
(1970), Dinerstein (1973), and Smith (1974), 
Smith (1976) found a mean maximum 
amplitude of 4.2 visual magnitudes with a 
dispersion of 1 . 1 mag among these three stars. 
The light curves of \x Cep and VV Cep suggest 
that they have something in common, but they 
are different from the other three supergiants 
which are rather inhomogeneously classified in 
GCVS: S Per as SRc, VY CMa as Lc, and VX 
Sgr as SRb. 
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It seems to be necessary to support this divi- 
sion into small and large amplitude variables 
by studying other stellar properties and obser- 
vations of other stars. Humphreys (1974) found 
that the three stars, S Per, VY CMa, and VX 
Sgr, share a few peculiarities, such as large ex- 
cess radiation between 0.7 and 1.5 /xm and be- 
tween 1.5 and 8 ^m. This last excess has a slope 
like that of free-free emission in VX Sgr and 
S Per. Moreover, these stars show absorption- 
line weakening in their near-infrared spectra, 
silicate emission features between 10 and 20 fix n, 
and strong OH and H 2 0 emissions. S Aur and 
Y Lyn can be added to this group of large 
amplitude variables. 

Among the observed variable supergiants, 
the number of large amplitude variables is very 
limited. Then, if we suppose that the 
supergiants cross through the small and large 
amplitude phase, we could conclude that the 
latter phase is very short. Moreover, no 
theoretical works give any information on the 
order and the duration of these two phases; 
only observations on the cluster, h and x Per, 
suggest that S Per (large amplitude variable) 
may lie at the tip of the red supergiant branch 
in the HR diagram because it appears redder 
than the small amplitude supergiants of the 
same cluster (Smith, 1976). 

Different Time-Scale Light Variations of Ir- 
regular Variables. Visual or photoelectric pho- 
tometry monitoring of a few irregular stars 
(e.g., VY CMa, and VX Sgr) indicates that their 
variations result more from chaotic or random 
processes than from harmonic processes. This 
could be interpreted as irregular surface 
features — spots, loops, or cells — that should 
have an effect on the light variation larger than 
have the pulsational effects. Spectroscopic 
observations in support of random variations 
in these stars appear in literature. As an exam- 
ple, we quote the variability in the emission 
lines in the IUE low-resolution spectra of the 
Lb carbon star, TX Psc, over a 3-year period, 
reported by Baumert and Johnson (1984). 
These authors show a probable time variation 
of the Mg II X2800 lines by a factor of 3 to 25 


and the possible variation of the C II X2325 line 
by 25 percent in the spectra of this object. On 
the basis of these data, Avrett and Johnson 
(1983) attempted to model a chromosphere for 
cool carbon stars (see de la Reza, this volume). 
A monitoring during 4 years (1981-1984) was 
jointly developed by the Bloomington and 
Toulouse groups on the low-resolution and 
long-spectral observations of TX Psc (Johnson 
et al., 1986). The intensity variations of Mg II, 
C II, and Fe II emission lines were confirmed 
with a factor of variability of 6 to 15, 5, and 
8 to 80, respectively. Moreover, all the lines ap- 
pear to vary together. The possible mechanisms 
for the production of these emission lines are 
reviewed: shock waves such as in Miras, mag- 
netic fields associated with plages, heating by 
short-period acoustic waves, and changes in ab- 
sorbing overlying clouds by their own motion. 
However, the continuous flux in the long-wave- 
length spectral range (LWR) is essentially con- 
stant. A recent high-resolution LWP spectrum 
of TX Psc confirms unusual profiles and line 
strengths for h and k lines. These lines are 
heavily absorbed by overlying Mg II, Mn I, and 
Fe I absorption features (Eriksson et al., 1985). 

As in semiregular variables, changes in light- 
curve shapes are observed: 

1 . Some irregular stars, like W CMi, have 
an irregular light curve, but during some 
periods of time, their light curve suggests 
a regular light period and during other 
periods a very constant brightness 
(Krempec, 1973). Some other stars have 
periodic variations during some intervals 
of time, and shortly thereafter, are not 
variable at all. BM Gem is a good exam- 
ple pointed out by Krempec (1973); RX 
Lep is the typical representative star. 

2 . Some stars have rapid changes of 
brightness. Krempec (1973) notes a 
change of 0.80 mag during 30 days on 
UY And and “very rapid light variations 
of the order of a few tenths of magnitude 
during a few days” on the light curve of 
SV Cyg. 
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3 . Other stars, such as TT CVn, have a 
period approximately defined at different 
epochs (Wendel, 1913; Krempec, 1973); 
however, during some time intervals, the 
brightness varies rather irregularly 
(Krempec, 1973). 

Ashbrook et al. (1954) describe the light 
curve of the SRc-Lc star, fi Cep, by a second 
stochastic chain. “The light variation is not ex- 
plicable by a simple pulsation; rather, it may 
be interpreted as arising from temporary, ran- 
dom surface disturbances on the star.” Sharp- 
less et al. (1966) found that the light curve is 
the result of the sum of periodic terms, with 
4836 days of fundamental period. Its light 
variations have short-period components and 
at least four different periods between 700 and 
14000 days (Polyakova, 1983). By spectral anal- 
ysis, Mantegazza (1982) shows that the light 
variation of fi Cep may be explained by the 
superposition of two periodic terms and their 
nonlinear couplings, the longer term having a 
cepheid-like shape. The results of Mantegazza 
tend to support the second hypothesis of Sto- 
thers and Leung (1971; see below) and are in 
agreement with the theoretical ones obtained 
by Uus (1976), who demonstrates that some red 
supergiants are pulsating in the radial fun- 
damental and first-overtone modes. 

In general, the Lb and Lc stars have not 
been sufficiently observed. To have a global 
point of view of their variations, some of them 
must be followed individually during many 
years. However, the analysis of the light curves 
of many of them should give indications to 
ascertain if chaotic variations are really possi- 
ble in the Lb stars. More observations are 
urgently needed to confirm this suggestion. 

Finally, we point out the observations of 
flares in these stars. In the chaotic point of view 
of the irregular star variations, flares should 
represent the shortest time variations. However, 
they are not really used in the analysis of the 
star behavior because, at present, a permanent 
followup is impossible with the conventional 
telescopes, photoelectric equipment, and adapt- 
ed filters, which prevent the study of con- 


secutive flares and their correlation with other 
variations. As an example of flares on such 
variable stars, we describe the recent event 
recorded on ^ Cep by Arsenijevic (1985). A sud- 
den fast (20-second duration) brightening of 
0.034 magnitude on V was observed on August 
3, 1981, at Oh 31.9 mn UT, with a standard 
deviation in the measurements of ± 0.01 1 mag; 
simultaneously, a shallow minimum in the 
polarization position angle appeared. The 
declining phase was 28 mn long. The similari- 
ty of this phenomenon with the red-dwarf flares 
strongly favors a stellar flare in the supergiant, 
fi Cep, although the declining phase was longer 
than in red dwarfs, maybe being correlated with 
a larger amount of irradiated energy. The flare 
appeared during the minimum light, such as in 
the case of the solar ones which are due to a 
large number of spots at the stellar surface. The 
total energy irradiated by the ft Cep flare is 
about 10 32 joules, fitting flare energies observed 
in the Sun and in stars. 

We also note flare activities in the radio 
wavelength range, for example, in a Ori and 
a Sco (Hjellming and Gibson, 1980; see also 
the section on Photometric Observations ); they 
infer the presence of magnetic fields (see M. 
Querci, this volume). 

To conclude the discussion on irregular 
stars, we mention the light variability of the 
brightest and the most famous irregular var- 
iable star, which has been observed more or less 
regularly over 60 years— a Orionis. 

The M2 lab supergiant, a Ori (Betelgeuse), 
was monitored by photoelectric photometry in 
the visual quite regularly from 1916 to 1931 by 
Stebbins (1932). Occasionally, Johnson et al. 
(1966) and Rucinski (as quoted by Goldberg, 
1984) have also made some observations. In re- 
cent years, Krisciunas (1982) systematically 
observed a Ori in B and V filters from 1979 to 
1982, and Guinan (1984) began measurements 
(around Ha and in the blue at X4530) in 1981 
that have continued to the present. Besides 
these data, AAVSO has recorded regular mea- 
surements for the last 60 years, and this exten- 
sive temporal coverage makes them very useful, 
especially in the absence of professional works. 
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However, these data must be viewed with great 
caution because the scatter is very large in a 
single measurement. All the above-quoted data 
are shown on Figures 1 and 3 through 5 in 
Goldberg (1984) and on Figure 2 (upper panel) 
in Guinan (1984). The characteristics of the a 
Ori light curves, particularly clearly seen on 
Guinan’s Figure 2, are long-term light varia- 
tions with a period of about 5.78 years, and 
with a mean amplitude of about 0.4 mag, and 
short-term irregular variations with rise or 
decline time scales of a few weeks to some 
months superposed on the long-term periodic 
light variations. The following systematic 
trends are noted: 

1.1916-1931: Stebbins (1932) found a 
period of 5.4 years after the use of a dras- 
tic smoothing process which cancels all 
the short-term fluctuations (see his Fig- 
ure 1). From radial velocity data, Jones 
(1928) and Sanford (1933) obtain a per- 
iod of 5.78 years. 

2. 1934-1938: The photoelectric data of a 
Ori show no evidence of periodicity, the 
points scatter widely, and the radial ve- 
locities are not connected to the 6-year 
period throughout the period of Adams’s 
(1956) observations: 1937-1947. 

3. 1939-1965: The brightness variations 
seem to follow the 5.78-year period and 
between 1949 and 1961, a phase lag 
relative to the mean curve previously 
found by Stebbins (1932) seems to be 
observed. 

4. 1965-1975: Random variations on 
brightness and radial velocities appear 
again. 

5. 1975-1980(7): The 6-year period was 
followed surprisingly well by the radial- 
velocity measurements of Boesgaard 
(1979) and Goldberg (1979). 

Moreover, short-lived abrupt drops in 
brightness are evident in 1926, 1934, 1942, 
1943, 1947, 1951, 1971, 1974, and 1981 and 
very short time variations appear in the 10-/tm 


region; the intense infrared emission surround- 
ing a Ori varies in shape, size, and intensity 
from night to night during September and 
December 1966. The stars fi Cep and a Tau 
have the same behavior (Low, 1965). 

As noted by M. Querci (this volume), the 
atomic lines giving the radial velocities are not 
listed in the text of some of the papers men- 
tioned above. Consequently, we suppose that 
the lines used are from various excitation poten- 
tials (i.e., from atmospheric layers of different 
behavior). Therefore, many possible correla- 
tions between the radial-velocity variations and 
drops or declines of brightness in large band- 
width filters, like UBV, cannot be fully inferred 
and exploited because we are unable to locate 
the layers from which the light variations come. 
Some examples of such unexploitable data are: 

1 . Between Julian days 2424538 and 
2424598, the radial velocity decreases by 
nearly 4 km/s, while the brightness in- 
creases by 0.2 mag in the same time 
(visual brightness). After 200 days, the 
radial velocity was found to be lower by 
3.5 km/s and the brightness was also 
declining during this time. 

2 . On Julian day 2431400, a Ori suffers a 
steep drop in radial velocity, followed 
some days later by a fall in brightness by 
about 0.6 mag. Goldberg (1984) argues 
that this rapid change takes place soon 
after the minimum of the “mean veloc- 
ity curve” like the changes seen in 1925, 
and shows that some other large distur- 
bances follow the two velocity minima of 
1972 and 1978. 

During the short-term variations, the change 
on the linear polarization of the visible light 
variation varies with the same time scale as the 
brightness. This suggests a high degree of asym- 
metry and a local origin in the features involved 
in these variations (Tinbergen et al., 1981). 
Hayes (1981) suggests that the change of po- 
larization follows the formation and growth of 
the local features on the surface of the star 
(lower chromosphere) and progressive changes 
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in their orientation. Viewing the entire phenom- 
ena described by Goldberg (1984), it seems 
“equally possible that a sharp drop in radial 
velocity was accompanied by an initial bright- 
ening followed by rapid fading,” which can be 
explained by a gaseous matter ejected from the 
star which diffuses and condensates into grains 
and becomes optically thick in the visible region 
(Goldberg, 1984). 

To conclude, it appears that the long-term 
light variations probably arise from the radial 
pulsation of the star, the short-term variations 
having been believed to be caused by large-scale 
photospheric convective cells (Schwarzschild, 
1975). Hayes (1982), Goldberg (1984), and An- 
tia et al. (1984) favor such an explanation, the 
two first authors for explaining polarization 
changes. However, Roddier and Roddier (1983, 
1985), Petrov (1983), Roddier et al. (1984), and 
Karovska (1984) argue against it (see M. Quer- 
ci, this volume), and Guinan (1984) notes a lack 
of correlation between the short-term variations 
in the TiO band strength and the short-term 
light changes. This fact suggests that “the 
mechanism responsible for the short-term light 
enhancement is not linked to temperature in- 
creases as would be expected from ascending 
giant convective cells... and that the brightness 
enhancements are produced above the stars’s 
surface.” 

The problem might be resolved if the recent- 
ly claimed companion objects to a Ori are con- 
firmed. That circumstance would also seem to 
explain some variations in chromospheric emis- 
sion-line radial velocities, as well as circum- 
stellar dusty clumps (see M. Querci, this vol- 
ume). (Anecdotally, Bottlinger (1910) mentions 
a Ori as a binary star for which he finds a 
6-year period. Other authors (Wilson and 
others) as reported by Karovska (1984), also 
classified a Ori as a binary star.) 

Eruptive Variables: The RCB Stars 

Generalities. The main characteristics of the 
eruptive RCB stars are the abrupt drops in 
brightness, followed by a longer climb back to 
normal light, probably due to the sudden ejec- 


tion of highly absorbent matter. They are oc- 
casionally thought to be the progenitors of type 
I supernovae (Wheeler, 1978), novae (Mc- 
Laughlin, 1935), planetary nebulae (Webster 
and Glass, 1974), or helium stars (Tutukov and 
Iben, 1985), and white dwarfs (Schonberner, 
1977). 

In 1795, E. Piggott discovered the variability 
of R CrB, the prototype of the well-known 
RCB stars. In the southern hemisphere, a few 
random observations of RY Sgr have been 
made since 1751 (Lacaille, 1847). From 1895 to 
1908, the brightness of UV Cas dropped twice 
by 1 to 1.5 mag (Florya, 1949); its variability 
was recognized with the 4-mag fall of 1913 
(Shenavrin, 1979). Innes(1903, 1907) published 
the first series of observations of some RCB 
stars. After monitoring R CrB for several years, 
Sterne (1935) concluded that this star is a 
“perfect irregular” with its visual variation 
from 5.8 at maximum to 14.8 for the deepest 
minimum, whereas the light of SU Tau falls 
almost annually. 

As in the other variables, the spectral type 
of the RCB stars changes during their bright- 
ness variations, and accordingly, changes in the 
spectra and colors of the RCB stars are noted. 
The latter are extremely complex because the 
star waxes and wanes. A general behavior of 
the stellar spectrum during one of the abrupt 
drops in brightness and recovery may be 
described as follows, mainly on the basis of R 
CrB and RY Sgr observations. 

At maximum, the spectrum is very similar 
to that of a supergiant of spectral type F or R, 
with narrow and sharp absorption lines, the 
main differences being in the weakness or 
absence of the Balmer lines and in the strong 
lines due to carbon. The cooler stars have 
strong molecular absorption bands of CN and 

C 2‘ 

As the star dims, we continue to see the 
absorption spectrum, but it is veiled. However, 
in general, the molecular spectrum becomes 
more and more intense as the minimum is 
reached. There is a sharp emission spectrum 
which appears during the fall and suffers 
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gradual changes. Some of these narrow emis- 
sion lines sometimes replace the atomic absorp- 
tion lines (Herbig, 1949). First, the sharp lines 
of Fe II appear (Payne-Gaposchkin, 1963). As 
the fading continues, the spectrum becomes 
dominated by Ti II lines. Occasional anomalies 
are detected. For example, Alexander et al. 
(1972) found the same intensity for the two Sc 
II lines at 4354.61 and 4320.75 A as in the 
laboratory. Spite and Spite (1979) showed that, 
in the next drop of RY Sgr, the 4354.61 A line 
has disappeared and the 4320.75 A line is ap- 
proximately 500 mA. In general, this sharp 
emission spectrum is slightly displaced to the 
blue (3 to 10 km/s) relative to the absorption 
lines at maximum. Its intensity increases rapidly 
during the initial part of the drop and then 
grows weaker as the star continues to fade. It 
is likely produced by an expanding region 
(chromosphere or large envelope), with decreas- 
ing expansion as the star approaches minimum. 
In the decline of RY Sgr during 1977, Spite and 
Spite (1979) note that the narrow emission lines 
and the absorption lines of the atmosphere itself 
have the same radial velocity. The behavior of 
this emission spectrum is apparently not related 
to the absolute magnitude of the star (Feast, 
1975), but seems to depend on the time from 
initial decline (Payne-Gaposchkin, 1963). Dur- 
ing the 1948 minimum of R CrB (Herbig, 1949), 
only the second part of the narrow emission 
spectrum was seen (Ti II). Perhaps this is con- 
nected to the observed low initial rate of 
decline. 

A very broad emission spectrum appears 
some days after the brightness begins to fade 
(40 days in the case of R CrB in 1948 and 70 
days in the case of RY Sgr). In the spectrum 
of R CrB, the lines of He I 3889, Ca II H and 
K, and Na I D appear with different widths 
(Rao, 1975). Using narrowband filters, Wing 
et al. (1972) discovered the presence of a strong 
He I 10830 emission line during the minimum 
of R CrB. As the minimum progresses, the pro- 
file of these lines changes, and the intensity re- 
mains constant at first, then increases, and 
finally decreases. In the spectrum of the cooler 
RCB stars, these broad emission lines are 


sometimes accompanied by absorption com- 
ponents with velocities of -250 to -130 km/s 
(Querci and Querci, 1978; Rao, 1980b). (See M. 
Querci, this volume, for details on observations 
of the He I 10830 A line in red giants and the 
implied constraints to future modeling.) Be- 
cause the UV spectrum fades much less than 
the visible and red spectrum, additional sources 
of emission shortward of 4000 A have to be 
taken into account. As the star approaches 
minimum light, the relative intensity of the 
emission lines of each element changes gradual- 
ly (implying variations in the conditions of ex- 
citation). Also, a variable polarization with 
variable wavelength dependence is observed in 
stars like R CrB and RY Sgr (Serkowski and 
Kruszewski, 1969; Coyne and Shawl, 1973). At 
the end of minimum phase, P Cygni profiles 
develop (Payne-Gaposchkin, 1963; Querci and 
Querci, 1978; Rao, 1980b). The bands of CN 
are observed in emission at this time (Wing et 
al., 1972). Seeds and Ignatuk (1973) have 
shown changes in strength of the C 2 emission 
bands in the R CrB minimum. These changes 
are probably due to the filling in of the 
photospheric absorption bands by chromo- 
spheric emission. Taking into account the nu- 
merous atomic lines in the violet, the chromo- 
spheric contribution to the total observed light 
increases to shorter wavelengths; consequently, 
the dependence of polarization with wavelength 
is not linked to the particle size only. 

As the star brightens again, the absorption 
lines gradually come back, while the emission 
spectrum fades with a varying rapidity in the 
various lines. Finally, the maximum spectrum 
is reached approximately 2 magnitudes before 
the light maximum itself; in R CrB, Rao et al. 
(1981) observe h 3 and k 3 Mg II components, 
which demonstrate that this star also has a con- 
tinuous mass loss during the maximum light 
through a permanent chromosphere (Payne- 
Gaposchkin, 1963; Rao, 1974, 1975). 

Can we draw a coherent qualitative model- 
ing of the RCB phenomenon? 

Many models have been proposed during 
the past 20 years. Payne-Gaposchkin (1963) 
suggested that the particles are formed in the 
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upper photosphere of the star, because the 
emission lines observed at minimum light are 
not themselves obscured. Wing et al. (1972) 
propose that there is a quasi-permanent blotchy 
cloud orbiting around the star. Humphreys and 
Ney (1974) associate the particle cloud with the 
atmosphere of a probable cool LPV compan- 
ion. These last two models are attractive, but 
pose problems in so far as they fail to explain 
why the decline is rapid and the rise is slower. 

Nowadays, it appears that the initial 
hypothesis of Loreta (1934) and O’Keefe 
(1939), in which the light minimum of the RCB 
stars is due to dust obscuration, gives the best 
account of the observed features. Feast (1979), 
sustaining this basic model of ejection of thick 
dust clouds, explains how it might act. The gas 
is ejected at the top of one or some of the large 
convective cells, and it crosses the deeper layers 
of the star’s atmosphere. This ejection is made 
radially through a fairly large area of the stellar 
surface in a semiangle of -20° and roughly at 
20 km/s. If the material is ejected along our 
line of sight, a major visual luminosity 
minimum occurs, whereas, if it is ejected at an 
angle to our line of sight, it causes a minor 
minimum (e.g., Forrest et al., 1972). The gas 
expands and cools, and the graphitic carbon 
condenses. The resultant dust clouds expand 
and cause an eclipse which allows the 
chromospheric sharp lines to be seen. The par- 
ticles move away by radiation pressure. When 
the new cloud collides with the circumstellar 
shells, strong and broad emission lines are pro- 
duced. Because the speed of the collision is 
greater than the escape velocity, the matter is: 
(1) partially ejected in space with a drop in den- 
sity and in optical depth of the medium caus- 
ing the emission lines to fade, (2) partially 
replenishing the circumstellar patchy envelope 
of dust and gas, giving rise to the infrared ex- 
cess observed in some RCB stars (e.g., in RY 
Sgr and R CrB). Recent IUE data confirm the 
Loreta-O’ Keefe cloud ejection model and are 
consistent with the dust to be composed of car- 
bon (Hecht et al., 1984). In addition, they show 
that 5 to 60 nm glassy or amorphous carbon 

rather than graphite is formed around the RCB 
stars. 


It is now worth turning to some observations 
that argue in favor of this model and complete 
it. 

Patterson et al. (1976) showed that their 
spectrophotometric measurements made on dif- 
ferent nights during the R CrB climb in 1974 
were consistent with extinction of particles that 
they thought to be spherical graphite ones of 
about 0.07-/un radius, expanding by radiation 
pressure. No variations of the particle size were 
detected. 

Orlov and Rodriguez (1974, 1981) find a 
microturbulent velocity of 8.9 km/s in XX Cam 
and 1 1 km/s in UV Cas, which supports the 
idea suggested by Howarth (1976) that the in- 
itial condensation may be due to turbulent mo- 
tions (seed creation). 

Forrest, as quoted by Rao (1980a), points 
out that the extinction during the decreasing 
light tends to be neutral without any color 
changes, but that the colors redden during the 
climb back to maximum. Using these confir- 
mations and the electron pressure computed 
from the broad emission lines, Rao (1980a) sug- 
gests that the gas that produces these broad 
emissions also produces the neutral extinction 
by electron scattering during the decreasing 
light. The grains appear later, and the cir- 
cumstellar shell is replenished by ejecta at 
minimum light. Consequently, the light from 
the photosphere is obscured, and the 
chromosphere (Payne-Gaposchkin, 1963; Alex- 
ander et al., 1972) or the circumstellar gas shell 
(Hartmann and Apruzese, 1976) is the obvious 
source of the narrow emission spectrum. As 
pointed out by Herbig (1949), the emitting gas 
has a relatively low level of excitation. 

The radial velocity of the chromosphere is 
the same as that of the star itself (Spite and 
Spite, 1979). When the emission region is sud- 
denly cut off from its source of excitation, it 
would be an ideal place for the production of 
an electron recombination spectrum. Feast 
(1969) suggests that the CN molecule is 
reasonably abundant in these stars and that its 
spectrum has a longward edge at 4000 A. This 
could explain the blue and UV continuum 
anomalies by: CN + e — CN“ + Ivy. The 
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emission excess on the V band observed by 
Shenavrin (1979) on UV Cas tallies with this ex- 
planation. Moreover, this model supports the 
changes in the wavelength-dependent polariza- 
tion, showing that different particle sizes exist 
successively. 

As noted above, some RCB stars have an 
excess of infrared. The first observations are 
recorded by Stein et al. (1969), who find the 
infrared emission much brighter than the ex- 
pected emission from the star itself, and who 
conclude that 40 percent of the star’s total 
luminosity is in the infrared; Shenavrin (1979) 
finds 60 percent for SU Tau. An excess of in- 
frared is also found by Lee and Feast (1969) 
in RY Sgr. However, XX Cam and UV Cas 
have no infrared excess at all (Shenavrin, 1979; 
Rao et al., 1980), and sometimes R CrB does 
not have any IR excess during some minima, 
as in 1972 (Rao et al., 1980). Either the physical 
conditions are not favorable for grain forma- 
tion in XX Cam and UV Cas or these stars are 
between the RCB and the nonvariable HdC 
stars. To explain the infrared emission, Herbig 
(1949) believes that grains have to condense far 
from the stellar surface, whereas Maron (1974) 
suggests the following way: the freshly ejected 
particles could have the form of “Platt” par- 
ticles of 3 to 30 A; they would absorb the visi- 
ble radiation better than graphite particles but 
would not reradiate in the infrared. Then, these 
particle grow by accretion up to sizes between 
10‘ 6 and 10“ 5 cm (the classical grain nuclei 
size). Therefore, the occurrence of the infrared 
emission waits for the growth of the particles. 
This growth period is a function of the physical 
conditions in the shell, mainly low density and 
velocity of escape, which vary from minimum 
to minimum. 

The broad emission lines which are observed 
after the beginning of the decline are due to 
matter which is ejected at high velocity and col- 
lides with the circumstellar envelope material 
with a phase lag of 30 to 70 days. From the 
radial velocity of the circumstellar lines and this 
phase lag, we conclude that the shell of gas and 
grains is at about 4 to 8 A. U. from the star. 
This agrees more or less with Lee and Feast 


(1969), giving 5 to 10 A. U., and with Pecker 
(1976), giving a dust-cloud envelope of approx- 
imately 350 R, at the IR maximum brightness 
to 600 R, at the minimum. Moreover, Pecker 
estimates the total mass of the envelope to 
5.10 -7 solar mass. The broad emission lines are 
noted violet-displaced by Payne-Gaposchkin 
(1963). However, during the end of the decline 
of RY Sgr observed by Spite and Spite (1979), 
the broad emission lines seem slightly displaced 
to the red relative to the star. (They extend from 
-100 to +240 km/s.) This last observation is 
explained by an ejection of matter in which only 
a small amount is ejected toward the observer 
and a larger one is ejected backward. Wdowiak 
(1975) assumes that 10 to 100 convection cells 
are responsible for the formation of dust blobs 
that are ejected by radiation pressure. Rao 
(1980a) indicates that the broad emission lines 
are formed in a gas with Ne = 10" - 10 12 per 
cm 3 and Te » 10 4 K. 

High-velocity shells and/or clumps are 
observed in MV Sgr at -200 km/s with a 
temperature of 800 to 900 K (Krishna Swamy, 
1972), although the star has a hotter shell about 
1300 K (Feast and Glass, 1973). 

The presence of the circumstellar shell in 
some RCB stars is also pointed out by IRAS. 
It discovers a bimodal distribution of RCB 
stars: (1) the majority of the RCB stars have 
a flux such as F v ~ v 1,6 ; they have an IRAS 
spectrum that can be understood as the spec- 
trum of a star with a constant ejection rate of 
dust where the emission is proportional to the 
frequency (Schaefer, 1985); (2) a fifth of RCB 
stars exhibit a qualitatively different IRAS spec- 
trum with a flux such as F v ~ v" 1 and a peak 
flux at 100 /im or beyond; they are brighter than 
the majority of the RCB’s. Their IRAS spec- 
trum is undistinguishable from normal plane- 
tary nebulae (known to have a large amount of 
dust at a great distance from the central star 
(Schaefer, 1985)). This fact is in agreement with 
one of Iben and Tutukov’s (1984) evolutionary 
scenarios that place the RCB stars in planetary 
nebulae. 

Light Curves. Surprisingly, the RCB stars are 
the only class of eruptive variables for which 
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maximum brightness is the normal state. The 
duration of this state is variable (up to several 
years). In addition, epochs exist in which the 
full maximum is not reached. For example, 
from 1898 to 1948, XX Cam had only one 
minimum in 1939-1940 (Yuin, 1948); R CrB 
gave flat maxima from 1925 to 1935, 1936 to 
mid- 1938, 1953 to 1956, 1969 to 1972, etc. The 
latter also shows a series of minima without full 
maxima being reached in the 1860’s (Mayall, 
1960b) and from mid-1962 to 1965. Howarth 
(1977) defines fades and minima, respectively, 
“as an initial drop of one magnitude from max- 
imum” and “a fall of one magnitude relative 
to local sub-maxima.” The events appear to 
obey Poisson’s statistics. For R CrB and SU 
Tau, the mean intervals between minima are 
532 ± 57 and 625 ± 104 days, respectively, and 
between fades, they are 1026 ± 156 and 1140 
± 220 days, respectively. The cooler star, S Aps 
(R3), has a mean time between fades of 1249 
days (Howarth, 1976). This should suggest a 
trend of decreasing activity with later spectral 
type. However, observations of more stars are 
needed to confirm this correlation. 

The speed of the decline has also been in- 
vestigated. In R CrB, Oberstatter (1972) points 
out a drop of 0.5 visual magnitude per day dur- 
ing the 1972 decline; amateurs have observed 

0.3 visual magnitude per day during the decline 
in 1983 (Proust and Verdenet, 1983), while the 
“normal” fade is about dm/dt * 0.1 mag per 
day (Howarth, 1976). 

The amplitude, the frequency, and the dura- 
tion of minima are unpredictable. The declines 
are much more abrupt than the rises to the max- 
imum. However, the shape of the decline and 
the rise is not the same for each minimum of 
the RCB stars; sometimes, “a star may rise 
halfway to maximum only, to fall again to an 
even deeper minimum” (Howarth, 1976). This 
could be caused by many ejection centers. 

Superposition of Light Curves. Different time- 
scale variations are very well defined in the 
RCB variable stars. They are: 

1 . The unpredictable decline of many 


magnitudes (already described) called 
“obscurational” minimum. 

2. The semiregular oscillations with a 
visual amplitude of 0.2 to 0.4 magni- 
tudes and with a pseudoperiod between 
19 and 120 days according to the star, 
called “pulsational” oscillations. 

3 . The very short nonpermanent oscilla- 
tions with a time scale of 1 or 2 hours 
detected in some RCB stars and prob- 
ably originating from the star itself. 

4 . A long-term periodicity in the infrared 
excess, perhaps due to a natural pulsa- 
tion of the circumstellar dust shell. 

Semiregular oscillations over a 38.6-day 
period were detected by Jacchia (1933) on RY 
Sgr. Many observations of RCB stars reveal 
that these small fluctuations appear at each 
phase of the variation cycle (Bateson, 1978). 
Mendoza (1978) detected them photometrical- 
ly in Ha and OI (X7774) filters. Pugach (1977) 
concludes from observations of RY Sgr that 
these pulsational variations do not depend on 
light fading. First, these oscillations are observed 
during the brightest, largest observable phase 
(i.e., the maximum). Fernie et al. (1972) report 
a visual variation of R CrB with an amplitude 
of 0. 15 mag and a period of 45 days. Totochava 
(1973a) confirms this amplitude variation, but 
deduces a variation period of 40 days from her 
1971-1973 observations. Alexander et al. (1972) 
confirm the oscillations of RY Sgr with a 
0.5-mag average amplitude and a 38.6-day 
period observed by Jacchia (1933); the light and 
color curves (Figure l-23a, b) show smooth 
variations: AV « 0.5 A (B-V) « 0.3 A(U-B) « 
0.5 mag. Sherwood (1976) finds tentative 
periods from 19 to 54 days for some RCB stars 
and a period greater than 90 days for S Aps, 
whereas Waters (1966) suggests 120 days and 
0.3 mag of amplitude variation for this star. 
Later, Kilkenny and Flanagan (1983) find 
evidence for a rapid period decrease, and 
Kilkenny (1983) shows that S Aps changes from 
the 120-day periodicity (present in 1960) to one 
near 40 days (around 1971). Bateson and Jones 
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Figure 1-23. R Y Sgr (a): V light curve for 1968 plotted against JD. The open triangle is an observa- 
tion from Lee and Feast (1969). The predicted maxima and minima in V are shown at the top of 
the diagram, (b): (U-B), (B-V), and spectral classification for 1968 plotted against JD. The same remarks 
as stated previously apply for the open triangles, maxima and minima. For the spectra, absorption 
(C) types are shown by filled circles (from Alexander et al., 1972). 


(1972) find such semiregular variations in two 
stars, UW Cen and GU Sgr, with 0.5 mag and 
42 days and 0.4 mag and 38 days of amplitude 
and pseudoperiod, respectively. (See the section 
Irreversible Changes in HdC Stars for details.) 

During the deep decline and the rise to max- 
imum, R CrB also has semiregular oscillations 


with a time scale of 40 days (Oberstatter, 1972). 
Isles (1973) confirms that this star occasional- 
ly displays oscillations, the amplitude of which 
decreases as the star brightens. Another RCB 
star, RY Sgr, shows oscillations during the 
climb back to maximum, similar to those de- 
tected during its maximum; these oscillations 


72 



support the idea that no major physical changes 
took place in the star itself during the minima. 
The greatest amplitude of the UW Cen varia- 
tions seems to be when the star is brightening 
(Bateson and Jones, 1972); it persists along the 
light curve except during the sudden falls to 
deep minima (Bateson, 1972). During maxi- 
mum light, XX Cam varies with the same amp- 
litude as that of R CrB, but the period of the 
former is more or less shorter than those of the 
latter (Totochava, 1973b). The light, color, and 
radial velocity of the RCB’s demonstrate that 
these stars pulsate. 

Research on light fluctuations with a smaller 
time scale requires continuous observations of 
the light. Herbig (1967) reports that Miskin ob- 
served rapid oscillations in R CrB near mini- 
mum, but very high frequency pulsations in R 
CrB were not detected by Horowitz et al. (1971) 
in 1969. Interesting observations were also 
made by Totochava (1975) on XX Cam from 


1972 to 1974. At first, Totochava conducted her 
observations with a single-channel photometer, 
alternatively in U, B, and V. Some variations 
appear in the yellow filter, with an amplitude 
of about 0.1 mag over several minutes. Fluc- 
tuations alternate with quiet stellar phases (e.g., 
no light variations are observed during the 1973 
autumn). The observations are simultaneously 
obtained with a three-channel spectrophotom- 
eter in the following ranges: 3350 to 3650 A, 
4155 to 4280 A, and 5120 to 5320 A. The ob- 
served curves obtained in this way show several 
synchronous fluctuations in the violet and blue 
regions during certain nights (Figure l-24a, b). 
The amplitude of the UV fluctuations are a lit- 
tle larger than in the blue; they amount to ap- 
proximately 0.1 mag and last for about 2 hours. 
During a period of several nights, the fluctua- 
tions become much less visible in all the col- 
ors. At the quiescent phase of XX Cam, 
Kolotilov et al. (1974) find interstellar polariza- 
tion characteristics, but a higher value of Pthan 
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Figure 1-24. Light curves of XX Cam obtained from simultaneous three-color observations (see text): 
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73 


for neighboring stars, suggesting that an intrin- 
sic polarization of the star makes some 
contribution. 

Following R CrB itself, Humphreys and Ney 
(1974) and Strecker (1975) observed a smooth 
variation between 1 .4 and 3.0 mag at L (3.5 /un) 
with a pseudoperiod of about 1 100 days. (The 
monitoring was done only between JD 2440000 
and 2442360.) In RY Sgr, besides the 38-day 
variation, Feast et al. (1977) also detected a long 
period variability at L, which is related to the 
variation of the star itself. They found a long- 
term change of 1.5 mag, but the period was not 
well analyzed for lack of observations. Feast 
and Glass (1973) and Feast (1975) show how 
the star can make only a minor contribution at 
L, and Feast et al. (1977) testify that 70 to 80 
percent of the L flux is due to the infrared ex- 
cess (in the case of RY Sgr), with an estimated 
dust shell radius of 5 to 10 A. U. (Lee and Feast, 
1969). “It seems natural to suppose that this 
is due to pulsation of the circumstellar dust shell 
itself” (Feast et al., 1977). 

It is useful to remember that Bergeat and 
Lunel (1980) state that “the contribution of 
shell emissions to the J values is usually small” 
for the Miras, SR, and L stars, while Walker 
(1980) concludes that there is no evidence of 
dust thermal emission in the near infrared for 
most of the N-type irregular variables, and 
Tsuji (1981c) uses the L flux for his 
temperature-effective determination of SRb 
and Lb carbon stars. Consequently, although 
the C spectral type Mira, SRb, Lb variables, 
and the RCB stars are all carbon stars, the ef- 
fectiveness of the circumstellar dust shells is 
very different at L wavelength among them. 

Additional observations are required to im- 
prove our knowledge of the RCB circumstellar 
dust shells. For example, the phase lag between 
V, 1104, and L magnitude curves could help us 
to define the thickness of the envelope. Some 
cool variable stars with extensive dust shells 
have already shown such long-term variability. 

Correlations Between Observed Quantities. The 

observations of the RCB stars at different 
epochs indicate that, in general, the infrared ex- 


cess compensates for the loss of light in the visi- 
ble and UV wavelengths (Alexander et al., 
1972); however, this is not true for all the stars. 

The oscillations of brightness of RY Sgr dur- 
ing the climb back to maximum, and during the 
maximum itself, are correlated with the radial- 
velocity variations of the absorption lines (Fig- 
ure 1-25) and have 38.6 and 39.0 days of period, 
respectively (Alexander et al., 1972), confirm- 
ing that RY Sgr pulsates as a helium star of two 
solar masses. 

The emission lines of RY Sgr do not have 
the same behavior pattern as the absorption 
lines (Alexander et al., 1972). We are not able 
to decide whether the emission lines follow the 
38.6-day pulsation or not; large differences on 
radial velocity between absorption and emission 
lines are nevertheless observed. 

Pugach (1977) found that the beginning of 
the obscurational decline appears between 
phases 0.24 and 0.37, and he chose the origin 
<p 0 of the pulsational phase <p at the minimum 
light of pulsation. If we suppose that these pul- 
sational oscillations are roughly sine waves and 
if we change the origin <fi 0 to the maximum, we 
see that the beginning of the obscurational 
decline is around the phases 0.74 and 0.87. This 
reminds us of the famous phase 0.8 of the Mira 
variables. Can we conclude that the obscura- 
tional decline is a result of one pulsational 
oscillation being stronger than the others? 
Pugach (1977) concludes that “the serial 
number of cycles when declines occur, is prob- 
ably accidental,” which seems erroneous today. 
For the decline of RY Sgr in 1967, the pulsa- 
tional phase of the drop of light does not con- 
firm the Pugach phase. Howarth (1976, 1977) 
notes that “the 1967 fade of RY Sgr took place 
at the time of minimum of the secondary varia- 
tion.” 

The occurrence of pulsational minima in RY 
Sgr can be well represented by a quadratic solu- 
tion of the cycle number n (JD min = T q + nP g 
+ n 2 k). The (O-C) residuals are shown in Fig- 
ure l-26a. The value of k is —(-51 ± 2) X 
10' 5 day per cycle. This value is in agreement 
with theoretical models for the evolution of 
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Figure 1-25 . /?Y Sg/* /ty: K light curve and radial velocity (km/s) curve during the rising of 1969; 
(b): during the maximum of 1970. Predicted maxima and minima in V and predicted maxima in radial 
velocities are marked above the curves (from Alexander et a l 1972). 


HdC stars (Schonberner, 1977). The compari- 
son of the RY Sgr light curve (Figure l-26b) and 
the depth estimates of pulsational minima ob- 
served on the same time scale (Figure l-26c), 
shows that the deeper pulsational minima tend 
to occur on or near rising branches of the 
obscurational minima. 

Bateson (1972) has noted a similar effect on 
UW Cen. This behavior remains unexplained 
(Kilkenny, 1982). Forrest et al. (1972) were the 
first to stumble on one of the most important 
clues to the RCB phenomenon: when the star 
drops into the deep minimum, the infrared ex- 


cess is not affected (no change in the emitting 
flux of the dust shell and no change in the phase 
of the L curve; Figure 1-27). The envelope and 
the atmosphere of these stars seem to be 
dynamically disconnected. This restricts the 
model of obscuration, and Hartmann and 
Apruzese (1976) suggest that the dust is formed 
above a small region of the stellar surface and 
is blown away by the radiation pressure. Con- 
sequently, comparative studies of the visible 
and infrared light curves come to the same con- 
clusions as those of the foregoing investiga- 
tions. Glass (1978) argues in the same vein 
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about the 1975/1976 visual minimum of S Aps: 
when J, mainly due to the star itself, drops to 
minimum, L is not so strongly affected. 

There is a pressing need to summon ob- 
servers’ help in making observations during 
light fades of different amplitude and duration 
to find correlations between observed features 
and to conclude which parameters are linked 
together and which are not. 

Nonvariable Stars 

After this extensive review of all the types 
of cool variable giants and supergiants and their 
probable connections with white dwarfs, 
planetary nebulae, or supernovae, we have to 
discuss stars which— though located in the same 
area of the HR diagram as the red variables — 
are surprisingly mentioned as nonvariable or 
normal stars. 


n 



Figure 1-26. RY Sgr (a): (O-C) residuals from 
the quadratic formula of the pulsational 
minima plotted against cycle n and Julian date 
for ephemeris with decreasing period; (b): 
sketch occurrence of obscurational minima on 
same time scale as (a); (c): estimates of depths 
of pulsational minima on same time scale as (a) 
(from Kilkenny , 1982). 
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Figure 1-27. SAAO observations of RY Sgr at L(3.5 fim) and J(1.2 fim) (from Feast , 1979). 
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Nonvariable HdC Stars. Knowing that the RCB 
stars have fluctuations of brightness around 0. 1 
mag or more with periods as large as 40 days 
(38.9 for RY Sgr and 44 for R CrB), Rao 
(1980b) searched for variability in the non- 
variable HdC stars. He observed the cool HdC 
star, HD 137613, eleven times in UBV during 
85 days and detected no periodic variations 
greater than 0. 1 mag in V . This variation range 
is larger than the expected RMS error ( ± 0.025 
in V, ±0.020 in B-V, and ±0.025 in U-B ), 
which may indicate a real variability in cool 
HdC stars. However, we conclude that this 
number of observations is too low to point out 
clearly any kind of variability. The total range 
of 0.1 mag in V encourages new observations 
with a time scale shorter than the week time 
scale used by Rao. As an example, such a long 
time scale should prevent us from discovering 
the short time-scale fluctuations in XX Cam by 
Totochava (1975), as well as the suspected 
nonradial-connected phenomena as discussed 
above. 

Feast and Glass (1973) demonstrate that the 
nonvariable HdC stars do not possess any IR 
excess, and Rao (1980a) points out that the non- 
variable star, HD 182040, has no envelope. This 
should be correlated with the low evidence for 
large mass loss on the nonvariable hot HdC 
stars (Schonberner and Hunger, 1978; Herber 
and Schonberner, 1980). The high-dispersion 
IUE observations of HD 182040 show no emis- 
sion lines (i.e., chromospheric indicators; 
Johnson et al., 1984). 

Helium Stars. Although all the following 
helium stars are known to be non variable: BD 
+ 37 °442, BD ±10°2179, HD 124448, HD 
160641, and HD 168476, Herbig (1967) stressed 
that possible long-period variations exist. Lan- 
dolt (1968) made intensive observations over 
periods of a few nights on two or three occa- 
sions each year for detecting short-term and 
middle-term variations; the set of stars are those 
given above, plus HD 264111 (Stephenson, 
1967a) and BD + 13 °3224 (Berger and Green- 
stein, 1963). Only short-term variations are 
noted on BD + 10° 2179: on February 3, 1968, 


the star changed by 0.04 mag in V between two 
consecutive nights, a change greater than the 
probable error of a single measurement. Rao 
(1980b) also analyzed the variability of this star 
with 13 observations made on 100 consecutive 
days. The star does not seem to have any 
periodic variations in light greater than 0. 1 mag 
in V. These variations are in agreement with 
those observed at the same epoch by Landolt 
(1973), who gives a total range of 0.09 mag in 
V. Here also, the observations are too few and 
lead to the same conclusions made for the low- 
temperature HdC stars. 

Observations like these made by Landolt 
(1973) must be carried out: the large scatter in 
the data of HD 160641 and the evident bright- 
ening of the star by -0.1 mag during 7 hours 
of observation demonstrate that short time- 
scale variations exist in helium stars as well as 
in the cooler HdC stars (Figure 1-28). On the 
other hand, Herbig (1964) concluded that: “A 
survey of the bright helium stars on a long-time 
series of ordinary patrol photographs would ap- 
pear to be the more promising approach.” To- 
day, we endeavor to collaborate with the 
IAPPP or GEOS members or to use automatic 
telescopes with automatic photometers or 
satellites to do such surveys. Correlations be- 
tween brightness and radial velocity are also 
seen in the extremely hydrogen-deficient hot 
star, BD + 13 ° 3224, by Hill et al. (1981). The 
correlations between observations and evolu- 
tionary models (Schonberner, 1977; Hill et al., 
1981) are to be analyzed carefully (Tutukov and 
Iben, 1985). 

Early R Stars. There are no systematic studies 
on the relative number of variables and non- 
variables among carbon (R and N) stars, such 
that it has been necessary to gather such infor- 
mation from catalogs. For example, among the 
122 carbon stars that are brighter than 10.0 mag 
in the DDO catalog (Lee et al., 1940, 1947; Lee 
and Bartlett, 1945), 68 percent are variable. 

Looking at the early R stars, we discover 
that only some are variables. (The variability 
is everywhere!) Among these early R variables, 
some, like BD +20°5071 and BD +69°417 
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Figure 1-28. Top: plot of the differential V magnitude for HD 160641 (the gap in 
data at JD-2441452.70 and 452.83 occurs at times when UBV standard stars were 
observed); bottom: V magnitude plot for the comparison star BD-17°4889. For the 
7-hour observing run, the large scatter in comparison star data is due to technical 
difficulties. The observations are made on May 15, 1972 UT (from Landolt, 1973). 


(RU Cam), are bluer than the nonvariables. RU 
Cam has a period of 22 days (W Virginis type). 
Its spectrum varies from RO at minimum to KO 
at maximum. BD + 57°2161 is also suspected 
to be variable, but it is redder than the non- 
variable R stars. A few further photometric 
analyses should be undertaken for a better 
knowledge of these stars. Among the R2 stars, 
BD + 19°3109 and BD + 02°3326 are variable 
and have a higher color index than the normal 
stars. 

Nonvariable M-Type Stars. Only a few papers 
are devoted to normal M stars. Generally, these 
stars are discussed in papers in which their 
characteristics are used as reference to the Mira, 
SR, or L variable stars. One of the most 
valuable proofs that we observe a normal M 


star is the nonvariability of the strength of ab- 
sorption molecular bands and atomic lines and 
of the adjacent flux peaks. 

The nonvariable M-type stars extend from 
MO to M8, with a luminosity class from lab to 
III. On the following stars, /3 And, 7r Aur, a 
Cet, n Gem, i Peg, 119 Tau, 83 UMa, and 8 
Vir, Maehara and Yamashita (1978) obtained 
repeated observations of their fluxes within 
about 0.05 mag for most wavelengths. The 
energy distribution varies gradually from class 
to class. For stars later than M5, the gradient 
of the 4000 to 5000 A region is mainly deter- 
mined by strong TiO absorption bands rather 
than by temperature effects on the flux peaks. 

We will not elaborate further on the non- 
variable M stars which are not especially within 
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the scope of this volume on variability and non- 
thermal phenomena. However, a red star can 
be accepted as nonvariable and, consequently, 
as a photometric standard, only if it is observed 
to be constant over a sufficient period of time. 
It often happens that red stars accepted as UBV 
standard stars (i.e., non variable a priori) prove 
to have highly variable chromospheric lines 
(Sherbakov, 1979). As a consequence, they can- 
not be considered as nonvariable stars in nar- 
rowband photometric systems (Wings, 1971; 
Querci and Querci, 1985b). 

IRREVERSIBLE CHANGES AND 
RAPID EVOLUTION 

Irreversible Changes in Miras — 

An Example: R Aqr 

Among the Mira variables, one evolved star 
that illustrates the purpose of this paragraph 
is the complex object, R Aqr. This binary 
star — a Mira (M7e) and a blue companion 
(Ofp, Bep, hot subdwarf, or white dwarf with 
an accretion disk)— is imbedded in a compact 
high-excitation nebula (Lampland, 1923a). 
These three components, which form the spa- 
tially unresolved central object of the R Aqr 
complex, are surrounded by a convex lens- 
shaped nebula of 2 arc-minutes (Lampland, 
1923b) with filamentary and clumpy structures. 
More recently, an inner jet with an apparent 
changing structure and discrete variable knots 
was observed within 10 arc-seconds of the star. 
A possible counter jet was recently pointed out 
with radio and optical mappings. 

Let us briefly review the main studies look- 
ing at the various causes of variability of this 
object and showing its rapid evolution. 

The variability of the classical Mira was dis- 
covered by K. L. Harding in 1811, and the first 
objective-prism spectrograms obtained at Har- 
vard as early as August 1893 show an Me spec- 
trum. One of them, taken on October 17, 1893, 
shows a faint nebular line at X5007 and bright 
hydrogen lines from H/3 to H£, with no trace 
of the Me spectrum (Merrill, 1928). The first 
slit spectrogram, taken at Mount Wilson in Oc- 
tober 1919, shows (superposed on the normal 


spectrum of the M7e star) several characteristic 
emission lines of the gaseous nebula, a com- 
bination never observed previously. 

As pointed out by Merrill (1940), the inner 
nebulosity is variable in position and bright- 
ness. The existence of both inner and outer 
nebulosities have been confirmed by numerous 
photographs taken in 1919 by Hubble (1940, 
1943) with the Mount Wilson 100-inch refrac- 
tor. From the measurements of a pair of plates 
taken at 16-year intervals, Baade (1943, 1944) 
estimates from the expansion rate that the onset 
of expansion was about 600 years ago. These 
observations could confirm the nova-like event 
noted by the ancient Japanese records in A.D. 
930 and which formed the outer extended neb- 
ulosity. 

During many years, the emission lines of the 
inner nebula have been photographed repeated- 
ly. Their intensities vary through a large range, 
apparently without any correlations with the 
phase of the Mira star. Table 20 in Merrill 
(1940) describes the behavior of the spectrum 
of the companion and of the nebular lines from 
1919 to 1939, but Merrill neglected the month- 
to-month variations in his description. 

During the 1920-1930 decades, the very 
high-excitation spectrum marked by bright lines 
of hydrogen, helium, and iron was ascribed to 
the companion, which varies in an irregular 
manner (Merrill, 1935; Campbell, 1938) from 
the 11th magnitude or fainter to the 8th mag- 
nitude in 1933. Slight brightnesses were seen 
during 1922-1925. In July and September 1925, 
the nebular lines and the companion spectrum 
were particularly weak, but by December 1926, 
a rapid and well-marked increase in brightness 
occurred; the continuous spectrum of the com- 
panion was abnormally strong, and the spec- 
tral features were completely different from 
anything previously seen. The bright hydrogen 
lines had P Cygni profiles, and numerous ion- 
ized iron lines were prominent, whereas the 
nebular lines, X4959 and X5007, of [O III] were 
just detectable. 

In 1927, the continuous spectrum and the 
nebular lines were much weaker, but the for- 
bidden iron lines were stronger. Until 1933, no 
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large variation was noticeable; all the lines were 
slightly diffuse but not of great width. During 
1931-1933, Payne-Gapochkin and Boyd (1946) 
observed an apparent correlation between the 
extremely bright active companion and the 
smeared-out light curve of the Mira in which 
the minima were bright and the maxima were 
abnormally low. 

In July 1934, the combined flux of the Me 
and hot companion appeared to be fainter than 
it had been for several years, and the spectrum 
was very different. The [Fe II] lines in the blue- 
violet and the [S II] lines near X4068 were well 
marked, while the Ca II H and K, the Fe II, 
and the nebular lines were very weak. After the 
Mira maximum of August 1934, the Me spec- 
trum became fainter, and in November, the 
post-maximum Mg I X4571 line appeared and 
became rapidly bright. The active companion 
was close to the 10th magnitude during the 
summer of 1935, and its forbidden lines domi- 
nated the spectrum during the Mira minimum 
of July 1935. These forbidden lines remained 
fainter during 1936-1939, few of them were 
recorded superposed to the Me spectrum. 

Some spectral features must be noted for 
understanding this complex object. For 
example: 

1 . The nebular line, X3967, of [Ne III] is not 
present (Wright, 1919); it could be ab- 
sorbed by the Ca II H line and could 
demonstrate the existence of Ca II in the 
line of sight (i.e., in the outer nebulosity 
or in a large absorption shell surround- 
ing it). 

2 . The P Cygni profile of the companion 
lines has often varied in intensity and 
structure, perhaps demonstrating impor- 
tant motions inside the inner nebula. 

Between 1936 and 1949, the Me component 
of R Aqr was varying like a “ normal” Mira 
(Merrill, 1950), while the spectrum of the ac- 
tive companion was subject to changes from an 
Ofp star spectrum with bright lines of N III and 
Fe II to the spectrum of a Bep star with nu- 


merous Fe II lines and a strong continuous 
spectrum. Its nebular lines, weak in 1936, in- 
creased in intensity until 1940 and remained 
relatively strong during a decade. On 10 A/mm 
spectra, their widths appear variable and some- 
times slightly unsymmetrical. Variable motions 
of matter are thus demonstrated in the forming 
layers of these lines. Merrill (1950) suggests 
either a probable orbital interpretation (a spec- 
troscopic binary with 27-year period) or large- 
scale pulsations in the inner nebulosity to inter- 
pret the radial-velocity displacement of the [O 
III] and [Ne III] nebular lines. However, these 
variations reflect a continuing activity in the in- 
ner nebula. 

The spectra obtained by Herbig (1965) in the 
late 1950’s show a very strong nebular spectrum 
that nearly disappears in the observations of Il- 
ovaisky and Spinrad (1966) made during the 
1964 and 1965 minima; present were only the 
[S II], X4068, and the doublet, [O II] XX 
3726-3729, from which an electron density N 
« 10 3 /cm 3 was deduced. However, the Mira 
spectrum at minimum is clearly present, and 
moreover, there were no traces of the active 
companion. 

The IR emission of R Aqr was detected in- 
dependently by Woolf (1969) and Stein et al. 
(1969). The U-pan emission of R Aqr is not 
brighter than that of o Cet or p, Cep, but it is 
much larger (from 8 to 12 pan). It is not as 
broad in wavelength as that observed in the 
planetary nebula, NGC 7027 (Gillett et al., 
1968). This peak shape should be correlated 
with the various sizes of the emission grains 
located in the nebula. Stein et al. (1969) also 
note a strong and narrow absorption feature at 
7.8 pm. Today, the l.S-pm feature found in HD 
44179 is attributed to coronene (C 24 H J2 ) by 
L£ger and Puget (1984). 

The IUE spectra show strong emission lines, 
which are probably formed in the dense com- 
pact nebula located around the binary system 
and close to it (d » 2 x 10 14 cm). The spec- 
tral type of the active companion is not conclu- 
sively resolved; it was identified as a bright 
white dwarf (T > 50000 K) or a subluminous 
central planetary-nebula star. Whatever it is, it 
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would have to be less than T „ * 65000 K in 
order to explain the weakness or absence of the 
He II 1640 A line as stressed by Michalitsianos 
et al. (1980). Its brightness is comparable to that 
of the solar one. As suggested by Merrill (1950), 
it can produce enough ionizing photons to ex- 
cite the continuum and the emission lines, but 
it is relatively faint for direct observations: 
“Such a star would photoionize the inner, high 
density nebula, but lacks a sufficient flux to 
photoionize the entire extended nebulosity.” 
The observed IUE line fluxes, such as those 
from He II, C II, C II], C III], C IV, [O II], 
O III, [O III], O IV, S II, Si III], etc., can be 
used to obtain the general parameters of the 
ionized nebula (T e ~ 15000 K, L « 2 X 10 14 
cm, N = 10 6 to 10 7 cm‘ 3 ). The observed UV 
continuum, essentially flat, is attributed to 
Balmer recombination. Michalitsianos et al. 
(1980) conclude that the compact nebula (in- 
ner nebula) could be entirely due to a mass-loss 
phenomenon from the primary M7e star with 
a mass-loss rate larger than 10' 7 M© /year. To 
explain the 1933 brightness of the companion 
( m v = 8 mag), they propose that it was pro- 
voked by mass transfer from the primary to the 
secondary. A set of parameters for the compan- 
ion was deduced in their Table 3. Moreover, 
from the emission lines such as O I, Mg II, and 
Si II, they deduce the presence of a warm 
chromosphere around the M primary star. 

A model of such a transfer of mass was de- 
veloped by Kafatos and Michalitsianos (1982). 
In this model, the high brightness could be trig- 
gered when the companion accretes matter 
from the Mira component when it crosses over 
the periastron of its highly elliptical orbit. This 
hypothesis was one of the two previously sug- 
gested by Wallerstein and Greenstein (1980) to 
explain the apparent correlation observed by 
Payne-Gaposchkin and Boyd (1946). As sug- 
gested by Usher (quoted by Payne-Gaposchkin 
and Boyd), the companion’s outburst might in- 
hibit the pulsation and suppress the maxima by 
a change in the boundary conditions of the 
Mira itself. However, Wallerstein and Green- 
stein (1980) prefer the possibility of a single 
Mira star with a flaring region made by com- 


plex magnetic fields. During such flares, the 
pulsation of the Mira is inhibited by the 
magnetic fields and its maxima are smeared 
out. Jacobsen and Wallerstein (1975) arrive at 
the same conclusion by radial-velocity mea- 
surements of the nebular lines on plates from 
1957-1965 and from 1970-1971. They do not 
find the 26.7-year period of Merrill (1940), and 
they conclude that a type of activity takes place 
in the outer layers of the Mira during the in- 
terval 1930-1940. 

Using the criterion for Roche formation, the 
value of the primary radius, /? * 2 X 10 13 
cm, and the half-value of the orbit axis, a = 
2.5 x 10 14 cm, Kafatos and Michalitsianos 
(1982) found 0.84 < e< 0.92. With e = 0.85, 
they obtain M ] + M 2 = 2.5 M© , P ~ 44 
years, a visible brightness of the disk of the 
order of = < 1/5 the Mira brightness, and a disk 
optically thick with an external temperature of 
= 2300 K and R . = 2 /?,. If these values are 
correct, it is probable that the gravitational field 
of the secondary could be more important than 
that of the Mira. Kafatos and Michalitsianos 
(1982) believe that the disk around the unseen 
companion is not as large as 2 R v but has its 
larger extent at periastron because a moderate 
dimension of the disk would be more consis- 
tent with the expected mass of the secondary 
(1 M 0 ). 

Viewing the visual light curves of the Mira, 
Willson et al. (1981) suggest that the 1934 and 
1978 low maxima are caused by an eclipse of 
the Mira itself by the accretion disk or by the 
gas cloud of the system and conclude that the 
orbital period is 44 years. 

The thermal relaxation time of the envelope 
of the Mira, as well as the free-fall time scale, 
is about some years. These times are compar- 
able in range to the 8.5-year eclipse duration 
observed between 1928 and 1935 and to the ap- 
pearance of the jet between 1970 and 1977, 
which occurs about 44 years later (i.e., binary 
period). Such observations may not represent 
an eclipse at the periastron ( = 1 year), but they 
would in fact infer the characteristic dynamical 
time scale required by the Mira to recover to 
its preperiastron quiescent state. 
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The first hypothesis of Wallerstein and 
Greenstein (1980), quoted previously, is sup- 
ported by Bath (1977), who proposes that the 
mass transfer from the Mira envelope to the 
companion at its periastron orbital phase builds 
a spatially thick accretion disk at supercritical 
accretion rates, accompanied by the formation 
of a jet, possibly driven by radiation pressure. 
To confirm this hypothesis, Herbig (1980) 
found such a spike on the inner nebulosity, ex- 


tending 10 arc-seconds from the star, with an 
angle position of approximately 22 ° on direct 
plates obtained with the Lick telescope. Tapia 
et al. (1982) reveal an 8" jet at position angle 
(P.A.) 26° formed by discrete variable knots 
and perhaps a diffuse component and ap- 
parently a very small extension at the opposite 
direction not observed before. They suggest 
that this violent ejection may have taken place 
during the current minimum of the variable. 
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Figure 1-29. A 6-cm map of the R Aqr complex showing the secondary source located 
at 196 arc-seconds and the neighborhood of the star (from Kafatos et al., 1983). 
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This optical jet was rapidly confirmed by a 
radio counterpart. Sopka et al. (1982) describe 
its observational properties from optical and 
very large array (VLA) mappings. They found 
extensions both north and south of the star, at 
right angles to the outer nebular arcs, and they 
confirm that the inner nebulosity is variable in 
both brightness and structure, as already men- 
tioned by Hubble (1940, 1943). These observa- 
tions must be made again. 

The mean radial velocity of the jet is -70 
km/s, whereas the mean radial velocity cen- 
tered on the star is -44 km/s (Sopka et al., 
1982). Kafatos et al. (1983) made VLA obser- 
vations of the complex R Aqr, with a resolu- 
tion of 1 arc-second. At 6.4" from the radio 
emission of R Aqr itself, they found a peak 
radio intensity at 29 °3 P.A. (jet B on Figure 
1-29). An unresolved radio source is also found 
at =196 arc-second away on the line defined 
by the previous jet and the star itself. This fea- 
ture may represent matter previously ejected by 
the complex. Besides these two sharp radio fea- 
tures, a new undetached one has been detected 
at =45° P.A. and =2.7 arc-second from the 
central star. 

The near-UV map (Figure 1-30) taken by 
Mauron et al. (1985) duplicates the 6-cm VLA 
map made by Kafatos et al. (1983), but the 
northern knots seem a little more resolved; on 
this near-UV map, the counterjet located in the 
opposite side of the jet is well identified (C in 
Figure 1-30). 

The optical and the radio structures of knot 
B are quite similar, suggesting that we see the 
same emitting region which produces the Bal- 
mer continuum in the near-UV and the free-free 
radiation at 6 cm. 

Knot A seems more radially elongated on 
the near-UV electronic camera photograph than 
in the Kafatos et al. (1983) 6-cm map. More- 
over, near knot A, on the isophotes of the un- 
resolved inner nebula, a bump is visible (E on 
Figure 2 in Mauron et al., 1985). 

Assuming that the R Aqr distance is about 
300 pc (Whitelock et al., 1983), Kafatos et al. 
(1983) deduce from their VLA observations 
that knot A is located about 1.2 x 10 16 cm 


and knot B is about 2.9 X 10 16 cm from the 
central cool Mira variable. Since the knot B 
structure has not changed much in approx- 
imately 1 year (from September 1982 to De- 
cember 1983), if we compare the 6-cm map 
from Kafatos et al. (1983) to the near-UV direct 
electronic camera photograph from Mauron et 
al. (1985), the latter authors deduce that this 
knot cannot have been ejected later than about 
20 years ago. However, the proper motions of 
the knots obtained by Kafatos et al. (1983) and 
Mauron et al. (1985) have very poor accuracy; 
new measurements, chiefly from space, are re- 
quired. 



Figure 1-30. The R Aqr inner nebulosity with 
knots A and B and the possible counter-jet 
feature (from Mauron et al., 1985). 
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Taking into account that Si II XI 81 5 and Si 
III XI 892 are seen in the spectrum of the cen- 
tral region and are absent in the jet spectrum, 
Michalitsianos and Kafatos (1982) suggest that 
the material in the jet is comparatively cooler 
or less dense than in the central region. Quoting 
the optical spectra of the O II and O III emis- 
sion lines of R. Fersen, Kafatos et al. (1983) 
confirm the cooler temperature of the jet. 

In conclusion, these knots may be formed 
in different ways: 

1 . In the Mira inhomogeneous stellar wind, 
being perhaps excited by the UV flux 
from the orbiting companion when they 
enter its Stromgren sphere (Spergel et al., 
1983), with the pressure radiation on dust 
in the knots making them move outward. 

2 . During the periastron phase, when the 
Mira fills in its Roche lobe (accretion at 
supercritical rates; Ferland et al., 1982). 

3 . In a probable polar ejection from the 
disk, giving them a roughly axial orien- 
tation relative to the elliptical outer 
nebula. 

Kafatos et al. (1983) write that: “if the jet 
and the nearby point source (feature A) are 
ejecta, their difference in position angles could 
be interpreted as precession of the system while 
it expels material. The morphology of the outer 
*2' nebula could be accounted for in this man- 
ner, thus explaining the characteristic lens- 
shaped filamentary structure.” This is to be 
confirmed by repeated high spatial resolution 
observations. 

The extended nebula has a mass about 0.2 
Mq , a kinetic energy larger than 2 x 10 46 
ergs. It radiates more than 5 x 10 44 ergs/year 
in Balmer and Lyman emission lines and con- 
tinua. Its cooling time scale is estimated to be 
*=2 years. It seems unlikely that the photo- 
ionization from the companion star is the 
power source of the extended nebula for the 
past 600 years because its excitation requires 
much more energy than the Eddington limit for 
a 1 Mq star. However, from purely energetic 


arguments, the jet itself could be the source of 
the excitation of the outer nebula. The jet 
parameters estimated by Sopka et al. (1982) are 
used by Kafatos and Michalitsianos (1982) to 
conclude to a jet kinetic energy of ^5 x 10 45 
ergs, which the authors estimate in reasonable 
agreement with the radiation cooling given pre- 
viously. 

The material of the extended nebula is 
trapped in rings, generally perpendicular to the 
axis of the jet. Such a structure is favored for 
the formation of the nebula in a single explosive 
event hundreds of years ago. This is confirmed 
by the high content of the nebula in He, and 
more particularly in N (Wallerstein and Green- 
stein, 1980), which is also in agreement with 
models of recurrent novae (Starrfield et al., 
1982). 

The counterjet feature situated on the op- 
posite side of knot A (C in Figure 1-30) was in- 
itially detected by Tapia et al. (1982) on B and 
V plates taken on September 23-24, 1982. By 
radio observations, Kafatos et al. (1983) de- 
tected a possible counterjet. Using the ratio of 
the UV C IV to C III] lines at various places 
on the R Aqr complex, Michalitsianos and Ka- 
fatos (1982) found that this ratio at the opposite 
side of the jet is much different from that 
elsewhere. 

What future projects are now planned for 
clarifying the R Aqr complex? To study the 
structure and the evolution of the inner nebula, 
the knots of the jet, and the counterjet, it is 
necessary to reach Mira phases and wavelengths 
when and where the Mira brightness is very low. 
Observations during the Mira minima and in 
the UV region could be powerful. 

In the aim to diminish or to eliminate the 
turbulence of the Earth atmosphere usually at- 
tached to ground-based observations, observa- 
tions in good seeing observatories by Speckle 
interferometry and Space Telescope (ST) in- 
vestigations are needed. The Faint Object 
Camera (FOC) was developed especially for 
high spatial resolution studies. Moreover, the 
image restoration techniques could help in the 
analysis of observations which are at the limit 
of resolution of the ST, such as for the inner 
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nebula of R Aqr. The proper motions and evo- 
lution of the shape of the elements of the inner 
and outer nebulas could be investigated by this 
instrumentation. The kinematics of various fea- 
tures of the R Aqr complex could be deduced. 
High-resolution spectra of the jet and the knots 
should be made to evaluate the richness of the 
extended nebula in helium and nitrogen — rich- 
ness suspected by Waller stein and Greenstein 
(1980). 

Irreversible Changes in SR and L Stars 

Many M and C class stars showed outstand- 
ing and nonreversing variations in one or more 
decades. We shall elaborate here on some ex- 
amples of these rapid and definitive changes. 

V1016 Cyg (MHa 328-1 16) was classified as 
a strong Ha emission object (Merrill and Bur- 
well, 1950) with an estimated V magnitude 
around 15 (Fitzgerald et al., 1966). Between 
July 1963 and August 1965, this object bright- 
ened by 4 mag (McCuskey, 1965). Since this 
outburst, its behavior has been monitored op- 
tically by Mammano and Ciatti (1975); in the 
infrared by Harvey (1974), Puetter et al. (1978), 
and Aitken et al. (1980); in the radio by Pur- 
ton et al. (1973) and Altenhoff et al. (1976); and 
in the UV-IUE by Flower et al. (1979) and 
Carpenter and Wing (1979). In the visible, 
many lines indicate a very hot source of radia- 
tion although some low-excitation lines are also 
present, and this, together with the absorption 
bands of TiO and VO, has been taken as evi- 
dence of a binary system with an LPV and a 
hot star exciting a nebula (Mammano and Ciat- 
ti, 1975). From their analysis of the IUE spec- 
trum, Carpenter and Wing (1979) conclude that 
there has been sudden shell ejection with a very 
large excitation which developed a rich 
emission-line spectrum, with lines of O I, Fe 
II, Mg II, C II, N V, O V, [Ne V], and [Mg 
V]. As these lines are also seen in the spectra 
of planetary nebulas such as NGC 7027, 
something like a planetary nebula was probably 
formed or excited around the star (see also 
Flower et al., 1979). Nussbaumer and Schild 


(1981) interpret that spectrum with a single star 
planetary model, the central star having a 
temperature of T = 160000 K and a radius of 
0.06 R © . The expansion velocity deduced from 
the line shape is around 100 km/s. The elec- 
tron density in the shell is approximately 3 x 
10 6 cm" 3 , and the electron temperature varies 
from 8000 to 18000 K. From the reddening 
determined in the IUE observed UV region, 
they derive a distance of 2.2 kpc. 

The carbon star, HD 59643, provides an- 
other example of nonreversing and rapid evolu- 
tion. Shane (1928) and Keenan and Morgan 
(1941) considered it as a cool, but otherwise 
rather ordinary, nonvariable carbon star and 
classified it as R8 or R9. After reduction of its 
spectrum, Wildt (1941) considered it as a nor- 
mal carbon star, with normal excess of heavy 
elements and with Ca II H and K lines in ab- 
sorption. After more than 25 years, Gordon 
(1967) and Yamashita (1967) noticed an in- 
crease in the equivalent widths of the absorp- 
tion lines of the heavy elements Ba, La, Sc, and 
Sr, as well as a deeper molecular G band. 
Greene and Wing (1971) examined the plates 
of HD 59643 taken by Keenan in 1949 and 
found no evidence of hydrogen emission. Ha 
was seen in emission for the first time in No- 
vember 1966 by Gordon (1967), when it was a 
little brighter than its nearby continuum. In 
1969, Utsumi (1970) saw H/3 quite clearly in ab- 
sorption. These data suggest that the activity 
is of recent origin but not easy to understand. 
Greene and Wing (1971) found an abnormal 
and variable spectrum below 4000 A with a fill- 
ing in of the K line and the (0,0) CN violet 
band. This activity of February 1970 was also 
observed near the hydrogen emission lines 
H6-H£, but the long ward spectrum was not 
significantly different from that analyzed by 
Wildt (1941). The main difference is that the 
Ca II H and K lines were completely invisible 
in 1970, although they were in absorption in 
1941. This absence cannot be interpreted as a 
calcium deficiency because the Ca I X4226 line 
appeared to be strong. Another difference is the 
appearance of a very weak Ha emission line in 
1966, along with the strong H5-H£ and the UV 
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Figure 7-57. Low-resolution IUE spectra of HD 59643 . The abscissas are in 
Angstroms, the ordinates are absolute fluxes . The lower spectra are obtained from 
the data bank . The observations are made on: Feb . 17, 1979 (SWP 4290); Feb . 14, 
1979 (LWR 3765); Oct . 25, 7957 75550, LIFT? 11846); March 23, 1982 (SWP 

16603, L WR 12845); April 27, 1983 (SWP 19829, L WP 1853). We note the high level 
of noise from 2000 to 2200 A for the L WP spectrum and a spike at 2180 A in the 
L WR spectra (from Querci and Querci, 1985b). 


continuous emission in 1970. Greene and Wing 
(1971) conclude that the UV continuous emis- 
sion does not come from the photosphere; the 
radiation must emanate either from a hot cir- 
cumstellar envelope or from the secondary star 
of the binary system. The latter interpretation 
disagrees with Wildt’s observations, in which 
no continuum or emission lines are reported. 
Therefore, the UV emission could be con- 
sidered to be an activity located above the layers 
responsible for the violet Ca II H and K lines, 
giving an IUE spectrum (Carpenter and Wing, 
1979) with ions of a great range in ionization: 


O I, Mg II, C II, and N IV. Querci and Querci 
(1985b) demonstrate that the emission lines are 
variable and that the continuum from 2200 to 
1200 A decreased significantly between 1979 
and 1983 (Figure 1-31). The Balmer lines in HD 
59643 are formed in the same outer layers be- 
cause the Balmer decrement is very small, 
whereas it is very steep in the Mira variables 
and consequently affected by the violet de- 
pression. Through IUE and visible high- 
resolution spectra, the Bloomington and 
Toulouse groups are following the evolution of 
this complex object. 


i 
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The ultraviolet continuum of the SRb car- 
bon star, UU Aur, decreases with time. Shajn 
and Struve (1947) measured the Ca II H and 
K lines to be seven times weaker than those of 
the gMO comparison star, HD 80943, while the 
Ca I 4226 A line is not weakened and the Sr 
II 4078 A line seems to be less affected. They 
suggest that the weakening of the continuous 
violet background may be attributed to the ef- 
fect of many overlapping unknown molecular 
bands rather than to continuous absorption, in 
any case not to some temperature effect. They 
also measured the wavelength of emission lines 
like X\ 4036, 4009, 3959, 3949, 3914, and 
others, whereas in the regions centered at XX 
3983, 4005, and 4052, there is no apparent emis- 
sion in the N stars. In the 200-inch plates ob- 
tained by Bowen (1951) in January 1951, the 
H and K lines, as well as the A1 I 3944 and 
3961.5 A lines, are not visible on a spectrum 
which has a very low continuous background 
without absorption features shortward to Mn 
I 4030 A multiplet. Only numerous emission 
lines are observed, especially at 3982.5 A, which 
are identified as Ti I, V I, or Zr I low-excitation 
lines (Gilra, 1976). We therefore suggest an in- 
crease of the optical depth of the absorbing 
layers situated at the temperature minimum, 
which smears out all the absorption features 
coming from layers situated beneath. The very 
extended layers situated above might give rise 
to emission lines (low-temperature chromo- 
sphere) and narrow absorption lines (extended 
shell) such as those from Mn I multiplets, which 
eat away the emission lines. It is suspected that 
the changes in UU Aur are more in the nature 
of variations than secular changes. 


Irreversible Changes in HdC Stars 

The irreversible evolution for the different 
RCB stars is not exactly the same. In the fol- 
lowing, we detail the observed evolutions of 
only a few chosen RCB stars, all the more so 
because, unfortunately, the sample of the cor- 
rectly observed RCB is very poor. 


A Slowly Evolving Star: RY Sgr. Using data 
from 1967 to 1970, Alexander et al. (1972) con- 
firm the cepheid-like behavior of RY Sgr by 
spectroscopic and photoelectric observations. 
Yet again, they find the cyclic variations of 0.6- 
mag amplitude and the 38.6-day period found 
by Jacchia (1933), with data covering the period 
1920-1932. They conclude that “neither evolu- 
tionary effects nor mass loss by ejection are of 
enough importance to significantly change the 
period on the time scale of 50 years.” However, 
observations with longer time intervals could 
give slightly different results. Pugach (1977) 
was the first to point out a shortening of 0.9 
days in 40 years, using observations from 1926 
to 1977, but with a gap of 168 pulsational 
cycles. However, he was unable to find a rela- 
tion between the period shortening and the 
time. This decrease of the pulsational period 
length was confirmed by the results of Marraco 
and Milesi (1980, 1982) and those of Kilkenny 
(1982). With observations of RY Sgr from 1897 
to 1977, Marraco and Milesi (1982) give the in- 
stantaneous periods of the pulsational oscilla- 
tions in 1897, 1926, and 1977 to be 39.3, 39.0, 
and 38.2 days, respectively. The periods of Jac- 
chia (1933) and Pugach (1977) are again found, 
and the period decrease of 1.1 day in 80 years 
is deduced. Kilkenny (1982) reexamines the ar- 
chive material from 1926 to 1978, material 
largely due to the efforts of amateur 
astronomers. He shows that the epochs of the 
observed minima agree with JD min = T 0 +P Q n 
+ kn 1 , where T Q , P Q are the “zero” epoch and 
period, respectively, and k represents the linear 
rate of change of the pulsational period. The 
decreasing period thus ties in with the general 
data rather well, but there appears to be some 
modulation of k , with a time scale of 100 or 
200 pulsation periods. Therefore, RY Sgr has 
k = -0.0005 (Kilkenny, 1982), a value which 
is close to that expected by Schonberner (1977) 
for a deficient hydrogen star of M = IM © and 
r ff = 6900 K, evolving rapidly from the red- 
giant to the white-dwarf stage. 

A Star Rapidly Changing Its Period: S Aps. 

Waters (1966) analyzed the visual observations 
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Figure 1-32. Visual observations (3-day means) 
of S Aps. Abscissas are Julian dates . Symbols 
indicate the number of observations in the 
mean: squares have five or more observations; 
filled circles three or four; and open circles , 
two; singles observations are omitted (from 
Kilkenny , 1983). 


made during the period 1922-1960 and found 
a period of about 120 days. Kilkenny (1983) in- 
vestigates 8000 individual observations cover- 
ing the interval 1960-1982. After the obscura- 
tional minimum of 1967, S Aps shows an 
instantaneous pulsational period of about 135 
days (Figure l-32a). Surprisingly, after the next 
obscurational minimum in 1971, S Aps appears 
to have developed an oscillation of a 37.5-day 
mean pulsational period (Figure l-32b). Obser- 
vations on the 1979 maximum show an instan- 
taneous period of about 40 days (Figure l-32c). 
Other observations of such a phenomenon are 
necessary before we can conclude that the 
period change is linked to the obscurational 
minima or if the deep minima facilitate the 
change of pulsational mode (from fundamen- 


tal to first or other overtones). We cannot ex- 
plain this observation in terms of decreasing 
pulsational rate, and there is no model of evolu- 
tion of hydrogen-deficient stars which allows 
such a rapid change of period (Roser, 1975), 
although the evolution of a hydrogen-deficient 
star to a white dwarf is rapid. Therefore, any 
changes in the observational material should be 
virtually undetectable on a time scale of a few 
years (Schonberner, 1977). In addition to these 
changes on pulsational period, there appears to 
be a fading of the amplitude of the pulsational 
variations in some RCB stars (Fernie et al., 
1972; Kilkenny, 1983). 

As for UW Cen, another example, the situa- 
tion is not as comprehensible as for S Aps 
because the reduction of the observations gives 
many solutions with +0.003 -0.006; con- 

sequently, the period could be constant or 
decreasing/increasing, but at a relatively low 
rate. We must point out that “if UW Cen has 
an increasing period, this could mean that it has 
not yet reached the top of the asymptotic giant 
branch” (Kilkenny and Flanagan, 1983). 

The helium stars could also have irrevers- 
ible changes. For example, HD 168476 had V 
= 9.37 mag in 1964 (Hill, 1964); the data of 
1969-1971 given by Landolt (1973) and com- 
municated to him by Hill in 1973, as well as the 
observations of mid- 1972 made by Landolt 
(1973), seem to indicate a slight secular bright- 
ening over the decade 1964-1973. 

The V mean value and the colors of the HdC 
star, HD 182040, obtained by Rao (1980b), 
agree with the data deduced from the one 
observation of this star made 25 years before 
by Mendoza and Johnson (1965). Consequent- 
ly, no middle-term variations seem to appear. 


PROSPECTS AND CONCLUSIONS 

At the end of this description of the slow 
or rapid, reversible or irreversible phenomena 
presented by the cool giants and supergiants, 
we should like to emphasize the observations 
to be performed for having a necessary solid 
basis of data to model these hydrodynamical 
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phenomena, and progress in knowledge of the 
evolutionary status of these stars. 

To evaluate the nonthermal energy varia- 
tions on various time scales involved in the dif- 
ferent layers of the entire stellar atmosphere, 
simultaneous observations in many spectral 
ranges are required. Consequently, extensive 
monitoring must be organized from the ground 
with the present and new well-adapted in- 
struments and with space facilities that will 
become more and more available in the not too 
distant future. 

Some of the ground telescopes and their 
focal instrumentation could be used for photo- 
metric, high-resolution spectroscopic, and 
polarimetric analysis, as well as in the visible 
and infrared ranges and in the radio range. 

To stimulate long-term variability observa- 
tions of the carbon-star energy distribution, to- 
day partially followed with non-well-adapted 
means, Querci and Querci (1985b) propose a 
photometric system covering the main atomic 
and molecular features over the largest wave- 
length range. This system is inspired by their 
previous studies of N-type spectra from the far- 
UV to the radio frequencies and by synthetic 
spectra deduced from their model-atmosphere 
calculations. Also selected are three points on 
the continuum: at 1 .5 mm, 3.8 fim, and 1 .04 /*m 
(Wing’s 1104 filter). Active observations in the 
visible are presently made in ESO-Chile by the 
Sterken group with some of the Querci’s filters. 
Variations in color indices such as in C 2 , CN, 
D-Na, Ha, 1104, etc., can be followed, and 
various correlations with the known physical 
parameters can be shown. Monitoring on other 
spectral intervals is urgently required, such as 
around C IV (if any, possibly at some phases) 
and CO features and at 3.8 jim, 11.3 jim, and 
1.5 mm, etc. This should allow a detailed anal- 
ysis of the behavior of the various stellar layers. 
In addition, the high spectral resolution pres- 
ently available, which gives fine details on the 
structure and dynamics of the red giants (see 
M. Querci, this volume), should let us envisage 
the modeling of the entire variable atmosphere. 

Collaboration with amateurs on photomet- 
ric monitoring cannot be neglected. The high 


instrumental and scientific levels of some of 
them permit productive cooperation. The In- 
ternational Amateur-Professional Photometric 
Photometry (IAPPP) members, mainly in the 
United States, and the Groupe d’Etudes et 
d’Observations Stellaires (GEOS), mainly in 
Europe, are experienced in the photoelectric 
photometry. In the United States, after the first 
telescopes and photometers that work on semi- 
automatic control or remote control, some 
Automatic Photometric Telescopes (APT’s) of 
40-cm aperture are already operating, mainly 
at the Fairborn Observatory in Phoenix, Ari- 
zona (Boyd et al., 1985). These telescopes are 
driven by a computer that has on memory the 
star program and the comparison stars. The 
tests on sky quality and the night observations 
are automatically made on-line (i.e., finding the 
stars and centering them on the aperture, re- 
cording the data, etc.). The computer makes 
decisions such as when to open the dome, in 
what order the stars should be observed, and 
when to shut down. The telescope calls for 
human assistance if some failures appear. The 
data are reduced at the end of the night. The 
accuracy of such measurements made with the 
UBV filters of the Johnson system is of -0.02 
magnitude. An APT service permits anyone to 
request differential photometry for his own 
program stars. 

Many efforts are also being made in Europe; 
among them, the twin-telescope technique is be- 
ing developed by the Section de Photometrie 
Photoelectrique of GEOS (Gregory and Quer- 
ci, 1985). Because the sky transparency is rapid- 
ly variable during a night in Western Europe 
and because the GEOS amateurs aim at an ac- 
curacy of 0.002 mag in UBV filters, human 
assistance during the night is preferred for 
various technical and sky-quality controls 
through the rms of the on-line data on the com- 
puter display. In the twin-telescope technique, 
two telescopes with 28-cm apertures are fixed 
on the same mounting: one follows the star pro- 
gram, and the other, which is also driven by 
computer, follows the comparison star. The 
angle between their optic axes is calculated by 
the computer, which drives stepper motors for 


89 


I 


k 

h 


a differential orientation of the second quoted 
telescope relative to the first. The two stellar 
lights are led into the two-channel photometer 
by liquid optic fibers. This new photometer is 
the second European Photoelectric Photometer 
(EPP-2) built by Fontaine (1985). It was derived 
from the EPP-1 (one-channel photometer) 
made by Fontaine (1984) and Walker (1985) in 
a collaboration decided at the IAPPP meeting 
at Herstmonceux-Castle during the fall of 1984. 
Full-sky photometry is made by using the stan- 
dard stars. The UBV filters are now being used. 
Soon, we plan to use narrower filters and larger 
telescopes. 

This photometric monitoring is the basis of 
research on variability. Narrow-filter photo- 
metry with 1-m telescopes and visible and in- 
frared high-resolution spectra with 2.0- or 
3.6-m telescopes could be made when interest- 
ing phases of variation are detected by photo- 
metric monitoring with small telescopes. Larger 
telescopes, such as VLTs and NNTTs, could 
be required when the stellar flux is very low. 

It is also worth noting that combined radio 
and infrared observations should be developed 
between radio observatories such as Onsala, 
IRAM, Eiffelsberg, and ESO and infrared tele- 
scopes such as UKIRT, CFHT, and ESO to 
provide constraints on models describing the 
molecular excitation mechanisms on the emit- 
ting layers, since it is well known that the in- 
frared continuum and lines can excite the radio 
molecular transitions (see Nguyen-Q-Rieu, this 
volume). New molecules could be detected in 
the millimetric wavelength. The very far-UV 
and submillimetric wavelength ranges (i.e., be- 
tween 200 and 900 fim) have not been exten- 
sively examined. There, photometry, together 
with spectroscopy, should give insight into the 
outer layers of the cool stars and special iso- 
topic abundance determinations. 

Already, the space gave a new scope on the 
cool giant research through the IUE and IRAS 
spacecrafts. (See the section Photometric Ob- 
servations.) However, the small aperture of 
IUE seriously limits the high-resolution analysis 
to a very small number of cold bright stars. 


The first EXOSAT investigations on cool 
giants and supergiants are presently in progress. 
(See the section Photometric Observations .) It 
seems to us that X-ray fluxes should be obtain- 
ed at some phases of Miras (i.e., when the 
shock waves are the most powerful), and some 
constraints could therefore be added on the 
shock-wave interpretation of the Mira 
variations. 

In the near future, new spacecraft will be 
launched: the Hubble Space Telescope (ST), the 
largest one with its 2.4-m diameter, by NASA 
in 1988, the space astrometry satellite, Hippar- 
cos, by ESA in July 1988, and the Infrared 
Space Observatory (ISO), another ESA space- 
craft, foreseen for 1992. 

The scientific community has received de- 
tailed information on the different focal in- 
struments of the Hubble Telescope via the 
“Space Telescope Science Institute” in Balti- 
more, Maryland (The Space Telescope Obser- 
vatory, 1982). With its planned mission over 15 
years, long-term programs on the variability of 
Miras and semiregular, irregular, and eruptive 
RCB stars are possible. Such an observational 
basis could permit great progress in the knowl- 
edge of the hydrodynamics of these stars. Many 
proposals have already been presented on the 
cool giants and supergiants, some of which are: 

• WF/PC observations for the spatial dis- 
tribution analysis of the clumpy ejected 
matter around individual stars as a func- 
tion of star distance (e.g., a Ori), 

• WF/PC observations of the R Aqr jet 
in the direction of the secondary radio 
source seen on the 6-cm map and analy- 
sis of the jet formation, 

• FOC observations of the R Aqr com- 
plex, mainly the high-excitation nebula 
surrounding the mass-losing Mira 
variable, and the compact hot compan- 
ion with its accretion disk for hydro- 
dynamical modeling of the whole 
object, 
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• FOC observations of extended chro- 
mospheres for an analysis of the spatial 
properties of the wind such as asym- 
metric mass loss, radial variations of 
density and temperature, interaction of 
cool wind with the cool companion 
(e.g., a Ori), or interaction of cool 
wind with a hot companion (e.g., a Sco 
and o Cet), 

• FOC observations of the rapidly evolv- 
ing stars such as V1016 Cyg, HM Sge, 
and HD 59643 to determine the loca- 
tion of the variable high-excitation 
region pointed out by IUE spectra, 

• FOC exploration of the nearby giants 
for the structure of their own surface 
(convective cells, spots, etc.) and for 
the determination of the spatial and 
particle-size distributions on their near- 
by neighborhood (e.g., a Ori and o 
Cet), 

• HRS Echelle spectra of supergiants for 
dynamical analysis of the various 
layers by the monitoring of radial 
velocities and profiles of the spectral 
lines (e.g., 5 Yir, a Ori, and R CrB). 

Hipparcos (High Precession Parallax Col- 
lecting Satellite) was named in memory of Hip- 
parcos, who was the first Greek astronomer to 
measure a parallax (that of the moon) and who 
cataloged 800 stars. During the 2.5 years of the 
Hipparcos mission, any particular point of the 
sky enters in both fields of the satellite several 
times a year, making accurate photometry of 
selected stars possible (Turon-Laccarieu, 1978). 
The absolute coordinates of each star, its par- 
allax, and each component of its annual prop- 
er motion will be measured with a maximum 
mean error of ±0.002 arc-seconds (for stars 
brighter than B = 12). Many astrophysical 
parameters can be deduced from these astro- 
nomical data: (a) with the parallax, we get the 
distance of the star; (b) with its apparent mag- 
nitude and distance, we obtain the absolute 


magnitude; (c) from already measured apparent 
angular diameter and distance, we deduce the 
stellar radius; and (d) with the effective tem- 
perature and stellar radius, we have the lumi- 
nosity and the mass of the stars. During the 
Hipparcos mission, the period of many Miras 
could be defined, the calibration of the period- 
luminosity relation could be improved, and the 
theory of pulsation of Miras could be checked. 
Knowing the center-of-mass radial velocity (ob- 
tained from the radio thermal lines) and the 
proper motions, the spatial motion of the var- 
ious populations of stars could be deduced. 
Moreover, the satellite will collect a consider- 
able quantity of accurate photometric data on 
the variability of the Miras and the semiregular 
and irregular stars; their analysis could provide 
the various pulsating periods and their temporal 
changes. Tests on stellar evolution and struc- 
ture theories could be developed. The variations 
on various time scales detected on some stars 
could also be analyzed. 

During the 18-month mission of the 60-cm 
cooled telescope named Infrared Space Obser- 
vatory (ISO), some interesting programs on 
late-type giants and supergiants should be 
developed: (a) temperature distribution on cir- 
cumstellar material and planetary dust clouds; 
(b) determination of the shape of various dif- 
fuse features (e.g., solid-state spectral features 
at 1 1 and 18 /un) and correlation of these shapes 
with the physical parameters of the star and its 
envelopes; (c) mapping of the neighborhood (1 0 
by 1 °) of some hypergiants and planetary 
nebulae at 40, 80, 120, and 160 fim and deter- 
mination of the shell diameters, energy balance, 
densities, etc.; and (d) determination of mass- 
loss rates. The problem described by Wesselius 
(1984) on the origin of the outer envelopes of 
a Ori should be solved by the ISO high- 
resolution spectroscopy (Barlow and Storey, 
1984). 

Finally, a large organized international col- 
laboration on photometry, spectroscopy, and 
polarimetry both from the ground in the visi- 
ble, near-IR, and radio ranges and from space 
in the X-ray, far-UV, near- and far-IR, and 
submillimetric frequencies seems unavoidable 
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if we want to effectively progress in the knowl- 
edge of these variable cool stars. 

Constructive comments by Hollis R. 
Johnson have improved the first part of this 
chapter. 
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